Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,440,173 T→C E286G (GAG→GGG)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,440,1730TC97.3% 375.4 / ‑2.5 111E286G (GAG→GGG) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (46/62);  total (48/63)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 111 total aligned reads displayed.

AACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACG  >  NC_000913_3_pae_pgi/3440026‑3440318
                                                                                                                                                   |                                                                                                                                                 
aaCCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCg                                                                                                                                                 <  1:1236410/150‑1 (MQ=255)
aaCCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCg                                                                                                                                                 <  2:146504/150‑1 (MQ=255)
   cTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTa                                                                                                                                              >  1:1057434/1‑150 (MQ=255)
    tgtgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTAc                                                                                                                                             >  2:10787/1‑150 (MQ=255)
    tgtgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTAc                                                                                                                                             >  2:787094/1‑150 (MQ=255)
    tgtgATCCGGTTACTCGTCAGCGATGCTCGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGGAAGAGATTTTTTACCAAGGCTAGGCGTTTTAAGGAGCACAACCCCGGTAc                                                                                                                                             >  4:334531/1‑150 (MQ=255)
      tgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCCGTCGCATGCCCAGAGGCAGTCCACGGGAAGCCGGCACGTCTTTAATCCCGGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGc                                                                                                                                           >  2:1157765/1‑150 (MQ=255)
      tgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGTGTTTTAAGGAGCTCAACCCCGGTACGc                                                                                                                                           <  1:918250/150‑1 (MQ=255)
       gATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCt                                                                                                                                          <  4:557430/150‑1 (MQ=255)
        aTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTg                                                                                                                                         >  1:303129/1‑150 (MQ=255)
        aTCCGGTTACTCGTCAGCGAGGCGTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTg                                                                                                                                         <  1:10787/150‑1 (MQ=255)
           ttggcccTCGTCAGCGATGCTTGCCGGTGGACAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTAc                                                                                                                                      <  1:827316/144‑1 (MQ=255)
           cGGTTACTCGTAAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGTAGCCAGCTCGTCTTTTATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTAc                                                                                                                                      <  2:1269835/150‑1 (MQ=255)
             gTTACTCGTCAGCGATGCTTGTCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCa                                                                                                                                    >  2:208099/1‑150 (MQ=255)
              ttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAg                                                                                                                                   >  2:529859/1‑150 (MQ=255)
                  tCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCa                                                                                                                               >  1:143938/1‑150 (MQ=255)
                  tCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCa                                                                                                                               >  1:764047/1‑150 (MQ=255)
                     tCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCg                                                                                                                            <  3:534084/150‑1 (MQ=255)
                      cAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGa                                                                                                                           <  1:373483/150‑1 (MQ=255)
                           aTGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAg                                                                                                                      <  1:1107224/150‑1 (MQ=255)
                            tGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGt                                                                                                                     <  4:135314/150‑1 (MQ=255)
                             gCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGGGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTg                                                                                                                    <  2:425827/150‑1 (MQ=255)
                              cTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTgg                                                                                                                   <  2:511350/150‑1 (MQ=255)
                               ttGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGa                                                                                                                  <  1:1183334/150‑1 (MQ=255)
                               ttGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGa                                                                                                                  <  2:448409/150‑1 (MQ=255)
                                  ccGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCAATATAGTGGATAg                                                                                                               >  2:1126030/1‑150 (MQ=255)
                                     gTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCtt                                                                                                            >  2:704644/1‑150 (MQ=255)
                                      tGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTATACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTc                                                                                                           <  2:686882/150‑1 (MQ=255)
                                       ggCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCt                                                                                                          <  3:696122/150‑1 (MQ=255)
                                         ccAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                                                                                                        >  4:554907/1‑150 (MQ=255)
                                          cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                                                                                                       <  1:529859/150‑1 (MQ=255)
                                          cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                                                                                                       >  1:426981/1‑150 (MQ=255)
                                          cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                                                                                                       >  2:599696/1‑150 (MQ=255)
                                          cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                                                                                                       >  3:427354/1‑150 (MQ=255)
                                           aGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt                                                                                                      >  4:614497/1‑150 (MQ=255)
                                             ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa                                                                                                    <  1:1133671/150‑1 (MQ=255)
                                                tCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg                                                                                                 <  4:713880/150‑1 (MQ=255)
                                                   aGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg                                                                                              >  1:525142/1‑150 (MQ=255)
                                                    ggCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                                                                                             <  2:847469/150‑1 (MQ=255)
                                                       gcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg                                                                                          <  1:569770/150‑1 (MQ=255)
                                                            cccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag                                                                                     >  4:396296/1‑150 (MQ=255)
                                                             aaagagACAGTCCACGGGAAGCAAGAACGTCATTAATCTCAGTAAGAGATTTTTTACCAAGGTTTGGCGTTTTAAGGAGCTCAACCCAGGTACGCTGTACCAGATTACCGATATAGTGGATAGTTTCTGCTTTAAGGCAGTTAGCAGAGc                                                                                    <  2:881697/148‑1 (MQ=255)
                                                              cAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCg                                                                                   <  1:555314/150‑1 (MQ=255)
                                                                     aGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAggc                                                                            <  4:487876/150‑3 (MQ=255)
                                                                     aGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAggc                                                                            >  4:460659/1‑148 (MQ=255)
                                                                       tCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCaa                                                                          >  2:85504/1‑150 (MQ=255)
                                                                         cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAAtt                                                                        <  3:62707/150‑1 (MQ=255)
                                                                          aCGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAc                                                                                                                              <  3:52411/95‑1 (MQ=255)
                                                                            gggAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCa                                                                     <  3:311235/150‑1 (MQ=255)
                                                                            gggAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCa                                                                     <  1:753290/150‑1 (MQ=255)
                                                                              gAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAATGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGa                                                                   >  4:72682/1‑150 (MQ=255)
                                                                              gAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGa                                                                   >  2:170822/1‑150 (MQ=255)
                                                                              gAAGCCAGCACGTCTTTAATCCCAGTAAGAGATTTTTTACCAAGGTCAGGCGCTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCCTTAAGGCAGTTAGTAGAGCGGACAGGCAATCCCAGa                                                                   >  2:26937/1‑150 (MQ=255)
                                                                                aGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATc                                                                 <  4:442407/150‑1 (MQ=255)
                                                                                aGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATc                                                                 <  1:599696/150‑1 (MQ=255)
                                                                                  ccAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGt                                                               >  1:431196/1‑150 (MQ=255)
                                                                                   cAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTc                                                              <  1:787094/150‑1 (MQ=255)
                                                                                    aGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCa                                                             >  4:248755/1‑150 (MQ=255)
                                                                                      cACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATc                                                                 >  3:740937/1‑144 (MQ=255)
                                                                                      cACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATc                                                                 <  4:740937/144‑1 (MQ=255)
                                                                                         gTCTTTATTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAgg                                                        >  4:673095/1‑150 (MQ=255)
                                                                                           cTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGgc                                                      >  1:1019841/1‑150 (MQ=255)
                                                                                           cTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGgc                                                      <  2:1194992/150‑1 (MQ=255)
                                                                                            tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACAGAtata                                                                                                                       <  1:66644/84‑1 (MQ=255)
                                                                                                 tctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGa                                                                   <  2:86045/131‑1 (MQ=255)
                                                                                                 tctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGa                                                                   >  1:86045/1‑131 (MQ=255)
                                                                                                   tcAGTAAGAGATTTTTTGCCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAgg                                              >  2:1250445/1‑150 (MQ=255)
                                                                                                     aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGt                                                                                                                                               <  1:459946/51‑1 (MQ=255)
                                                                                                     aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGt                                                                                                                                               >  2:459946/1‑51 (MQ=255)
                                                                                                     aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggat                                            <  3:673095/150‑1 (MQ=255)
                                                                                                        aaGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgg                                         <  4:427354/150‑1 (MQ=255)
                                                                                                          gagaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggat                                       <  2:303129/150‑1 (MQ=255)
                                                                                                          gagaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggat                                       >  4:563520/1‑150 (MQ=255)
                                                                                                           agaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                      <  2:525142/150‑1 (MQ=255)
                                                                                                           agaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGGACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                      <  3:417472/150‑1 (MQ=255)
                                                                                                               tttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgga                                        <  2:1138922/144‑1 (MQ=255)
                                                                                                               tttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgga                                        >  1:1138922/1‑144 (MQ=255)
                                                                                                                     ccAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGt                            <  4:594537/150‑1 (MQ=255)
                                                                                                                      cAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGAt                                                                  <  4:128625/111‑1 (MQ=255)
                                                                                                                      cAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGAt                                                                  >  3:128625/1‑111 (MQ=255)
                                                                                                                       aaGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGttt                          >  1:96399/1‑150 (MQ=255)
                                                                                                                       aaGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGttt                          <  2:426981/150‑1 (MQ=255)
                                                                                                                       aaGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGttt                          <  1:234821/150‑1 (MQ=255)
                                                                                                                          gTTAGGCGTTGTAAGGAGCGCGACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCGTCTGCTTTAAGGCAGTTAGCAGAGTGGACAGGCAATTCCAGTTCGTCAACAGGGCGTTGCAGGATCGGATCGAACTCTGGTTtctc                       <  2:110188/150‑1 (MQ=255)
                                                                                                                             aGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttc                    <  2:764047/150‑1 (MQ=255)
                                                                                                                              ggCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGCGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGCTCGGATCGAACTCTGGTTTCtcttct                   >  1:1263207/1‑150 (MQ=255)
                                                                                                                              ggCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATAACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                   <  2:886719/150‑1 (MQ=255)
                                                                                                                               gCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCgca                                                   <  4:103064/117‑1 (MQ=255)
                                                                                                                               gCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCgca                                                   >  3:103064/1‑117 (MQ=255)
                                                                                                                                cGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                 <  3:185119/150‑1 (MQ=255)
                                                                                                                                cGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                 <  3:396296/150‑1 (MQ=255)
                                                                                                                                  ttttAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCa               <  2:31551/150‑1 (MQ=255)
                                                                                                                                   tttAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGCTATAGTGGATAGCTTGTGCTTTAAGGCGGTTAGCAGAGCGGACAGGCCATTCCAGATCGTCAACCGGGCGCAGCAGGATCGGATCGAACTCTGGTTTATCTTCTTTCAc              <  2:568134/150‑1 (MQ=255)
                                                                                                                                   tttAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCAc              >  1:754647/1‑150 (MQ=255)
                                                                                                                                       aGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCa          >  4:182774/1‑150 (MQ=255)
                                                                                                                                       aGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCa          <  1:26937/150‑1 (MQ=255)
                                                                                                                                         gAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAgg        >  4:423833/1‑150 (MQ=255)
                                                                                                                                            cTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTg     >  4:446656/1‑150 (MQ=255)
                                                                                                                                             tCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGa    <  1:909695/150‑1 (MQ=255)
                                                                                                                                               aaCCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACg  >  3:799654/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACG  >  NC_000913_3_pae_pgi/3440026‑3440318

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑