Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,424,646 C→A 5.2% noncoding (2138/2904 nt) rrlD ← 23S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,424,6460CA 5.2% 310.3 / 2.1 129noncoding (2138/2904 nt)rrlD23S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base C (71/48);  new base A (2/5);  total (73/56)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

CCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGACAGC  >  NC_000913/3424504‑3424795
                                                                                                                                              |                                                                                                                                                      
cctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTa                                                                                                                                                 >  2:1054498/1‑150 (MQ=31)
cctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACACGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTa                                                                                                                                                 >  2:480567/1‑150 (MQ=14)
  tccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCCACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCCCGGTCGGCTCCATGCAGCCTGGCGTCCCAACTTCTAAg                                                                                                                                              >  2:143905/1‑150 (MQ=11)
  tccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTAACCACCAGACACTGTCCGCAACCCGGAATACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACATCTAAg                                                                                                                                              >  2:458628/1‑150 (MQ=12)
   ccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGATCATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGc                                                                                                                                             <  1:2113360/150‑1 (MQ=17)
    cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGcc                                                                                                                                            >  2:1595100/1‑150 (MQ=32)
     gTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCt                                                                                                                                           >  1:826693/1‑150 (MQ=32)
        aCTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTccc                                                                                                                                        <  1:1883741/150‑1 (MQ=32)
        aCTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTccc                                                                                                                                        <  1:1218844/150‑1 (MQ=32)
        aCTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTccc                                                                                                                                        >  2:1994702/1‑150 (MQ=32)
         ctctTTAGGAGGAGACCGCCCGAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCa                                                                                                                                       <  1:256034/150‑1 (MQ=17)
         ctctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCa                                                                                                                                       <  2:1101706/150‑1 (MQ=32)
           ttcagggaggaGACCGCCCCAGTCAAACTACCCACCAGCCACTGTCCGCAACCCGGATTACGGGTCAACGTTATAACATCAAACATTAAAGGGTGGTATTTCAATGTCGGCTCCATGCAGACTGGCGTCCACACATCTAAGCCTCCCAcc                                                                                                                                     >  2:184196/6‑150 (MQ=11)
           ctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGTGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                                     <  2:1314269/150‑1 (MQ=17)
           ctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                                     <  2:1665223/150‑1 (MQ=32)
               aggaggAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGCTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTaaa                                                                                                                                 >  1:1246446/1‑148 (MQ=17)
                ggaggaGACCGTCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCCCCGATcc                                                                                                                                >  2:1921692/1‑150 (MQ=12)
                ggaggaGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATcc                                                                                                                                >  1:1321776/1‑150 (MQ=32)
                   ggagACCGCCCCAGCCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTac                                                                                                                             <  2:1359873/150‑1 (MQ=17)
                    gagaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTaca                                                                                                                            <  1:1307429/150‑1 (MQ=32)
                      gaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTAcaca                                                                                                                          >  2:1648387/1‑150 (MQ=32)
                      gaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTAcaca                                                                                                                          >  1:2018672/1‑150 (MQ=32)
                         cGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGCGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCa                                                                                                                       >  1:1897157/1‑150 (MQ=17)
                          gCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCACACAGTAAAGGGTGGTATTTCTAGGTCGGCTCCATGCAGAATGCCGTCCACACTTCTAAGCCGCCCAAATATCCTACACATCaa                                                                                                                      >  2:1338105/1‑150 (MQ=11)
                          gCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACGCATCaa                                                                                                                      >  1:1433878/1‑150 (MQ=17)
                           ccccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAg                                                                                                                     >  2:1085618/1‑150 (MQ=32)
                            cccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAgg                                                                                                                    <  1:633907/150‑1 (MQ=32)
                               aGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACATCTAAGCCTCCCACCTATCCTACACAACAACGCTc                                                                                                                 >  2:424458/1‑150 (MQ=12)
                               aGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACAGCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACCCTTCTAAGCCTCCCACCTATCCTACCCATCAAGGCTc                                                                                                                 >  1:554028/1‑150 (MQ=12)
                                gTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCa                                                                                                                <  1:1994702/150‑1 (MQ=32)
                                   aaaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGTGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATg                                                                                                             <  1:940112/150‑1 (MQ=17)
                                   aaaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATg                                                                                                             <  1:940116/150‑1 (MQ=32)
                                      cTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTc                                                                                                          <  1:1842146/150‑1 (MQ=32)
                                      cTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTc                                                                                                          >  2:1732670/1‑150 (MQ=32)
                                        aCCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCCCCTATCCTACCCATCAAGGCTCAATGTTCAg                                                                                                        >  1:1012426/1‑150 (MQ=17)
                                         cccAGCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTTGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAgt                                                                                                       <  1:1689291/150‑1 (MQ=17)
                                          ccaccaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAgtg                                                                                                      >  2:772995/1‑150 (MQ=32)
                                            accaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCCTGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTc                                                                                                    >  1:1817574/1‑150 (MQ=17)
                                            accaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTTCACACTTCTAAGCCTCCCACCTTTGCTACACATCAAGGCTCAATGTTCAGTGTc                                                                                                    <  2:876179/150‑1 (MQ=12)
                                            accaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTc                                                                                                    >  1:661651/1‑150 (MQ=32)
                                            accaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGCCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTc                                                                                                    <  1:842108/150‑1 (MQ=17)
                                                 acacTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCt                                                                                               >  1:2143976/1‑150 (MQ=32)
                                                   acTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCtat                                                                                             >  2:210373/1‑150 (MQ=32)
                                                   acTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCtat                                                                                             >  1:1343761/1‑150 (MQ=32)
                                                     tGTCCGCAACCCGGATTACGGGTCAACGTTGGAACATCAAACATTAAAGGGTGGTATTTTAAGGTCGTATTAATTATTTTTTTTTTTTTAATTTTAATGCGTCCAACGTTTTCTACACATTAAGGCTCAATGTTCAGTGTCAAGCTATAg                                                                                           <  2:306141/150‑1 (MQ=0)
                                                      gTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGt                                                                                          >  1:1124173/1‑150 (MQ=32)
                                                       tCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGCGGCGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTa                                                                                         <  1:375979/150‑1 (MQ=12)
                                                        ccGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTaa                                                                                        >  2:1291636/1‑150 (MQ=32)
                                                        ccGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTaa                                                                                        >  1:670627/1‑150 (MQ=32)
                                                            aaCCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCAAACTTCTAAGCCTCACACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGt                                                                                    >  1:1414387/1‑150 (MQ=17)
                                                                          gTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTc                                                                      >  2:1359766/1‑150 (MQ=32)
                                                                          gTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTc                                                                      <  1:454118/150‑1 (MQ=32)
                                                                           tCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTcc                                                                     >  2:1216798/1‑150 (MQ=32)
                                                                           tCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTcc                                                                     >  2:1456151/1‑150 (MQ=32)
                                                                              aCGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTc                                                                  <  1:1239948/150‑1 (MQ=32)
                                                                               cGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCt                                                                 <  1:2109979/150‑1 (MQ=32)
                                                                                 ttAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGTCTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGATCAATTTTCAGTGTCAAGCTATAGTAAAGGTTCAAGGGTTCTTTCCTTCTTg                                                               >  2:1187171/1‑150 (MQ=11)
                                                                                 ttAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTg                                                               <  1:424458/150‑1 (MQ=32)
                                                                                   aGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGcc                                                             >  2:708611/1‑150 (MQ=32)
                                                                                   aGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGcc                                                             <  1:2018598/150‑1 (MQ=32)
                                                                                    gAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcg                                                            >  1:1871691/1‑150 (MQ=32)
                                                                                      aCATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcgcg                                                          >  2:1215928/1‑150 (MQ=32)
                                                                                       cATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCgg                                                         >  2:215752/1‑150 (MQ=32)
                                                                                        aTCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATTCAGATTTTTTTTTTATTGTTAATGCGTCCAACGATTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCggg                                                        <  2:1528986/150‑1 (MQ=0)
                                                                                          cAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTa                                                      >  1:1187879/1‑150 (MQ=32)
                                                                                          cAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTa                                                      >  2:101646/1‑150 (MQ=32)
                                                                                          cAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGAGGTTa                                                      >  2:1540534/1‑150 (MQ=17)
                                                                                           aaaCATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTac                                                     >  2:2174330/1‑150 (MQ=32)
                                                                                              cATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCTGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACt                                                  <  2:569686/150‑1 (MQ=17)
                                                                                                ttAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGGCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGc                                                <  1:593345/150‑1 (MQ=17)
                                                                                                 tAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGTATTCCACAATTCTAAGCGTCCAACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCa                                               <  2:1730220/150‑1 (MQ=11)
                                                                                                 tAAAGGGTCGTATTTCAAGGTCGGCTCCATGCATACTTGTTTTTAAAATTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCa                                               <  2:1096579/150‑1 (MQ=11)
                                                                                                  aaaGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCAt                                              >  2:1586464/1‑150 (MQ=32)
                                                                                                    aGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCt                                            <  1:1648387/150‑1 (MQ=32)
                                                                                                     gggTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCtt                                           >  2:1417008/1‑150 (MQ=32)
                                                                                                        gggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTcac                                        <  1:1254916/149‑1 (MQ=17)
                                                                                                         ggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTcaca                                       <  1:1437797/150‑1 (MQ=32)
                                                                                                         ggtATTTCAAGGTCGGATCCATGCCGACTGTCGTCCACACTTCTAAGCATCCCACCTATAATACACAAAAAGGCTCAATGTTCAGTGTCAAGCTATAGTTAAGGTTCCCGGGGTCTTTCAGTCTTTCCGCGGGTACACTGCATCTTcaca                                       >  2:2062297/1‑150 (MQ=11)
                                                                                                         ggtATATCAAGGTCGGCTCCATGCAGTTTGTTTTTATAATTTTAATGCGTCCAACGTTATCTACACATAAAGGCTCTATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTcaca                                       <  2:424510/150‑1 (MQ=0)
                                                                                                           tATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                     >  1:965214/1‑150 (MQ=32)
                                                                                                            aTTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAATCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCg                                    <  2:829612/150‑1 (MQ=17)
                                                                                                            aTTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCg                                    <  1:101646/150‑1 (MQ=32)
                                                                                                              ttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAg                                  <  1:708611/150‑1 (MQ=32)
                                                                                                                cAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTCCACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAggg                                >  1:1103307/1‑148 (MQ=17)
                                                                                                                cAAGGTCGGCTCCATGCAGAATGGTTTTTAAATTTTTAAGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGtt                                <  2:2057007/150‑1 (MQ=11)
                                                                                                                  aGGTCGGCTCCATGCAGACTTTTTTTTAAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCa                              <  2:704920/150‑1 (MQ=11)
                                                                                                                    gTCGGCTCCATGCAGACTGGCGTCCACACTTCTTAGCCTCCCACCTATCCTACAACTCCAGGCTCAATGTTCAGTGTCAATCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAAt                            <  1:208408/150‑1 (MQ=11)
                                                                                                                    gTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAAt                            >  2:211792/1‑150 (MQ=32)
                                                                                                                       ggCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGTCTCAATGTTCAGTGTCAAGCTGTAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCAGCGGGGACACTGCATCTTCACAGCGAGTTCAAttta                         >  2:275561/1‑149 (MQ=255)
                                                                                                                       ggCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCACACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTc                         <  2:1124173/150‑1 (MQ=17)
                                                                                                                         cTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCAc                       >  1:2082737/1‑150 (MQ=32)
                                                                                                                              tGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGt                  <  1:1054498/150‑1 (MQ=32)
                                                                                                                                      ggCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTg          >  2:1623864/1‑150 (MQ=32)
                                                                                                                                         gTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCg              <  1:2055244/143‑1 (MQ=32)
                                                                                                                                         gTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCg              >  2:2055244/1‑143 (MQ=32)
                                                                                                                                         gTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGgag       >  1:1779276/1‑150 (MQ=32)
                                                                                                                                          tCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGgaga      >  2:665128/1‑150 (MQ=32)
                                                                                                                                          tCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGgaga      <  1:1590112/150‑1 (MQ=32)
                                                                                                                                              cacTTCTATGCCTCCCACCTATACTACACATCAAGGCTCAATGTTGAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGGGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGACAGc  <  2:347017/150‑1 (MQ=255)
                                                                                                                                              cacTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCGATAGTAAAGGTGCACGGGGTCGTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGGTCTATTTCACTGAGTCTCGGGTGgagcctgc  >  2:127727/1‑145 (MQ=11)
                                                                                                                                              |                                                                                                                                                      
CCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGACAGC  >  NC_000913/3424504‑3424795

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