Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,313,765 A→T 5.4% A750A (GCT→GCA infB ← translation initiation factor IF‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,313,7650AT 5.4% 542.1 / 13.4 168A750A (GCT→GCAinfBtranslation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base A (82/76);  new base T (6/3);  total (88/80)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

Only 100 of 169 total aligned reads displayed.

CTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCGGAGATCGCTTCGAC  >  NC_000913/3313617‑3313914
                                                                                                                                                    |                                                                                                                                                     
cTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGAGAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAg                                                                                                                                                      <  2:1734852/150‑1 (MQ=255)
cTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAg                                                                                                                                                      >  1:1366308/1‑150 (MQ=255)
  tCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCa                                                                                                                                                    >  2:1997970/1‑150 (MQ=255)
  tCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCa                                                                                                                                                    >  1:168200/1‑150 (MQ=255)
    cGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCacg                                                                                                                                                  <  2:1896831/150‑1 (MQ=255)
    cGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCacg                                                                                                                                                  <  2:1503925/150‑1 (MQ=255)
       cAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCTGATGCCTCTGCacgtaa                                                                                                                                               <  2:1644893/150‑2 (MQ=255)
         gACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCTGATGCCTCTGCacgtaaag                                                                                                                                             <  2:1954817/150‑4 (MQ=255)
         gACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCacgtacgt                                                                                                                                             <  1:83477/150‑1 (MQ=255)
         gACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCacgtacgt                                                                                                                                             >  2:416629/1‑150 (MQ=255)
         gACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCacgtacgt                                                                                                                                             >  2:1869726/1‑150 (MQ=255)
         gACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCacgtacgt                                                                                                                                             >  1:1714071/1‑150 (MQ=255)
           ccGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTa                                                                                                                                           >  2:1072218/1‑150 (MQ=255)
           ccGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTa                                                                                                                                           >  2:1288052/1‑150 (MQ=255)
            cGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTaa                                                                                                                                          >  2:490858/1‑150 (MQ=255)
             gATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTaaa                                                                                                                                         >  1:397518/1‑150 (MQ=255)
              aTAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAg                                                                                                                                        <  1:1264465/150‑1 (MQ=255)
                aaTCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGcc                                                                                                                                      <  1:1913294/150‑1 (MQ=255)
                    tgctgTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAAcc                                                                                                                                  >  1:2332217/1‑150 (MQ=255)
                     gctgTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCa                                                                                                                                 >  1:2274623/1‑150 (MQ=255)
                     gctgTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCa                                                                                                                                 >  1:2284091/1‑150 (MQ=255)
                     gctgTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCa                                                                                                                                 >  1:1415954/1‑150 (MQ=255)
                      ctgTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAg                                                                                                                                >  2:1088467/1‑150 (MQ=255)
                              gTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAACGCCAACCAGGATGGCGt                                                                                                                        >  1:685807/1‑150 (MQ=255)
                                tCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTg                                                                                                                      <  1:2248670/150‑1 (MQ=255)
                                tCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTACTAACGCATATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTg                                                                                                                      >  2:2279459/1‑150 (MQ=255)
                                 ccGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTgg                                                                                                                     >  1:1256031/1‑150 (MQ=255)
                                 ccGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTgg                                                                                                                     <  1:1431435/150‑1 (MQ=255)
                                  cGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGa                                                                                                                    <  1:1757032/150‑1 (MQ=255)
                                  cGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGa                                                                                                                    <  2:244791/150‑1 (MQ=255)
                                       aCAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGAcgccg                                                                                                               <  2:1689452/150‑1 (MQ=255)
                                         aGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCa                                                                                                             >  1:214174/1‑150 (MQ=255)
                                              aCCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAg                                                                                                        <  2:2382121/150‑1 (MQ=255)
                                               ccGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAgg                                                                                                       <  2:10835/150‑1 (MQ=255)
                                               ccGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAgg                                                                                                       >  1:1368105/1‑150 (MQ=255)
                                                cGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAggg                                                                                                      <  1:1110393/150‑1 (MQ=255)
                                                 gCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGATAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGt                                                                                                     >  2:2051087/1‑150 (MQ=255)
                                                   tCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTgg                                                                                                   <  1:167139/150‑1 (MQ=255)
                                                    cATCGCCTCTTTCACTTCGTCAATCAGGGTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGc                                                                                                  <  1:1580216/150‑1 (MQ=255)
                                                     aTCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCg                                                                                                 >  2:32689/1‑150 (MQ=255)
                                                         ccgcTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCgg                                                                                             <  2:88810/150‑1 (MQ=255)
                                                              ttCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCg                                                                                        >  1:43406/1‑150 (MQ=255)
                                                               tCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCgg                                                                                       <  1:993130/150‑1 (MQ=255)
                                                                 aCTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTg                                                                                     <  2:1561068/150‑1 (MQ=255)
                                                                 aCTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTg                                                                                     <  2:920173/150‑1 (MQ=255)
                                                                  cTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGa                                                                                    >  2:663659/1‑150 (MQ=255)
                                                                   ttCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGAt                                                                                   <  2:1963204/150‑1 (MQ=255)
                                                                        cAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATaccac                                                                              >  2:2261766/1‑150 (MQ=255)
                                                                        cAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATaccac                                                                              >  1:192692/1‑150 (MQ=255)
                                                                         aaTCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGCCGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATaccacc                                                                             >  2:177904/1‑150 (MQ=255)
                                                                         aaTCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATaccacc                                                                             >  1:1046246/1‑150 (MQ=255)
                                                                         aaTCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATaccacc                                                                             >  2:746293/1‑150 (MQ=255)
                                                                            cAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTAc                                                                          <  1:1430099/150‑1 (MQ=255)
                                                                                 tataGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAg                                                                     >  2:191597/1‑150 (MQ=255)
                                                                                    aGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAAc                                                                  >  2:715526/1‑150 (MQ=255)
                                                                                      aTGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCg                                                                >  1:2328305/1‑150 (MQ=255)
                                                                                      aTGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCg                                                                <  1:378244/150‑1 (MQ=255)
                                                                                        gACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCgat                                                              <  2:1577751/150‑1 (MQ=255)
                                                                                          cGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCgatga                                                            <  1:1869726/150‑1 (MQ=255)
                                                                                          cGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCgatga                                                            >  2:253334/1‑150 (MQ=255)
                                                                                           ggAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCgatgat                                                           <  2:1173613/150‑1 (MQ=255)
                                                                                           ggAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCgatgat                                                           >  1:482081/1‑150 (MQ=255)
                                                                                            gAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATc                                                          <  2:2169033/150‑1 (MQ=255)
                                                                                            gAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATc                                                          <  1:1317599/150‑1 (MQ=255)
                                                                                               tagtaACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTc                                                       <  1:1560815/150‑1 (MQ=255)
                                                                                                agtaACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCa                                                      <  1:1072218/150‑1 (MQ=255)
                                                                                                   aaCGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACtt                                                   >  1:762074/1‑150 (MQ=255)
                                                                                                     cGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTa                                                 >  2:1237978/1‑150 (MQ=255)
                                                                                                       cAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTTATAAAACATACGCCATAAACGATGATCTTCACTTTAAc                                               <  2:2274623/150‑1 (MQ=255)
                                                                                                       cAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAAc                                               >  2:1600414/1‑150 (MQ=255)
                                                                                                       cAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAAc                                               <  2:2284091/150‑1 (MQ=255)
                                                                                                         gATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACtt                                             >  1:2072895/1‑150 (MQ=255)
                                                                                                         gATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACtt                                             >  1:2064332/1‑150 (MQ=255)
                                                                                                           tCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCg                                           >  2:1528248/1‑150 (MQ=255)
                                                                                                            ccAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGt                                          >  1:903806/1‑150 (MQ=255)
                                                                                                             cAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTc                                         <  2:1835765/150‑1 (MQ=255)
                                                                                                             cAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTc                                         <  1:655472/150‑1 (MQ=255)
                                                                                                              aGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCa                                        <  2:1628644/150‑1 (MQ=255)
                                                                                                               ggCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAg                                       >  2:610582/1‑150 (MQ=255)
                                                                                                                 cTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGTGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTa                                     <  2:1345213/150‑1 (MQ=255)
                                                                                                                  tttCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAg                                    <  1:507627/150‑1 (MQ=255)
                                                                                                                    tCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGCACCGATGATCTTCACTTTAACTTCGTCAGTAGAc                                  >  2:92457/1‑150 (MQ=255)
                                                                                                                      cGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAg                                >  2:971423/1‑150 (MQ=255)
                                                                                                                      cGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAg                                >  1:1445415/1‑150 (MQ=255)
                                                                                                                             aTCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGc                         >  2:761579/1‑150 (MQ=255)
                                                                                                                                aacgtACGTGCTGATGCCTCTGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAg                      >  1:233540/5‑150 (MQ=255)
                                                                                                                                 cTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAgg                     <  2:647095/150‑1 (MQ=255)
                                                                                                                                 cTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAgg                     <  2:829905/150‑1 (MQ=255)
                                                                                                                                  cgtACGTGCTGATGCCTCTGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGa                    >  1:1964740/3‑150 (MQ=255)
                                                                                                                                  cgtACGTGCTGATGCCTCTGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGa                    >  1:1609733/3‑150 (MQ=255)
                                                                                                                                   ttACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAg                   >  1:1980090/1‑150 (MQ=255)
                                                                                                                                      cGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCg                >  1:1835718/1‑150 (MQ=255)
                                                                                                                                       gTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCgg               <  1:655389/150‑1 (MQ=255)
                                                                                                                                       gTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCgg               <  2:1415954/150‑1 (MQ=255)
                                                                                                                                        tGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCGga              >  2:1030013/1‑150 (MQ=255)
                                                                                                                                          cAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCGgaga            <  2:1256031/150‑1 (MQ=255)
                                                                                                                                           agagGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCGGAGAt           >  1:802461/1‑150 (MQ=255)
                                                                                                                                            gagGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCGGAGATc          <  2:1718412/150‑1 (MQ=255)
                                                                                                                                              ggCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGTCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCGGAGATCGc        >  1:478169/1‑150 (MQ=255)
                                                                                                                                                 aTCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCGGAGATCGCTTc     >  1:30802/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTCCGCCAGACCGATAATCTGCTGTTTCAGTTCCGGAGACAGCATACCGCTCATCGCCGCTTTCACTTCGTCAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGCAGAGGCATCAGCACGTACGTTAAAGCCAACCAGGATGGCGTTGGACGCCGCAGCCAGGGTGGCGTCGGTTTCGGTGATACCACCTACGCCAGAACCGATGATCTTCACTTTAACTTCGTCAGTAGACAGTTTCAGCAAGGAGTCGGAGATCGCTTCGAC  >  NC_000913/3313617‑3313914

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