Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,5450CT72.4% 62.1 / 65.5 76noncoding (787/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 = 1.25e-01

GGCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATC  >  NC_000913/226454‑226693
                                                                                           |                                                                                                                                                     
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:44732/1‑151 (MQ=31)
                                                            aGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGAt                               <  2:783043/151‑1 (MQ=31)
                                                              gggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                             <  1:60973/117‑1 (MQ=0)
                                                              gggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                             >  2:60966/2‑118 (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                                                             >  2:242605/1‑117 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                             <  1:242630/117‑1 (MQ=0)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCa                            >  2:230610/1‑150 (MQ=1)
                                                               ggATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAGGGGGTCATCCCGACTTACCAAACCGATGCa                            >  1:110064/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                                                        <  2:334662/121‑1 (MQ=0)
                                                                gATGACTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTAGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAggg                                                        >  1:334698/1‑121 (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                                                           <  2:352854/113‑3 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGcta                                                           >  1:352892/1‑111 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                             <  1:123993/111‑1 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGc                                                             >  2:123981/1‑111 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgg                                                         <  2:32326/115‑1 (MQ=0)
                                                                     cTTGTGGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCGGCTAgg                                                         >  1:32331/1‑115 (MQ=0)
                                                                       tgtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCgg                                                              >  1:354704/1‑108 (MQ=0)
                                                                       tgtgGATCGGAGTGAAAGGCTAATCAAGCTCGGAGATAGCTGGTTCTCCTCGAAAGCTATTTAGGTAGCGCCTCATGTATCACTCTGGGGGGTAGAGCACTGTTTCgg                                                              <  2:354666/108‑1 (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                                                        <  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)
                                                                                           |                                                                                                                                                     
GGCAGGTTGAAGGTTGGGTAACACTAACTGGAGGACCGAACCGACTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATC  >  NC_000913/226454‑226693

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: