Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,097,480 A→C 8.0% intergenic (+651/+85) ycdU → / ← serX putative inner membrane protein/tRNA‑Ser

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,097,4800AC 8.0% 246.7 / 13.6 81intergenic (+651/+85)ycdU/serXputative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (42/32);  new base C (3/4);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTC  >  NC_000913/1097332‑1097628
                                                                                                                                                    |                                                                                                                                                     
tgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGaa                                                                                                                                                       >  3:204541/1‑149 (MQ=11)
tgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGaa                                                                                                                                                       <  1:767593/149‑1 (MQ=11)
tgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGaa                                                                                                                                                       >  3:601396/1‑149 (MQ=11)
tgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGaa                                                                                                                                                       <  1:1460352/149‑1 (MQ=11)
  cAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAAcc                                                                                                                                                     >  1:458805/1‑149 (MQ=32)
   aaTATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGaa                                                                                                                                                       <  1:124211/146‑1 (MQ=11)
   aaTATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGaa                                                                                                                                                       >  2:124211/1‑146 (MQ=11)
      ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCgg                                                                                                                                                 <  4:1192907/149‑1 (MQ=31)
      ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCgg                                                                                                                                                 >  1:913120/1‑149 (MQ=31)
       taCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGt                                                                                                                                                >  2:917680/1‑149 (MQ=11)
                                                                             gtatattgcacacaccggggggatgagaagccccgaccggggttcgacaactggctgtcgccagttgTCGAAc                                                                                                                                                      >  1:1287198‑M2/68‑73 (MQ=255)
                                                                             gtatattgcacacaccggggggatgagaagccccgaccggggttcgacaactggctgtcgccagttgTCGAAc                                                                                                                                                      <  2:1287198‑M2/6‑1 (MQ=255)
                                                                                                        cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTc                                              <  2:1016040/149‑1 (MQ=11)
                                                                                                         tcTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGCGGGGATTCg                                             <  4:1352876/148‑1 (MQ=14)
                                                                                                         aCTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCg                                             >  3:933836/1‑149 (MQ=31)
                                                                                                           tGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGaa                                           >  2:685007/1‑149 (MQ=32)
                                                                                                           tGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGaa                                           <  2:458805/149‑1 (MQ=32)
                                                                                                            gCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAc                                          >  1:439812/1‑149 (MQ=33)
                                                                                                            gCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAc                                          <  3:681923/149‑1 (MQ=33)
                                                                                                               ttCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAcccc                                       <  2:1064191/149‑1 (MQ=35)
                                                                                                               ttCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAcccc                                       <  1:1543021/149‑1 (MQ=35)
                                                                                                                tCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAccccc                                      <  3:44137/149‑1 (MQ=37)
                                                                                                                  aCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGG‑TTCGAAtccccccc                                   <  2:913120/149‑4 (MQ=12)
                                                                                                                     gTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATAc                                 <  3:1544258/149‑1 (MQ=40)
                                                                                                                      tCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACg                                >  1:645316/1‑149 (MQ=255)
                                                                                                                       cTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGt                               <  4:105675/149‑1 (MQ=255)
                                                                                                                          tCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTc                                              >  4:1531657/1‑131 (MQ=11)
                                                                                                                          tCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTc                                              <  3:1531657/131‑1 (MQ=11)
                                                                                                                          tCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGc                            >  1:16322/1‑149 (MQ=255)
                                                                                                                          tCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGc                            >  4:491525/1‑149 (MQ=255)
                                                                                                                          tCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGc                            >  1:1267009/1‑149 (MQ=255)
                                                                                                                          gacaactggcgacagccagttgTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAAGCCCGTa                                                                                            <  2:390781‑M2/64‑1 (MQ=255)
                                                                                                                          gacaactggcgacagccagttgTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAAGCCCGTa                                                                                            >  1:390781‑M2/23‑86 (MQ=255)
                                                                                                                          gacaactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGc                            >  1:1230799‑M2/23‑149 (MQ=255)
                                                                                                                           ccAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGcc                           <  3:1339552/149‑1 (MQ=255)
                                                                                                                           acaactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTATGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGcc                           >  2:27679‑M2/22‑149 (MQ=255)
                                                                                                                           acaactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCtcttc                                                                          >  1:851918‑M2/22‑102 (MQ=255)
                                                                                                                           acaactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGcc                           >  3:1393406‑M2/22‑149 (MQ=255)
                                                                                                                           acaactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTAGAACCCCCGATACGTTGcc                           >  3:83753‑M2/22‑149 (MQ=255)
                                                                                                                           acaactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCAGTGTGTGCGATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCtcttc                                                                          <  2:851918‑M2/81‑1 (MQ=255)
                                                                                                                            caactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCg                          >  1:236514‑M2/21‑149 (MQ=255)
                                                                                                                             aactggcgacagccagttgTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAggggggg                                                  <  4:854158‑M2/105‑1 (MQ=255)
                                                                                                                             aactggcgacagccagttgTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAggggggg                                                  >  3:854158‑M2/20‑124 (MQ=255)
                                                                                                                             aactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGgtgt                                                                                                                        <  4:468059‑M2/36‑1 (MQ=255)
                                                                                                                             aactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGgtgt                                                                                                                        >  3:468059‑M2/20‑55 (MQ=255)
                                                                                                                             aactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGaa                                           >  4:522928‑M2/20‑131 (MQ=255)
                                                                                                                             aactggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGaa                                           <  3:522928‑M2/112‑1 (MQ=255)
                                                                                                                             aaCTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGt                         >  4:1536060/1‑149 (MQ=255)
                                                                                                                              aCTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGta                        <  2:645316/149‑1 (MQ=255)
                                                                                                                               ctggcgacagccagttgTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACG‑aaaaaaaa                                                                                                   >  3:967501‑M2/18‑73 (MQ=255)
                                                                                                                               ctggcgacagccagttgTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACG‑aaaaaaaa                                                                                                   <  4:967501‑M2/56‑1 (MQ=255)
                                                                                                                               ctggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGccc                                                                                               >  3:676397‑M2/18‑77 (MQ=255)
                                                                                                                               ctggcgacagccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGccc                                                                                               <  4:676397‑M2/60‑1 (MQ=255)
                                                                                                                                tGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtat                       >  4:1390932/1‑149 (MQ=255)
                                                                                                                                tGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTTGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtata                      <  2:1366786/149‑1 (MQ=255)
                                                                                                                                tGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtata                      >  2:727464/1‑149 (MQ=255)
                                                                                                                                tGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtata                      <  1:1286611/149‑1 (MQ=255)
                                                                                                                                tGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtata                      <  3:408857/149‑1 (MQ=255)
                                                                                                                                  gCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATaca                    >  1:1424019/1‑149 (MQ=255)
                                                                                                                                  gCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACATTGCCGTATaca                    <  3:115916/149‑1 (MQ=255)
                                                                                                                                    tGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATacaca                  >  1:753232/1‑149 (MQ=255)
                                                                                                                                      tCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACt                >  1:76380/1‑149 (MQ=255)
                                                                                                                                      tCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACt                >  2:404340/1‑149 (MQ=255)
                                                                                                                                        gccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAat                                                                   >  1:2436‑M2/9‑96 (MQ=255)
                                                                                                                                        gccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAat                                                                   <  2:2436‑M2/88‑1 (MQ=255)
                                                                                                                                        gccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTc                                              <  2:89204‑M2/109‑1 (MQ=255)
                                                                                                                                        gccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTc                                              >  1:89204‑M2/9‑117 (MQ=255)
                                                                                                                                        gccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACttt              >  3:259436‑M2/9‑149 (MQ=255)
                                                                                                                                        gccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACttt              >  1:654248‑M2/9‑149 (MQ=255)
                                                                                                                                         ccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTc             >  2:630280‑M2/8‑149 (MQ=255)
                                                                                                                                         ccagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTc             >  3:618338‑M2/8‑149 (MQ=255)
                                                                                                                                          cagttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGACGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTc             >  3:140603‑M2/7‑149 (MQ=255)
                                                                                                                                           agttgTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGa                                                 <  2:1109911‑M2/106‑1 (MQ=255)
                                                                                                                                           agttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCa           <  4:933836‑M2/144‑1 (MQ=255)
                                                                                                                                           agttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCa           <  2:304392‑M2/144‑1 (MQ=255)
                                                                                                                                            gttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAg          <  2:16322‑M2/145‑1 (MQ=255)
                                                                                                                                             ttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCtctt                                                                           >  3:922919‑M2/4‑83 (MQ=255)
                                                                                                                                             ttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCtctt                                                                           <  4:922919‑M2/80‑1 (MQ=255)
                                                                                                                                             ttgTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAgg         >  3:503092‑M2/4‑149 (MQ=255)
                                                                                                                                               gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGt                                                           >  2:136047‑M2/2‑97 (MQ=255)
                                                                                                                                               gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGt                                                           <  1:136047‑M2/96‑1 (MQ=255)
                                                                                                                                               gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCg       >  1:1049333‑M2/2‑149 (MQ=255)
                                                                                                                                               gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCg       >  3:166871‑M2/2‑149 (MQ=255)
                                                                                                                                                tCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGt      <  2:880359/149‑1 (MQ=255)
                                                                                                                                                  gAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtata                      >  2:147667/1‑131 (MQ=255)
                                                                                                                                                  gAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtata                      <  1:147667/131‑1 (MQ=255)
                                                                                                                                                  gAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGc    >  1:1133266/1‑149 (MQ=255)
                                                                                                                                                   aaCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGaa                                           <  2:1263112/109‑1 (MQ=33)
                                                                                                                                                   aaCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGaa                                           >  1:1263112/1‑109 (MQ=33)
                                                                                                                                                    aCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTc  <  1:774113/149‑1 (MQ=255)
                                                                                                                                                    aCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTc  >  1:1470968/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACG‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTC  >  NC_000913/1097332‑1097628

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