Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,5740CT69.6% 84.1 / 81.5 92noncoding (816/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base C (10/18);  new base T (34/30);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GGCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTC  >  NC_000913/226454‑226722
                                                                                                                        |                                                                                                                                                       
ggCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAAc                                                                                                                           <  1:603770/151‑1 (MQ=11)
        gAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTcc                                                                                                                  >  1:366758/1‑151 (MQ=31)
              taggTAACACTGACTGGAGGACCGAACCCACTCCCGTTGAAAAGGTAGGGGATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATgtatcactctggg                                                                                                                  <  2:209539/144‑13 (MQ=0)
              taggTAACACTGACTGGAGGACCGAACCCACTCCCGTTGAAAAGGTAGGGGATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATgtatcactctggg                                                                                                                  >  1:209556/3‑134 (MQ=0)
               gggTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGt                                                                                                            <  2:274170/150‑1 (MQ=31)
                ggTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTag                                                                                                          <  1:442925/151‑1 (MQ=31)
                   aaCACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGc                                                                                                       >  1:687420/1‑151 (MQ=31)
                     cacTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCAc                                                                                                     <  2:176382/151‑1 (MQ=31)
                          aCTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGttt                                                                                                >  2:300773/1‑151 (MQ=31)
                           cTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTc                                                                                               >  2:37725/1‑151 (MQ=31)
                             ggaggaCCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTc                                                                                               <  2:729531/149‑1 (MQ=31)
                             ggaggaCCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCgg                                                                                             >  1:702472/1‑151 (MQ=31)
                             ggaggaCCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCgg                                                                                             <  1:70505/151‑1 (MQ=31)
                                  accgaaccgaCTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAgg                                                                                        <  2:699605/151‑1 (MQ=31)
                                  accgaacccaCTCCCTTTGAAAAGGTAGGGGATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATgtatcactctg                                                                                                                    >  1:66878/1‑114 (MQ=0)
                                  accgaacccaCTCCCTTTGAAAAGGTAGGGGATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATgtatcactctg                                                                                                                    <  2:66871/124‑11 (MQ=0)
                                   ccgaacccaCTCCCGTTGAAAAGGTAGGGGATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATgtatcactctggg                                                                                                                  <  2:361341/125‑13 (MQ=0)
                                   ccgaacccaCTCCCGTTGAAAAGGTAGGGGATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATgtatcactctggg                                                                                                                  >  1:361379/1‑113 (MQ=0)
                                             tAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCg                                                                             >  1:685312/1‑151 (MQ=31)
                                                       aaaTTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAAc                                                                   <  2:720420/151‑1 (MQ=31)
                                                            aGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGAt                                                              <  2:783043/151‑1 (MQ=31)
                                                            aGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGAt                                                              >  2:44732/1‑151 (MQ=31)
                                                              gggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                                                            >  2:60966/2‑118 (MQ=0)
                                                              gggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                                                            <  1:60973/117‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGt                                                                                                  >  2:312736/1‑111 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGt                                                                                                  <  1:312767/111‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCgg                                                                                             >  2:383920/1‑116 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCgg                                                                                             <  1:383960/116‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                                                            <  1:242630/117‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                                                            >  2:242605/1‑117 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCa                                                           >  1:110064/1‑150 (MQ=1)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCa                                                           >  2:230610/1‑150 (MQ=1)
                                                               ggATAACTTGTGGAACGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCa                                                           >  1:282477/1‑150 (MQ=0)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGct                                                                                           <  1:144954/117‑2 (MQ=0)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCaa                                                          >  1:724366/1‑150 (MQ=1)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCaa                                                          >  2:173217/1‑150 (MQ=1)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTAGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAggg                                                                                       <  2:334662/121‑1 (MQ=0)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTAGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAggg                                                                                       >  1:334698/1‑121 (MQ=0)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgggggcc                                                                                   <  1:398972/124‑3 (MQ=0)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgggggcc                                                                                   >  2:398932/1‑122 (MQ=0)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAAc                                                        <  1:173230/151‑1 (MQ=1)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTACTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAAc                                                        >  2:275068/1‑151 (MQ=1)
                                                                   gACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACtc                                                      >  1:391795/1‑150 (MQ=1)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGcta                                                                                          >  1:352892/1‑111 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGcta                                                                                          <  2:352854/113‑3 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                                                            <  1:123993/111‑1 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                                                            >  2:123981/1‑111 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgg                                                                                        >  1:32331/1‑115 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgg                                                                                        <  2:32326/115‑1 (MQ=0)
                                                                       tgtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCgg                                                                                             <  2:354666/108‑1 (MQ=0)
                                                                       tgtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCgg                                                                                             >  1:354704/1‑108 (MQ=0)
                                                                       tgtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGa                                                   >  2:137510/1‑150 (MQ=1)
                                                                        gtgGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGa                                                   <  1:526095/150‑1 (MQ=31)
                                                                        gtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGaa                                                  >  1:255553/1‑150 (MQ=1)
                                                                        gtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAAt                                                 >  2:343317/1‑151 (MQ=1)
                                                                          ggatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTACCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAc                                               >  1:775449/9‑151 (MQ=11)
                                                                           gatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAGGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCaaac                                                                  >  2:1174/8‑129 (MQ=0)
                                                                           gatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAGGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCaaac                                                                  <  1:1174/124‑3 (MQ=0)
                                                                           gatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCgg                                                                                             >  1:141159/8‑104 (MQ=0)
                                                                           gatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCgg                                                                                             <  2:141147/97‑1 (MQ=0)
                                                                           gCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATAc                                               <  2:189193/151‑1 (MQ=31)
                                                                            atcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGAGGTCATCCCGACTTACCAAacc                                                                 >  2:319663/7‑131 (MQ=0)
                                                                            atcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGAGGTCATCCCGACTTACCAAacc                                                                 <  1:319696/125‑1 (MQ=0)
                                                                            atcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGAACTGTTTCGGCTAgg                                                                                        >  2:51347/7‑108 (MQ=0)
                                                                            atcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGAACTGTTTCGGCTAgg                                                                                        <  1:51353/102‑1 (MQ=0)
                                                                              gggggTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCg                                             <  1:529172/150‑1 (MQ=31)
                                                                              cggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccag                                           <  1:137522/147‑4 (MQ=1)
                                                                               ggggTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCgg                                            >  2:295880/1‑150 (MQ=31)
                                                                                  gTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagt                                       >  1:120653/1‑144 (MQ=1)
                                                                                  gTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaag                                        >  1:162012/1‑144 (MQ=1)
                                                                                   tGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCCAggg                                                                                       >  2:56283/1‑102 (MQ=0)
                                                                                   tGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCCAggg                                                                                       <  1:56290/102‑1 (MQ=0)
                                                                                    gAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATg                                     <  2:52124/151‑1 (MQ=31)
                                                                                       aGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagtgccga                                  <  1:230628/151‑13 (MQ=0)
                                                                                        ggCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagtgccgag                                 <  2:120641/151‑14 (MQ=0)
                                                                                        ggCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagtgccgag                                 <  1:343354/151‑14 (MQ=0)
                                                                                                   ccGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGaca                     <  2:216082/151‑1 (MQ=31)
                                                                                                   ccGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGaca                     <  2:264072/151‑1 (MQ=31)
                                                                                                    cGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGacac                    <  2:28914/151‑1 (MQ=31)
                                                                                                      ggAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGacacac                  >  1:736034/1‑151 (MQ=31)
                                                                                                       gagaTAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACAcg                 <  2:350340/151‑1 (MQ=31)
                                                                                                            aGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTACCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATACCCAGAAGTGCCGAGCATGGGAGACACAcggcgg             <  2:775367/150‑1 (MQ=2)
                                                                                                            aGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATACCCAGAAGTGCCGAGCATGGGAGACACAcggcgg             <  1:188242/150‑1 (MQ=0)
                                                                                                            aGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATACCCAGAAGTGCCGAGCATGGGAGACACAcggcgg             <  2:110053/150‑1 (MQ=0)
                                                                                                            aGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATACCCAGAAGTGCCGAGCATGGGAGACACAcggcgg             <  2:282449/150‑1 (MQ=0)
                                                                                                            aGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATACCCAGAAGTGCCGAGCATGGGAGACACACGGCggg            <  2:724293/151‑1 (MQ=0)
                                                                                                            aGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATACCCAGAAGTGCCGAGCACGGGAGACACACGGCggg            <  2:255527/151‑1 (MQ=0)
                                                                                                             gCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATACCCAGAAGTGCCGAGCATGGGAGACACACGGCggg            >  2:295485/1‑150 (MQ=0)
                                                                                                             gCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCggg            >  1:695421/1‑150 (MQ=31)
                                                                                                             gCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGt           <  2:721155/151‑1 (MQ=31)
                                                                                                              cTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATACCCAGAAGTGCCGAGCATGGGAGACACACGGCGGGTg          >  1:365667/1‑151 (MQ=0)
                                                                                                                    atcCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAGTACCCAGAAGTGCCGAGCATGGGAGACACACGGCGGGTGCTAAcg    >  2:309347/2‑151 (MQ=0)
                                                                                                                      ccccGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAAcgtc  <  1:154779/151‑1 (MQ=31)
                                                                                                                      ccccGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAAcgtc  <  1:446578/151‑1 (MQ=31)
                                                                                                                      ccccGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAAcgtc  <  2:395903/151‑1 (MQ=31)
                                                                                                                        |                                                                                                                                                       
GGCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTC  >  NC_000913/226454‑226722

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