Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,923,719 G→C 5.5% intergenic (‑39/+25) rsmG ← / ← mnmG 16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein/5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,923,7190GC 5.5% 476.1 / 10.7 147intergenic (‑39/+25)rsmG/mnmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein/5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base G (85/52);  new base C (5/3);  total (91/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 147 total aligned reads displayed.

GCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAA  >  NC_000913/3923571‑3923865
                                                                                                                                                    |                                                                                                                                                  
gCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGt                                                                                                                                                   >  2:636792/1‑150 (MQ=255)
        ccAATTATGCAGCATATTCACGTAGGCAATAAGCAGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTTATGCGTTGCCTGGTAAGCACCCGCTTACCa                                                                                                                                           <  1:191556/150‑1 (MQ=255)
         cATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCCGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCACCCGCTTACCAg                                                                                                                                          <  1:169760/150‑1 (MQ=255)
         cATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAg                                                                                                                                          <  2:1436668/150‑1 (MQ=255)
         cATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATAACTGATTCATTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGAGCTTACCAg                                                                                                                                          >  2:1304115/1‑150 (MQ=255)
          aTTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGACGGTGCTTACCAgg                                                                                                                                         >  1:207674/1‑150 (MQ=255)
           tttATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCATCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGc                                                                                                                                        >  2:253738/1‑150 (MQ=255)
              aTGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGTTCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCAtt                                                                                                                                     >  1:662936/1‑150 (MQ=255)
                 cagcaTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGATTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCAttttt                                                                                                                                  <  1:1664114/150‑1 (MQ=255)
                 cagcaTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTTCTTACCAGGCAttttt                                                                                                                                  <  2:1051960/150‑1 (MQ=255)
                    caTATTCACGTAGGCAATTTGCTGGTTTTTCTGGTGATCGGTTAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTGGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAAt                                                                                                                               <  2:1355521/150‑1 (MQ=255)
                    caTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAAt                                                                                                                               <  2:687536/150‑1 (MQ=255)
                    caTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTATTGATGCGTTGCCTGGTAATCGGGTGCTTACCAGGCATTTTTAAt                                                                                                                               >  1:294820/1‑150 (MQ=255)
                     atatTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATg                                                                                                                              <  2:355526/150‑1 (MQ=255)
                     atatTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGGCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATg                                                                                                                              <  1:135653/150‑1 (MQ=255)
                              tAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATgcg                                                                                                                     >  1:1629529/1‑150 (MQ=255)
                               aGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATgcgc                                                                                                                    <  1:1730791/150‑1 (MQ=255)
                                   aaTAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTacg                                                                                                                <  2:587766/150‑1 (MQ=255)
                                     tAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATTCCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGACCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTacgac                                                                                                              <  1:363988/150‑1 (MQ=255)
                                       aGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGTAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgc                                                                                                            <  2:425992/150‑1 (MQ=255)
                                       aGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgc                                                                                                            >  2:317157/1‑150 (MQ=255)
                                         cTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcag                                                                                                          >  1:1356712/1‑150 (MQ=255)
                                          tGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcagc                                                                                                         >  1:289520/1‑150 (MQ=255)
                                           ggTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcagca                                                                                                        >  1:406113/1‑150 (MQ=255)
                                            gTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCAt                                                                                                       >  1:439347/1‑150 (MQ=255)
                                            gTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCGTACCAGGCATTTTTTATGCGTTATGCGCTACGACGCAGCAt                                                                                                       >  1:1268357/1‑150 (MQ=255)
                                              ttttCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATAc                                                                                                     <  2:1153735/150‑1 (MQ=255)
                                                 tCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCACTACGACGCAGCATACCCt                                                                                                  >  2:465081/1‑150 (MQ=255)
                                                    ggTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGtt                                                                                               >  2:398003/1‑150 (MQ=255)
                                                       gATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGttttt                                                                                            >  2:1547803/1‑150 (MQ=255)
                                                       gATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGttttt                                                                                            <  1:1035763/150‑1 (MQ=255)
                                                         tCGGTAAGCGAACTACCTGCGTCTTTCGGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTc                                                                                          <  1:302654/150‑1 (MQ=255)
                                                          cGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATTCGCTACGACGCAGCATACCCTGTTTTTTCa                                                                                         >  1:145759/1‑150 (MQ=255)
                                                           ggTAAGCGAAATACCTGCGTCATTCAGAAGTAAGGAGAGTTTGTTGAGCACTGTGATAACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATATTTAATGCGTTAAGAGCTACTCCGCTCCATACCCTCTTTATTCAg                                                                                        >  1:959485/1‑150 (MQ=255)
                                                            gTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGGTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGc                                                                                       <  1:981519/150‑1 (MQ=255)
                                                               aGCGAAATACCTTCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCcac                                                                                    <  1:246905/150‑1 (MQ=255)
                                                                gCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCATCcaca                                                                                   >  1:220524/1‑150 (MQ=255)
                                                                  gAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACAcc                                                                                 <  1:738887/150‑1 (MQ=255)
                                                                  gAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACAcc                                                                                 >  2:1470339/1‑150 (MQ=255)
                                                                     aTACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACcagc                                                                              <  1:1971962/150‑1 (MQ=255)
                                                                       aCCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACcagcag                                                                            <  1:1695855/150‑1 (MQ=255)
                                                                           gCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATg                                                                        >  1:26803/1‑150 (MQ=255)
                                                                           gCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATg                                                                        >  2:582195/1‑150 (MQ=255)
                                                                            cGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATgg                                                                       >  2:153377/1‑150 (MQ=255)
                                                                            cGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATgg                                                                       <  2:516905/150‑1 (MQ=255)
                                                                              tCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGgag                                                                     >  2:301954/1‑150 (MQ=255)
                                                                               cTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCCGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCCGCCACCCCAGCAGAATGgaga                                                                    >  2:1211396/1‑150 (MQ=255)
                                                                               cTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGgaga                                                                    >  2:1640782/1‑150 (MQ=255)
                                                                                 ttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCCGCCACACCAGCAGAATGGAGATg                                                                  >  2:391523/1‑150 (MQ=255)
                                                                                 ttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATg                                                                  >  1:733215/1‑150 (MQ=255)
                                                                                   cagcagTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATgac                                                                <  1:235470/150‑3 (MQ=255)
                                                                                    agcagTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGcc                                                               >  2:847493/1‑150 (MQ=255)
                                                                                     gcagTAAGGAGAGTTTGGTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCAGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCg                                                              <  2:214223/150‑1 (MQ=255)
                                                                                      cagTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGATATTGCCGc                                                             >  1:1886409/1‑150 (MQ=255)
                                                                                        gTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAg                                                           <  1:999720/150‑1 (MQ=255)
                                                                                          aaGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGc                                                         >  2:1127250/1‑150 (MQ=255)
                                                                                          aaGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGAATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATTGCCGCAGGc                                                         >  2:793547/1‑150 (MQ=255)
                                                                                           aGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCg                                                        >  2:398864/1‑150 (MQ=255)
                                                                                           aGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGAAGAATGGAGATGGCCGCAGGCg                                                        >  1:10224/1‑150 (MQ=255)
                                                                                            ggAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGt                                                       >  2:2009535/1‑150 (MQ=255)
                                                                                             gagaGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTg                                                      <  2:969337/150‑1 (MQ=255)
                                                                                               gagTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTGGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGAc                                                    <  1:1832224/150‑1 (MQ=255)
                                                                                                     gttgAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATTCGCTACGACGCAGCATACCCAGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGa                                              >  1:851729/1‑150 (MQ=255)
                                                                                                      ttgAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGaa                                             <  1:950495/150‑1 (MQ=255)
                                                                                                      ttgAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGaa                                             <  2:1657091/150‑1 (MQ=255)
                                                                                                       tgAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTAGCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGaaa                                            <  2:176348/150‑1 (MQ=255)
                                                                                                       tgAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGaaa                                            >  2:766053/1‑150 (MQ=255)
                                                                                                        gAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAAt                                           >  1:647533/1‑150 (MQ=255)
                                                                                                           cACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATAcg                                        >  2:1683365/1‑150 (MQ=255)
                                                                                                             cGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATAcgcg                                      <  1:1198451/150‑1 (MQ=255)
                                                                                                             cGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATAcgcg                                      <  1:1362874/150‑1 (MQ=255)
                                                                                                                    tACCTGTTCTTGATGCGTTGCCAGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAAAACGCGAAGCGTg                               >  2:943120/1‑150 (MQ=255)
                                                                                                                      ccTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCGTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATAACCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGc                             >  1:1002021/1‑150 (MQ=255)
                                                                                                                       cTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGcc                            <  1:1943807/150‑1 (MQ=255)
                                                                                                                       cTGTTCTTGAGGCGTTTCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGcc                            <  2:835928/150‑1 (MQ=255)
                                                                                                                        tGTACTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCGAGAAATACGCGAAGCGTGGCCg                           >  1:89446/1‑150 (MQ=255)
                                                                                                                         gTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGa                          >  1:1110947/1‑150 (MQ=255)
                                                                                                                         gTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGa                          >  2:1837916/1‑150 (MQ=255)
                                                                                                                          ttCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGAt                         >  1:1893416/1‑150 (MQ=255)
                                                                                                                             ttGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATaga                      >  1:759647/1‑150 (MQ=255)
                                                                                                                             ttGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATaga                      <  1:1223753/150‑1 (MQ=255)
                                                                                                                              tGATGCGTTGCCTGGTGAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTTCGACGCAGCATACCATGTTTTTTCAGCCACACAAGAAGTATGGTGATGTCCGCAGGCGTGACGCCAGAAATACGCGTAGCTTGGCCGATagag                     >  1:1283585/1‑150 (MQ=255)
                                                                                                                               gATGCGTTGCCTGGTAAGCGGGTGCATACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTAAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAgg                    <  2:279697/150‑1 (MQ=255)
                                                                                                                                aTGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGc                   <  1:1357886/150‑1 (MQ=255)
                                                                                                                                 tGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACGCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCt                  >  1:304444/1‑150 (MQ=255)
                                                                                                                                   cGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTgg                >  2:35914/1‑150 (MQ=255)
                                                                                                                                     ttGGCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGtt              <  2:182073/150‑1 (MQ=255)
                                                                                                                                       gCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtg            >  2:1869432/1‑150 (MQ=255)
                                                                                                                                        ccTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgt           >  2:226624/1‑150 (MQ=255)
                                                                                                                                        ccTGGTAAGCACCCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgt           >  2:1417079/1‑150 (MQ=255)
                                                                                                                                        ccTGGTAAGCACCCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATAACCTGAGTTTTCACCCACACCAGCAGAATGGAGATGGCCGCAGCCGTGACGCCAGAAATACGCGCAGCTTGGCCGATAGAGGCTGGTTtgt           >  1:976746/1‑150 (MQ=255)
                                                                                                                                         cTGGTAAGCGGGTGCTTACCAGGCATTTTTAATTCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgtg          <  1:1386667/150‑1 (MQ=255)
                                                                                                                                         cTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgtg          <  2:274561/150‑1 (MQ=255)
                                                                                                                                         cTGGTAAGCACCCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgtg          >  2:496688/1‑150 (MQ=255)
                                                                                                                                          tGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGa         >  2:1283047/1‑150 (MQ=255)
                                                                                                                                          tGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGa         >  1:1090400/1‑150 (MQ=255)
                                                                                                                                             tAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCg      <  2:1125769/150‑1 (MQ=255)
                                                                                                                                             tAAGCGGGTGCTTAACAGTCATTTTTAATGCGTTATGCGCTACGACTCAGCATACCCTGTTTTTTCAGCCACACCTGCAGAATGGAGATGGCCGCTGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCg      <  1:334297/150‑1 (MQ=255)
                                                                                                                                              aaGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGt     >  1:893619/1‑150 (MQ=255)
                                                                                                                                                 cGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTaa  >  2:1297936/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAA  >  NC_000913/3923571‑3923865

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