Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,5540CT72.4% 57.9 / 56.4 76noncoding (796/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base T (31/24);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GGCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACAC  >  NC_000913/226454‑226704
                                                                                                    |                                                                                                                                                       
ggCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAAc                                                                                                       <  1:603770/151‑1 (MQ=11)
        gAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTc‑c                                                                                              >  1:366758/1‑151 (MQ=31)
              taggTAACACTGACTGGAGGACCGAACCCACTCCCGTTGAAAAGGTAGGGGATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATgtatcactctggg                                                                                              <  2:209539/144‑13 (MQ=0)
              taggTAACACTGACTGGAGGACCGAACCCACTCCCGTTGAAAAGGTAGGGGATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATgtatcactctggg                                                                                              >  1:209556/3‑134 (MQ=0)
               gggTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGt                                                                                        <  2:274170/150‑1 (MQ=31)
                ggTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTag                                                                                      <  1:442925/151‑1 (MQ=31)
                   aaCACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGc                                                                                   >  1:687420/1‑151 (MQ=31)
                     cacTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCAc                                                                                 <  2:176382/151‑1 (MQ=31)
                          aCTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGttt                                                                            >  2:300773/1‑151 (MQ=31)
                           cTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTc                                                                           >  2:37725/1‑151 (MQ=31)
                             ggaggaCCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTc                                                                           <  2:729531/149‑1 (MQ=31)
                             ggaggaCCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCgg                                                                         <  1:70505/151‑1 (MQ=31)
                             ggaggaCCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCgg                                                                         >  1:702472/1‑151 (MQ=31)
                                  accgaaccgaCTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAgg                                                                    <  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)
                                             tAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCg                                                         >  1:685312/1‑151 (MQ=31)
                                                       aaaTTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAAc                                               <  2:720420/151‑1 (MQ=31)
                                                            aGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGAt                                          <  2:783043/151‑1 (MQ=31)
                                                            aGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGAt                                          >  2:44732/1‑151 (MQ=31)
                                                              gggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGc                                                                        >  2:60966/2‑118 (MQ=0)
                                                              gggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGc                                                                        <  1:60973/117‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGt                                                                              >  2:312736/1‑111 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGt                                                                              <  1:312767/111‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCgg                                                                         <  1:383960/116‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCgg                                                                         >  2:383920/1‑116 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGc                                                                        <  1:242630/117‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGc                                                                        >  2:242605/1‑117 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCa                                       >  1:110064/1‑150 (MQ=1)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCa                                       >  2:230610/1‑150 (MQ=1)
                                                               ggATAACTTGTGGAACGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCa                                       >  1:282477/1‑150 (MQ=0)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGct                                                                       <  1:144954/117‑2 (MQ=0)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCaa                                      >  2:173217/1‑150 (MQ=1)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCaa                                      >  1:724366/1‑150 (MQ=1)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTAGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAggg                                                                   >  1:334698/1‑121 (MQ=0)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTAGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAggg                                                                   <  2:334662/121‑1 (MQ=0)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAgggggcc                                                               >  2:398932/1‑122 (MQ=0)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAgggggcc                                                               <  1:398972/124‑3 (MQ=0)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAAc                                    <  1:173230/151‑1 (MQ=1)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTACTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAAc                                    >  2:275068/1‑151 (MQ=1)
                                                                   gACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACtc                                  >  1:391795/1‑150 (MQ=1)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGcta                                                                      >  1:352892/1‑111 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGcta                                                                      <  2:352854/113‑3 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGc                                                                        >  2:123981/1‑111 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGc                                                                        <  1:123993/111‑1 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAgg                                                                    <  2:32326/115‑1 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAgg                                                                    >  1:32331/1‑115 (MQ=0)
                                                                       tgtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCgg                                                                         <  2:354666/108‑1 (MQ=0)
                                                                       tgtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCgg                                                                         >  1:354704/1‑108 (MQ=0)
                                                                       tgtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGa                               >  2:137510/1‑150 (MQ=1)
                                                                        gtgGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGa                               <  1:526095/150‑1 (MQ=31)
                                                                        gtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGaa                              >  1:255553/1‑150 (MQ=1)
                                                                        gtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAAt                             >  2:343317/1‑151 (MQ=1)
                                                                          ggatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTACCACTCT‑‑GGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAc                           >  1:775449/9‑151 (MQ=11)
                                                                           gatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAGGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCaaac                                              <  1:1174/124‑3 (MQ=0)
                                                                           gatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAGGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCaaac                                              >  2:1174/8‑129 (MQ=0)
                                                                           gatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCgg                                                                         >  1:141159/8‑104 (MQ=0)
                                                                           gatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCgg                                                                         <  2:141147/97‑1 (MQ=0)
                                                                           gCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATAc                           <  2:189193/151‑1 (MQ=31)
                                                                            atcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGAGGTCATCCCGACTTACCAAacc                                             >  2:319663/7‑131 (MQ=0)
                                                                            atcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGAGGTCATCCCGACTTACCAAacc                                             <  1:319696/125‑1 (MQ=0)
                                                                            atcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGAACTGTTTCGGCTAgg                                                                    <  1:51353/102‑1 (MQ=0)
                                                                            atcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGAACTGTTTCGGCTAgg                                                                    >  2:51347/7‑108 (MQ=0)
                                                                              gggggTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCg                         <  1:529172/150‑1 (MQ=31)
                                                                              cggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccag                       <  1:137522/147‑4 (MQ=1)
                                                                               ggggTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCgg                        >  2:295880/1‑150 (MQ=31)
                                                                                  gTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagt                   >  1:120653/1‑144 (MQ=1)
                                                                                  gTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaag                    >  1:162012/1‑144 (MQ=1)
                                                                                   tGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCCAggg                                                                   <  1:56290/102‑1 (MQ=0)
                                                                                   tGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCCAggg                                                                   >  2:56283/1‑102 (MQ=0)
                                                                                    gAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATg                  <  2:52124/151‑1 (MQ=31)
                                                                                       aGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagtgccga              <  1:230628/151‑13 (MQ=0)
                                                                                        ggCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagtgccgag             <  1:343354/151‑14 (MQ=0)
                                                                                        ggCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTA‑TCA‑CTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagtgccgag             <  2:120641/151‑14 (MQ=0)
                                                                                                   ccGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGaca   <  2:216082/151‑1 (MQ=31)
                                                                                                   ccGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGaca   <  2:264072/151‑1 (MQ=31)
                                                                                                    cGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGacac  <  2:28914/151‑1 (MQ=31)
                                                                                                    |                                                                                                                                                       
GGCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTC‑CGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACAC  >  NC_000913/226454‑226704

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