Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,845,750 T→C 5.1% intergenic (+20/+26) adeD → / ← uhpT cryptic adenine deaminase/hexose phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,845,7500TC 5.1% 442.6 / 7.1 134intergenic (+20/+26)adeD/uhpTcryptic adenine deaminase/hexose phosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (72/55);  new base C (4/3);  total (76/58)
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.93e-01

Only 100 of 135 total aligned reads displayed.

CCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCC  >  NC_000913/3845605‑3845898
                                                                                                                                                 |                                                                                                                                                    
 ccGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCg                                                                                                                                                  >  4:808353/1‑149 (MQ=255)
 ccGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCg                                                                                                                                                  >  2:842914/1‑149 (MQ=255)
  cGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCgg                                                                                                                                                 <  3:324645/149‑1 (MQ=255)
   gATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCggg                                                                                                                                                <  2:760333/149‑1 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCaaaa                                                                                                                                           <  2:605512/149‑1 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCaaaa                                                                                                                                           <  2:134357/149‑1 (MQ=255)
         ccGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                                                                                                          >  4:1194722/1‑149 (MQ=255)
         ccGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCATGACGCCGGGcttttt                                                                                                                                          <  3:842073/149‑6 (MQ=255)
         ccGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCATGACGCCGGGcttttt                                                                                                                                          <  2:72827/149‑6 (MQ=255)
           gTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGTCGTCATGCCGGGCAAAAGTc                                                                                                                                        >  2:188887/1‑149 (MQ=255)
            tttATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCa                                                                                                                                       <  1:139774/149‑1 (MQ=255)
             ttATTCAGATGGCGTTTCTTTCTCCGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCAc                                                                                                                                      <  1:1050503/149‑1 (MQ=255)
                  cAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGtt                                                                                                                                 >  2:1469254/1‑149 (MQ=255)
                   aGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTa                                                                                                                                >  1:1121191/1‑149 (MQ=255)
                   aGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTa                                                                                                                                >  1:495806/1‑149 (MQ=255)
                          gTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTAt                                                                                                                         >  2:1237688/1‑149 (MQ=255)
                           tttctttcTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATg                                                                                                                        >  2:880685/1‑149 (MQ=255)
                           tttctttcTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATg                                                                                                                        <  2:30128/149‑1 (MQ=255)
                                ttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACt                                                                                                                   >  2:236438/1‑149 (MQ=255)
                                ttctCTGCCAGAGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACt                                                                                                                   >  3:867315/1‑149 (MQ=255)
                                 tctcTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTg                                                                                                                  >  4:117854/1‑149 (MQ=255)
                                   tctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTc                                                                                                                >  3:1198175/1‑149 (MQ=255)
                                   tctGCCAGTAATCCCCGCGCTAAAACTAGCCAGTCGGGGGCTATTTGATGGCGAGACGTTTGCCTTCACTACGCTAGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGta                                                                                                                >  4:1071709/1‑148 (MQ=255)
                                      gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAc                                                                                                             >  1:560395/1‑149 (MQ=255)
                                       ccAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                                                                                            >  3:30755/1‑149 (MQ=255)
                                             aTCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGGCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGaa                                                                                                      <  1:1364018/149‑1 (MQ=255)
                                             aTCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCa                                                                                                                                       >  4:479712/1‑116 (MQ=255)
                                             aTCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGaa                                                                                                      <  1:658841/149‑1 (MQ=255)
                                             aTCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGGCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGaa                                                                                                      <  2:829537/149‑1 (MQ=255)
                                              tCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAt                                                                                                     <  1:1495786/149‑1 (MQ=255)
                                                cccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttt                                                                                                   >  1:623207/1‑149 (MQ=255)
                                                 ccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACTGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtttt                                                                                                  >  4:1212679/1‑149 (MQ=255)
                                                 ccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtttt                                                                                                  <  3:879059/149‑1 (MQ=255)
                                                  cgcgCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttttt                                                                                                 <  3:808353/149‑1 (MQ=255)
                                                  cgcgCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttttt                                                                                                 <  2:623207/149‑1 (MQ=255)
                                                  cgcgCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttttt                                                                                                 <  2:495806/149‑1 (MQ=255)
                                                       tAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtctc                                                                                            <  2:577750/149‑1 (MQ=255)
                                                            cTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggc                                                                                       >  3:1485183/1‑149 (MQ=255)
                                                             tAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggcg                                                                                      >  3:1065160/1‑149 (MQ=255)
                                                               aCCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGa                                                                                    <  2:101292/149‑1 (MQ=255)
                                                                    tCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTg                                                                               <  1:7212/149‑1 (MQ=255)
                                                                       ggggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGt                                                                            <  3:960639/149‑1 (MQ=255)
                                                                         ggCTATTTGATGGCGAGAAGTTTGCCTTCACTTCGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGgtc                                                                          <  2:556728/149‑1 (MQ=255)
                                                                         ggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttc                                                                          >  4:591448/1‑149 (MQ=255)
                                                                            tATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttcttc                                                                       >  2:1416731/1‑149 (MQ=255)
                                                                            tATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttcttc                                                                       >  4:1436151/1‑149 (MQ=255)
                                                                             aTTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTcc                                                                      <  1:1237688/149‑1 (MQ=255)
                                                                             aTTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTcc                                                                      <  3:1436151/149‑1 (MQ=255)
                                                                              tttGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCa                                                                     <  4:1375829/149‑1 (MQ=255)
                                                                                  aTGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTAc                                                                 >  1:1110626/1‑149 (MQ=255)
                                                                                  aTGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTAc                                                                 <  3:710653/149‑1 (MQ=255)
                                                                                    ggCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgc                                                                                          <  2:333675/122‑1 (MQ=255)
                                                                                    ggCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgc                                                                                          >  1:333675/1‑122 (MQ=255)
                                                                                       gcgaAGTTTGCCTTCACTACGCTGGAAGTCCCGGAATGATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCa                                                            >  3:595314/3‑149 (MQ=255)
                                                                                             tttGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCg                                                      >  2:1196772/1‑149 (MQ=255)
                                                                                             attGCCATCACTACGCTCGAGGACACGAAAAAATAAAAAAAGCCTGGCGTCATGCCGGGCAAACGTCACCAGTTACGTTGATGCCACAGACCACTGCTGAATTTTTTTCTCGCGGCGGATATTGCGTTCTTCCATTACCGCCACTATCg                                                      >  4:172143/2‑149 (MQ=255)
                                                                                                gCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAGGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCa                                                   <  2:173004/149‑1 (MQ=255)
                                                                                                  cTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATc                                                 >  3:14403/1‑149 (MQ=255)
                                                                                                   ttCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCa                                                <  2:517455/149‑1 (MQ=255)
                                                                                                   ttCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCa                                                <  1:929014/149‑1 (MQ=255)
                                                                                                   ttCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCa                                                >  3:625682/1‑149 (MQ=255)
                                                                                                   ttCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCa                                                >  1:1392024/1‑149 (MQ=255)
                                                                                                   ttCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCa                                                >  1:1048300/1‑149 (MQ=255)
                                                                                                     cACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAgg                                              <  2:1033198/149‑1 (MQ=255)
                                                                                                      aCTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttcttc                                                                       >  2:818021/1‑123 (MQ=255)
                                                                                                      aCTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGc                                             >  4:696273/1‑149 (MQ=255)
                                                                                                       cTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCa                                            <  2:1048300/149‑1 (MQ=255)
                                                                                                        tACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttc                                                                          >  2:736899/1‑118 (MQ=255)
                                                                                                        tACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttc                                                                          <  1:736899/118‑1 (MQ=255)
                                                                                                        tACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAg                                           >  1:773665/1‑149 (MQ=255)
                                                                                                        tACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAg                                           >  2:1189945/1‑149 (MQ=255)
                                                                                                           gCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATa                                        >  4:855375/1‑149 (MQ=255)
                                                                                                             tGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACa                                      >  3:512242/1‑149 (MQ=255)
                                                                                                             tGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACa                                      >  2:1399986/1‑149 (MQ=255)
                                                                                                              ggAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACaa                                     <  4:1002430/149‑1 (MQ=255)
                                                                                                              ggAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACaa                                     <  2:1121191/149‑1 (MQ=255)
                                                                                                              ggAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTAAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACaa                                     <  3:890030/149‑1 (MQ=255)
                                                                                                                     aCGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATcg                              >  1:1012238/1‑149 (MQ=255)
                                                                                                                     aCGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATcg                              >  3:122091/1‑149 (MQ=255)
                                                                                                                     aCGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGCCAGATACAACCAATcg                              >  1:149327/1‑149 (MQ=255)
                                                                                                                       ggAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATcgcg                            >  1:426121/1‑149 (MQ=255)
                                                                                                                             ataaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGata                      >  2:998425/1‑149 (MQ=255)
                                                                                                                              taaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGatat                     >  2:900439/1‑149 (MQ=255)
                                                                                                                              taaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGatat                     >  3:280606/1‑149 (MQ=255)
                                                                                                                              taaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGatat                     >  3:28396/1‑149 (MQ=255)
                                                                                                                              taaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGatat                     >  3:335096/1‑149 (MQ=255)
                                                                                                                               aaaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATc                    >  3:1159133/1‑149 (MQ=255)
                                                                                                                               aaaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATc                    >  3:535356/1‑149 (MQ=255)
                                                                                                                                aaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATcc                   <  3:723889/149‑1 (MQ=255)
                                                                                                                                aaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATcc                   >  1:342009/1‑149 (MQ=255)
                                                                                                                                aaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATcc                   <  4:1485183/149‑1 (MQ=255)
                                                                                                                                aaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATcc                   <  3:1194722/149‑1 (MQ=255)
                                                                                                                                 aaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCa                  >  1:929697/1‑149 (MQ=255)
                                                                                                                                  ttttgCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggc                          >  1:348136/5‑140 (MQ=255)
                                                                                                                                  ttttgCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggc                          <  2:348136/136‑1 (MQ=255)
                                                                                                                                  ttttgCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAg                 >  1:1425960/5‑149 (MQ=255)
                                                                                                                                  ttttgCCCGGCATGACACCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAg                 >  2:285560/5‑149 (MQ=255)
                                                                                                                                        ccGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGa                                          <  3:170643/118‑1 (MQ=255)
                                                                                                                                        ccGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGa                                          >  4:170643/1‑118 (MQ=255)
                                                                                                                                              tCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGt     >  4:355258/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCC  >  NC_000913/3845605‑3845898

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