Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,923,377 C→T 5.3% V102M (GTG→ATG) ‡ rsmG ← 16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,923,3770CT 5.3% 504.0 / 11.3 151V102M (GTG→ATG) ‡rsmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base C (78/65);  new base T (2/6);  total (80/71)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 159 total aligned reads displayed.

GCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATG  >  NC_000913/3923234‑3923524
                                                                                                                                               |                                                                                                                                                   
gCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGtt                                                                                                                                                <  4:807327/149‑1 (MQ=255)
gCTCACCATATCGTTCAGAGAGGCAAAAGCGCGACTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGtt                                                                                                                                                >  3:538930/1‑149 (MQ=255)
   cACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCGTGCGtttcc                                                                                                                                             <  1:828741/149‑3 (MQ=255)
   cACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCGTGCGtttcc                                                                                                                                             <  2:400139/149‑3 (MQ=255)
    aCCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTAc                                                                                                                                             <  3:223826/148‑1 (MQ=255)
    aCCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTAc                                                                                                                                             >  4:223826/1‑148 (MQ=255)
      cATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCGTGCGtttccttc                                                                                                                                          <  3:632812/149‑6 (MQ=255)
      cATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCaa                                                                                                                                          >  2:617391/1‑149 (MQ=255)
      cATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCaa                                                                                                                                          <  4:632394/149‑1 (MQ=255)
       atatCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAg                                                                                                                                         <  3:738387/149‑1 (MQ=255)
         atCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGc                                                                                                                                       >  3:985124/1‑149 (MQ=255)
         atCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGc                                                                                                                                       >  2:336/1‑149 (MQ=255)
             ttCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATc                                                                                                                                   <  4:1002865/149‑1 (MQ=255)
                agagagGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCaa                                                                                                                                >  4:451066/1‑149 (MQ=255)
                agagagGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCaa                                                                                                                                >  1:606020/1‑149 (MQ=255)
                agagagGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCaa                                                                                                                                >  2:999305/1‑149 (MQ=255)
                  agagGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACa                                                                                                                              >  1:1035595/1‑149 (MQ=255)
                   gagGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACag                                                                                                                             >  2:249826/1‑149 (MQ=255)
                      gCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGt                                                                                                                          >  1:440739/1‑149 (MQ=255)
                        aaaaGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGa                                                                                                                        >  2:316843/1‑149 (MQ=255)
                          aaGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGaaa                                                                                                                      <  2:864753/149‑1 (MQ=255)
                          aaGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGaaa                                                                                                                      <  4:902663/149‑1 (MQ=255)
                                ggCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCt                                                                                                                >  2:516148/1‑149 (MQ=255)
                                 gCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCtt                                                                                                               <  2:812221/149‑1 (MQ=255)
                                 gCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCtt                                                                                                               <  1:177958/149‑1 (MQ=255)
                                 gCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCtt                                                                                                               <  2:501496/149‑1 (MQ=255)
                                  cTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTc                                                                                                              >  4:46109/1‑149 (MQ=255)
                                   tAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCa                                                                                                             >  4:171573/1‑149 (MQ=255)
                                      ttACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGa                                                                                                          <  3:46109/149‑1 (MQ=255)
                                          gCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCa                                                                                                      >  2:930760/1‑149 (MQ=255)
                                           ccATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCAc                                                                                                     <  1:954938/149‑1 (MQ=255)
                                             aTCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGa                                                                                                   >  1:680638/1‑149 (MQ=255)
                                             aTCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGa                                                                                                   >  3:17303/1‑149 (MQ=255)
                                              tCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGAt                                                                                                  >  4:762166/1‑149 (MQ=255)
                                                 aaTGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCCGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATaga                                                                                               <  4:8288/149‑1 (MQ=255)
                                                     gcggcTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGt                                                                                           >  2:270890/1‑149 (MQ=255)
                                                          tctGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAg                                                                                                                                         <  1:101800/98‑1 (MQ=255)
                                                              aaGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATTGAGAGTGGAATGCCt                                                                                  >  2:864111/1‑149 (MQ=255)
                                                              aaGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCt                                                                                  >  3:41382/1‑149 (MQ=255)
                                                               aGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTg                                                                                 <  3:291701/149‑1 (MQ=255)
                                                                ggAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTgg                                                                                >  1:403035/1‑149 (MQ=255)
                                                                      tcttctACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTcc                                                                          <  1:617391/149‑1 (MQ=255)
                                                                        ttctACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTg                                                                        >  2:619603/1‑149 (MQ=255)
                                                                         tctACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTgg                                                                       <  1:295307/149‑1 (MQ=255)
                                                                         tctACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTgg                                                                       <  3:568891/149‑1 (MQ=255)
                                                                           tACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTc                                                                     >  2:881617/1‑149 (MQ=255)
                                                                            aCCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTcc                                                                    <  1:784918/149‑1 (MQ=255)
                                                                             cccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCg                                                                   >  3:751136/1‑149 (MQ=255)
                                                                              ccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGa                                                                                                   <  4:683941/116‑1 (MQ=255)
                                                                              ccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCgg                                                                  >  2:594539/1‑149 (MQ=255)
                                                                              ccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCgg                                                                  >  3:23698/1‑149 (MQ=255)
                                                                              ccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCgg                                                                  >  2:358301/1‑149 (MQ=255)
                                                                               cTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGt                                                                 >  1:344875/1‑149 (MQ=255)
                                                                               cTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGt                                                                 <  2:505917/149‑1 (MQ=255)
                                                                                tGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGtt                                                                                                                                                <  4:583171/69‑1 (MQ=255)
                                                                                   tctGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCg                                                             >  1:939148/1‑149 (MQ=255)
                                                                                        acaggTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCGGACATc                                                        <  2:368701/146‑1 (MQ=255)
                                                                                        aCTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATc                                                        <  2:295457/149‑1 (MQ=255)
                                                                                             ttCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATaa                                                   >  2:701613/1‑149 (MQ=255)
                                                                                              tCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCt                                                                                  >  3:811066/1‑117 (MQ=255)
                                                                                              tCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATaaa                                                  <  1:993966/149‑1 (MQ=255)
                                                                                                aaTATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAAcc                                                >  3:576212/1‑149 (MQ=255)
                                                                                                 atatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCg                                               >  3:72873/1‑149 (MQ=255)
                                                                                                  tatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGt                                              >  3:457459/1‑149 (MQ=255)
                                                                                                    ttCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTc                                            <  3:371374/149‑1 (MQ=255)
                                                                                                      ctcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCAc                                          <  3:421456/149‑1 (MQ=255)
                                                                                                       tcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATagaga                                                                                             <  4:160858/97‑1 (MQ=255)
                                                                                                       tcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATagaga                                                                                             >  3:160858/1‑97 (MQ=255)
                                                                                                           gTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCg                                                                   >  1:499613/1‑119 (MQ=255)
                                                                                                           gTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGc                                     >  3:670239/1‑149 (MQ=255)
                                                                                                            tttAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATa                                                    >  3:969428/1‑133 (MQ=255)
                                                                                                                 gCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATAc                               <  2:1035595/149‑1 (MQ=255)
                                                                                                                   tCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCg                                                       >  2:69509/1‑123 (MQ=255)
                                                                                                                   tCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCg                                                       <  1:69509/123‑1 (MQ=255)
                                                                                                                    cATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGt                            <  1:358301/149‑1 (MQ=255)
                                                                                                                    cATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGt                            <  1:881617/149‑1 (MQ=255)
                                                                                                                      tGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGc                          >  3:420758/1‑149 (MQ=255)
                                                                                                                       gTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCtt                                                                                                               >  4:100168/1‑63 (MQ=255)
                                                                                                                       gTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCtt                                                                                                               <  3:100168/63‑1 (MQ=255)
                                                                                                                        ttGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATc                                                        <  3:858093/117‑1 (MQ=255)
                                                                                                                        ttGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATc                                                        >  4:858093/1‑117 (MQ=255)
                                                                                                                        ttGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGcca                        <  2:43446/149‑1 (MQ=255)
                                                                                                                         tGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGccac                       <  2:4647/149‑1 (MQ=255)
                                                                                                                            aCCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCcaccac                    >  3:292698/1‑149 (MQ=255)
                                                                                                                               tGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAAt                 <  4:73190/149‑1 (MQ=255)
                                                                                                                                   gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGaa                                                                                                                       <  1:529/43‑1 (MQ=255)
                                                                                                                                   gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGaa                                                                                                                       >  2:529/1‑43 (MQ=255)
                                                                                                                                   gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGc                                                               >  1:981986/1‑99 (MQ=255)
                                                                                                                                   gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGc                                                               <  2:981986/99‑1 (MQ=255)
                                                                                                                                   gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTa             >  4:376583/1‑149 (MQ=255)
                                                                                                                                   gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTa             >  4:360144/1‑149 (MQ=255)
                                                                                                                                     aGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc           >  4:804091/1‑149 (MQ=255)
                                                                                                                                      ggAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCg          >  3:363932/1‑149 (MQ=255)
                                                                                                                                       gAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCga         <  3:168073/149‑1 (MQ=255)
                                                                                                                                          acgcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGaa      <  2:403035/149‑1 (MQ=255)
                                                                                                                                            gcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCa                                    <  4:489958/117‑1 (MQ=255)
                                                                                                                                            gcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAata    >  3:802904/1‑149 (MQ=255)
                                                                                                                                              acgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTcc                                                                          >  2:46790/1‑77 (MQ=255)
                                                                                                                                              acgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTcc                                                                          <  1:46790/77‑1 (MQ=255)
                                                                                                                                              acgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATg  >  1:384777/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATG  >  NC_000913/3923234‑3923524

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: