Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,397,885 A→G 5.1% N164N (AAT→AAC)  yhdE ← dTTP/UTP pyrophosphatase; m(5)UTP/m(5)CTP/pseudo‑UTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,397,8850AG 5.1% 281.5 / 2.5 110N164N (AAT→AAC) yhdEdTTP/UTP pyrophosphatase; m(5)UTP/m(5)CTP/pseudo‑UTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (47/54);  new base G (3/3);  total (51/59)
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.61e-01

Only 100 of 111 total aligned reads displayed.

AATATACGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAG  >  NC_000913/3397737‑3398031
                                                                                                                                                    |                                                                                                                                                  
aaTATACGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCCCACGCAGTGCTTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCa                                                                                                                                                    <  2:147412/149‑1 (MQ=255)
 atatACGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCAt                                                                                                                                                   >  1:1027409/1‑149 (MQ=255)
  tataCGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGCCTACGCGGTGCTAGCTCCCAtt                                                                                                                                                  >  2:689002/1‑149 (MQ=255)
  tataCGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCAtt                                                                                                                                                  >  1:338531/1‑149 (MQ=255)
  tataCGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCAtt                                                                                                                                                  >  2:768940/1‑149 (MQ=255)
  tataCGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGGGATGGGTGCCCtt                                                                                                                                                  >  2:1151250/1‑149 (MQ=255)
    taCGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTa                                                                                                                                                <  1:820693/149‑1 (MQ=255)
     accccACTCGCGTTTCCGAAGGCGGTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTAt                                                                                                                                               <  2:387465/146‑1 (MQ=255)
      cGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATc                                                                                                                                              <  1:1125900/149‑1 (MQ=255)
         cccTCGCGTCTCCGAAGGCTGTACGTTTACTAACAATCCAGCCGTCGTGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTc                                                                                                                                           <  1:298763/147‑1 (MQ=255)
         cACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTc                                                                                                                                           <  1:607881/149‑1 (MQ=255)
            tCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCt                                                                                                                                         <  1:861208/148‑1 (MQ=255)
            tCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTAAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCt                                                                                                                                         >  2:861208/1‑148 (MQ=255)
                gTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACaa                                                                                                                                    <  2:390736/149‑1 (MQ=255)
                gTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACaa                                                                                                                                    <  1:817502/149‑1 (MQ=255)
                 tttCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTCAAATTACTTAATAATTCATACTTTTCAACCAGCGGGGAGCCGACTACGGCGTGATAGCTGACATTTATCCTCCTGACaaa                                                                                                                                   >  1:299122/1‑149 (MQ=255)
                   tgcGAAGGCGTCACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAAc                                                                                                                                 <  2:372906/147‑1 (MQ=255)
                    ccGAAGGCGTTACGTTTACTAAAAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACa                                                                                                                                <  1:10192/149‑1 (MQ=255)
                     cGAAGGCGTTACGATTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATACCTGCCATTTATCTTCCTGACAAAACAg                                                                                                                               >  1:884709/1‑149 (MQ=255)
                       aaGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGcc                                                                                                                             >  1:682583/1‑149 (MQ=255)
                        aGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCa                                                                                                                            <  1:746383/149‑1 (MQ=255)
                          ctgTTACGTTTACTAACAATTCATCCGTCATGTTTATCCCTTTTCTCACGCAGCGCGTTAAAATTACTTCATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCAcc                                                                                                                          <  1:213543/147‑1 (MQ=255)
                             ttACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAg                                                                                                                       >  2:214725/1‑149 (MQ=255)
                                 gTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGcccc                                                                                                                   <  2:1103836/149‑1 (MQ=255)
                                 gTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGACATTTATCTTCCTGACAAAACAGCCACCCAGccac                                                                                                                   >  1:476930/1‑147 (MQ=255)
                                  tttACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCt                                                                                                                  >  2:1210867/1‑149 (MQ=255)
                                  tttACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCt                                                                                                                  >  1:412122/1‑149 (MQ=255)
                                       tAACAATTCAGCCGTCGTGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTCCGGCGTGATAGATGCCGTTTATCTTCCTGGCAAAACAGCCACCCAGCCCCTGGATa                                                                                                             >  2:492626/1‑149 (MQ=255)
                                       tAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCGTTTATCTTCCTGGCAAAACAGCCACCCAGCCCCTGagtc                                                                                                             >  2:606922/1‑146 (MQ=255)
                                       tAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCGTTTATCTTCCTGGCAAAACAGCCACCCAGCCCCTGaaaa                                                                                                             >  2:1293639/1‑147 (MQ=255)
                                       tAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATa                                                                                                             >  1:1015936/1‑149 (MQ=255)
                                       tAACAAGTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAATTTACTTAATAATGCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTTCTGACAGAACAGCCACCCAGCCCCTGAATa                                                                                                             <  2:890742/149‑1 (MQ=255)
                                        aaCAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCACCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGGCAAAACAGCCACCCAACCCCTGAATAc                                                                                                            >  2:692202/1‑149 (MQ=255)
                                        aaCAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATAc                                                                                                            >  1:1069358/1‑149 (MQ=255)
                                         aCAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCGGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATAcc                                                                                                           >  1:213320/1‑149 (MQ=255)
                                         aCAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATAcc                                                                                                           >  1:201917/1‑149 (MQ=255)
                                         aCAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATAcc                                                                                                           >  2:190018/1‑149 (MQ=255)
                                         aCAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCATTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCGGACAAAACAGCCACCCAGCCCCTGAATAcc                                                                                                           >  2:438734/1‑149 (MQ=255)
                                          cAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCg                                                                                                          <  1:687675/149‑1 (MQ=255)
                                          cAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGGGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCg                                                                                                          <  2:1035432/149‑1 (MQ=255)
                                           aaTTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAAAAGCCACCCAGCCCCTGaaacacgt                                                                                                         >  2:598607/1‑143 (MQ=255)
                                               ccgCCGTCATGTTTATCCCTTTTCTCCCGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGc                                                                                                     >  1:1162837/3‑149 (MQ=255)
                                               cAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTTATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGc                                                                                                     <  1:756438/149‑1 (MQ=255)
                                               cAGCCGTCATGTTTATCCCTTTTCCCACGCAGTGCGTTACAATTACTTCATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGc                                                                                                     <  1:756165/149‑1 (MQ=255)
                                                   cGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCCCCt                                                                                                 >  2:714711/1‑149 (MQ=255)
                                                   cGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCt                                                                                                 <  2:1251999/149‑1 (MQ=255)
                                                     tCATGTTTATCCCTTTTCTCACGCTGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGc                                                                                               <  1:800901/149‑1 (MQ=255)
                                                     tCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGc                                                                                               <  1:1286232/149‑1 (MQ=255)
                                                       aTGTTTATCCCTTGTCCCACGCAGTGCGTTAAAAGTACTTAATAATTCATACGTTTCACCCCGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCtt                                                                                             <  1:624723/149‑1 (MQ=255)
                                                       aTGGTTATCCCTTTTCTCACGCAGTCCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGATAAGCCGACTACGGCGTGATGGCTGCCATGGATCTTCCTCGCAAAACAGCCACCCAGCCCCTGAATACCGTCTGCCCCGGCtt                                                                                             >  2:426481/1‑149 (MQ=255)
                                                           ttATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATCGTTCATACGCTTCAGCCATCGGGAAGCCGACTACGTCGTGATAGCTGCCATTTATCTTGCTGTGAAAACAGCTACCCAGTCCCTGAATTCCGTATGCACCTGCTTTATc                                                                                         <  2:496089/149‑1 (MQ=255)
                                                                  ttttCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATCCCGTATGCCCCTGCTTTATCTAACGGt                                                                                  >  2:449353/1‑149 (MQ=255)
                                                                  ttttCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATTCCCCTGCTTTATCTAACGGt                                                                                  >  1:262258/1‑149 (MQ=255)
                                                                  ttttCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGt                                                                                  >  1:473403/1‑149 (MQ=255)
                                                                  ttttCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGt                                                                                  >  1:1154021/1‑149 (MQ=255)
                                                                            cAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTACCGATTCATCGCTgg                                                                        >  2:553258/1‑149 (MQ=255)
                                                                                    tAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAg                                                                >  2:1254668/1‑149 (MQ=255)
                                                                                      aagttACTTAATAATTCATACGTTTCAACCAGCTGTGAAACCAGGACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCCCCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGcc                                                              <  2:912257/146‑1 (MQ=255)
                                                                                      aaaTTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGcc                                                              <  2:27639/149‑1 (MQ=255)
                                                                                         ttccTTAATCATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCACTTTGTCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcgc                                                           <  2:832108/146‑1 (MQ=255)
                                                                                         ttACTTAAAACTTCCCACGCTTCACCCACCGGTAAGCCGACTACGGCCTGATAACGTCCAGTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcgc                                                           <  1:1005586/149‑1 (MQ=255)
                                                                                         ctacttactAATTCATACTCTTCAACCCGCGGTCAGCCGACTACGGCGTGATAGCTGCCACTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcgc                                                           <  2:832374/148‑1 (MQ=255)
                                                                                           aCTTAATAATTCATACCTTTCAACCGGCGGTAAGCCGACTACGGCGTGATAGCTGCCGTTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGa                                                         <  2:104471/149‑1 (MQ=255)
                                                                                            cTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTTACAAAACAGCCTCCCGGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACTTAGCCCGCGAt                                                        >  1:83826/1‑149 (MQ=255)
                                                                                            cTTAATAATGCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGAt                                                        <  2:775842/149‑1 (MQ=255)
                                                                                             ttAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGGGACATAGCCCGCGagc                                                       >  2:842323/1‑147 (MQ=255)
                                                                                                 taatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCAGTTATCTTCCTGACAAAACAGTCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCtt                                                   >  1:45205/1‑149 (MQ=255)
                                                                                                 taatTCATACGCTTCACCCAGCGGTAAGCCGTCTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACGGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCtt                                                   <  1:578822/149‑1 (MQ=255)
                                                                                                    ttCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGCCTTCGt                                                >  1:806303/1‑149 (MQ=255)
                                                                                                         accctacAACCAGCGGTAACCCGCCTACGCCGTGATAGCTGCCGTTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTTGCGACATATCCCGCGATGTCTTCGTCTGtt                                           <  1:677589/143‑1 (MQ=255)
                                                                                                         aCGTTACAACCGGCGTTAACACGCCTACGGCCGAACACCTTCCATGTATTATCCTGACAAAACAGCCCCCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGtt                                           <  1:30848/149‑1 (MQ=255)
                                                                                                         aCGATTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAAACCGTATGCACCTGCTTTATCTAACGGTCCCTCGCTGGCGGCATAGCCACCGATGTCTCCGTCTgct                                           >  2:1275484/1‑147 (MQ=255)
                                                                                                           gTTTCACCCGGCGGTAATCCGCCTACGGCGTGATAGCGGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGCCATAGCCCGCGATGTCTTCGTCGGTTaa                                         <  2:909245/149‑1 (MQ=255)
                                                                                                            tttCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAAc                                        >  1:114463/1‑149 (MQ=255)
                                                                                                               cAACCAGCGGGAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACCAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGtt                                     <  2:746287/149‑1 (MQ=255)
                                                                                                               cAACCAGCGGCAAGCCGACTACGGCGTGATAGCTCCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGCATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGtt                                     <  1:524009/149‑1 (MQ=255)
                                                                                                                aaCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGACGTCTTCGTCTGTTAACGTTc                                    >  2:361432/1‑149 (MQ=255)
                                                                                                                  ccAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTg                                  >  2:827052/1‑149 (MQ=255)
                                                                                                                   cAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCCGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGGCTTCGCCTGGTAACGTTCtca                                 >  2:1038876/1‑147 (MQ=255)
                                                                                                                    aGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaa                                <  2:1180277/149‑1 (MQ=255)
                                                                                                                     ggccgacgCCCACCACGGAGTGCTCGCTGCCGTCCATCTTCCGTAGAATACAGCCCCTCCGCCCCGGAATGCCGGTTGCTCCTCCTTTATCTTACGGTCCATCGCCGGCGACATTTCCCGCGATTTCTGCGTCTGTTAACGTTCTGaaa                               <  2:521963/142‑1 (MQ=255)
                                                                                                                     gCGGTAAGCCGACTCCGGCGTTATAGCTGCCATTTATCTTCCGGCCAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaaa                               <  2:63845/149‑1 (MQ=255)
                                                                                                                     gCGGTAAGCCCAGTACGGCGCTATATCTGCCAGTTAGCTTCCTGACAAATCAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTCCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaaa                               <  2:647067/149‑1 (MQ=255)
                                                                                                                     acGGTAAGCCGAATACGGCGTGATAGCTGCCATTTATCTTCCGGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTGATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaaa                               <  2:801801/148‑1 (MQ=255)
                                                                                                                         tacgccgccTTCGGCGTGATCGCTCCCTTTTTTCTTCCTCACAAACCAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTca                           <  1:190018/146‑1 (MQ=255)
                                                                                                                         tAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTca                           >  1:710963/1‑149 (MQ=255)
                                                                                                                            gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATATAACGGTTCAAGACTAGCGACATAACCAGCGATGTCTTCGTCAGTTAACGGTCTGAAAGTAACAt                        >  2:232031/1‑149 (MQ=255)
                                                                                                                               gACTACGTCTTGCTCGCTGCTATTTCTCTCCCTCACAAACCAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCgg                     <  1:438734/149‑1 (MQ=255)
                                                                                                                               gACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCgg                     <  2:201917/149‑1 (MQ=255)
                                                                                                                                aCTACGTCGTGACAGCTGCCATCTACCTTCCCGCCATAACCGCCACCCTCCCCCTGAATACCCTATGCACCTCCTTTATCTCACGGTCCACCGTTGGCGACATCGCCCGCGATGTCTTCGTCTTTTAACGTTCTGAAAGTCACATCGGt                    <  2:677790/149‑1 (MQ=255)
                                                                                                                                aCTACGGCGTGATAGCTGCCATTTGTCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGt                    <  2:1265070/149‑1 (MQ=255)
                                                                                                                                  tACGGCTTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGa                  <  1:1177007/149‑1 (MQ=255)
                                                                                                                                   aCGGGGTGCTAGCTTCCATTTATCTTCCTGACAAAACAGCCCCCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac                 <  2:1064244/149‑1 (MQ=255)
                                                                                                                                   aCGGCGTTACAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac                 <  2:1063830/149‑1 (MQ=255)
                                                                                                                                   aCGGCGTGATAGCTGCCAGTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac                 <  2:1096047/149‑1 (MQ=255)
                                                                                                                                    cGGCGTGATAGCTGCCATTCATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacc                <  2:1147785/149‑1 (MQ=255)
                                                                                                                                     ggCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacca               >  1:455959/1‑149 (MQ=255)
                                                                                                                                          gATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAgg          >  1:849485/1‑149 (MQ=255)
                                                                                                                                                 gCCATTTCTCTTCCTGACAAACCAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCga   <  2:468642/149‑1 (MQ=255)
                                                                                                                                                  ccATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCgag  >  1:92121/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATATACGCCACTCGCGTTTCCGAAGGCGTTACGTTTACTAACAATTCAGCCGTCATGTTTATCCCTTTTCTCACGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAG  >  NC_000913/3397737‑3398031

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