Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,5550GC69.6% 73.7 / 80.9 79noncoding (797/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base G (8/16);  new base C (31/24);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 8.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

GGCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACAC  >  NC_000913/226454‑226704
                                                                                                     |                                                                                                                                                      
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:70505/151‑1 (MQ=31)
                             ggaggaCCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCgg                                                                         >  1:702472/1‑151 (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                                                                         <  1:383960/116‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCgg                                                                         >  2:383920/1‑116 (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                                      >  2:173217/1‑150 (MQ=1)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCaa                                      >  1:724366/1‑150 (MQ=1)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTAGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAggg                                                                   >  1:334698/1‑121 (MQ=0)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTAGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAggg                                                                   <  2:334662/121‑1 (MQ=0)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgggggcc                                                               >  2:398932/1‑122 (MQ=0)
                                                                 aTGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgggggcc                                                               <  1:398972/124‑3 (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                                                                        >  2:123981/1‑111 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                                        <  1:123993/111‑1 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgg                                                                    <  2:32326/115‑1 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgg                                                                    >  1:32331/1‑115 (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                                              <  1:1174/124‑3 (MQ=0)
                                                                           gatcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAGGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCaaac                                              >  2:1174/8‑129 (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                                                                    <  1:51353/102‑1 (MQ=0)
                                                                            atcggagTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGAACTGTTTCGGCTAgg                                                                    >  2:51347/7‑108 (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                                                                   <  1:56290/102‑1 (MQ=0)
                                                                                   tGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCCAggg                                                                   >  2:56283/1‑102 (MQ=0)
                                                                                    gAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATg                  <  2:52124/151‑1 (MQ=31)
                                                                                       aGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagtgccga              <  1:230628/151‑13 (MQ=0)
                                                                                        ggCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagtgccgag             <  1:343354/151‑14 (MQ=0)
                                                                                        ggCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCAAACTCCGAATAcccagaagtgccgag             <  2:120641/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)
                                                                                                     |                                                                                                                                                      
GGCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACAC  >  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: