Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,3630GC63.1% 83.3 / 135.4 103noncoding (605/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base C (41/24);  total (62/41)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

AGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTG  >  NC_000913/226215‑226442
                                                                                                                                                    |                                                                               
aGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGtt                                                                               <  2:976894/151‑1 (MQ=31)
aGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGtt                                                                               >  2:480357/1‑151 (MQ=31)
 gTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTa                                                                              <  1:334912/151‑1 (MQ=31)
  tACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTa                                                                              >  2:702895/1‑150 (MQ=31)
   aCCGTGAGGGGAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAAc                                                                            <  2:438755/150‑1 (MQ=11)
   aCCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCAACTTATATTCAGTGGCAAGCTTAAc                                                                            <  2:412182/150‑1 (MQ=11)
         aGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAg                                                                     <  1:380918/151‑1 (MQ=11)
           ggAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAgg                                                                    >  1:226132/1‑151 (MQ=31)
           ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAggg                                                                   >  1:177381/1‑151 (MQ=11)
           ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGATCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAggg                                                                   >  1:42341/1‑151 (MQ=11)
                  cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgt                                                            >  1:449045/1‑150 (MQ=11)
                  cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCg                                                             >  1:502222/1‑150 (MQ=11)
                   gAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGTAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGTCgt                                                            >  2:171468/1‑149 (MQ=11)
                    aaaaGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGa                                                            >  1:134803/1‑150 (MQ=31)
                    aaaaGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGa                                                            >  2:184176/1‑150 (MQ=31)
                          aCCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGaa                                                      >  1:463435/1‑150 (MQ=31)
                            cccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAAcc                                                   <  2:666489/151‑1 (MQ=31)
                             ccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                  >  1:608285/1‑151 (MQ=31)
                             ccGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGT‑GGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAg                                                <  2:412149/151‑1 (MQ=2)
                              cGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                  >  1:921466/1‑150 (MQ=31)
                              cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAg                                                <  1:460448/151‑1 (MQ=2)
                               ggagAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGGGAGt                                               >  1:413122/4‑151 (MQ=2)
                                  gAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCt                                             >  1:957184/1‑151 (MQ=31)
                                   aGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCtt                                            <  2:461930/151‑1 (MQ=31)
                                    ggggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTa                                           >  1:277448/1‑150 (MQ=2)
                                     gggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaat                                         >  1:368470/1‑150 (MQ=2)
                                      ggAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAAc                                         <  1:953607/151‑1 (MQ=31)
                                      ggAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAAc                                         >  1:839260/1‑151 (MQ=31)
                                       gAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaatag                                       <  2:48701/151‑4 (MQ=2)
                                       gAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACt                                        >  2:64389/1‑151 (MQ=31)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGGGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     >  1:268761/1‑146 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCg                                                             <  1:1997/127‑1 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCg                                                             >  2:1997/1‑127 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAg                                                <  2:582455/140‑1 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAg                                                >  1:582393/1‑140 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     <  2:402604/151‑6 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     >  1:102586/1‑146 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     >  2:161260/1‑146 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     >  1:500062/1‑146 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     >  2:97809/1‑146 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCATACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCg                                                             <  1:391361/127‑1 (MQ=2)
                                         gTGAAATAGAACCTGAAACCGTATGCATACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCg                                                             >  2:391395/1‑127 (MQ=2)
                                         gTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTgg                                      <  2:747198/151‑1 (MQ=31)
                                          tGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTggg                                     >  1:19965/1‑151 (MQ=31)
                                           gAAATATAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTATCGAAAGCGAGTCTTAATAGGGc                                    <  1:396295/150‑1 (MQ=2)
                                           gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggacg                                   <  2:449091/151‑8 (MQ=2)
                                           gAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTggg                                     <  1:571886/150‑1 (MQ=31)
                                           gAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTggg                                     <  1:274643/150‑1 (MQ=31)
                                            aaaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCg                                   >  2:972596/1‑151 (MQ=31)
                                             aataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGtt                                 >  2:92399/4‑151 (MQ=2)
                                              ataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGCAg                                                                >  2:45751/3‑119 (MQ=2)
                                              ataGAACCTGAAACCGGATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTGCCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGGGAGTCTTAATAGGGca                                   <  2:413158/146‑2 (MQ=2)
                                               taGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGttt                                >  1:122522/2‑149 (MQ=2)
                                                  aaCCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTATGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAg                              >  2:311113/1‑150 (MQ=14)
                                                   acttGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAAc                                         <  1:190034/134‑1 (MQ=2)
                                                   acttGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAAc                                         >  2:190045/4‑137 (MQ=2)
                                                    ccTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTg                           >  1:137449/1‑151 (MQ=31)
                                                    ccTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGc                          >  2:185312/1‑151 (MQ=2)
                                                     cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTTCCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAAcc                                                                           <  2:293683/101‑1 (MQ=2)
                                                     cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAAcc                                                                           >  1:293668/1‑101 (MQ=2)
                                                     cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGc                          >  1:678793/1‑150 (MQ=2)
                                                       gAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAAcc                                                                           <  1:328041/99‑1 (MQ=2)
                                                       gAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAAcc                                                                           >  2:328069/1‑99 (MQ=2)
                                                        aaaCCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    <  1:825302/137‑1 (MQ=2)
                                                        aaaCCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    >  2:825387/1‑137 (MQ=2)
                                                         aaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAggg                      >  1:19260/1‑151 (MQ=31)
                                                          aCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGt                     <  2:385631/151‑1 (MQ=31)
                                                            cgtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                   >  2:37337/7‑150 (MQ=2)
                                                            cGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                   >  2:47677/1‑151 (MQ=31)
                                                            cGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                   <  2:1028032/151‑1 (MQ=31)
                                                            cGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGta                    >  2:273044/1‑150 (MQ=31)
                                                              tgtACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                   >  1:994466/1‑149 (MQ=31)
                                                               gtacgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAg                 <  2:550086/150‑1 (MQ=31)
                                                               gtacgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGa                <  2:307001/151‑1 (MQ=31)
                                                                tgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaatagg                                      <  2:360449/125‑5 (MQ=2)
                                                                tgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaatagg                                      >  1:360417/3‑123 (MQ=2)
                                                                tgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGAcc              >  1:264104/3‑151 (MQ=2)
                                                                tgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGAcc              >  2:520531/3‑151 (MQ=2)
                                                                  cgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGTGTCTTAATATGGCGTTTAGTCGCTGGGTATAGACCCg            >  2:136539/1‑151 (MQ=2)
                                                                  cgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                   >  1:173032/1‑128 (MQ=2)
                                                                  cgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                   <  2:173041/128‑1 (MQ=2)
                                                                  cgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            >  1:340846/1‑151 (MQ=2)
                                                                  cgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            <  2:264119/151‑1 (MQ=2)
                                                                  cgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            >  1:406865/1‑151 (MQ=2)
                                                                  cgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            >  2:443267/1‑151 (MQ=2)
                                                                  cgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGAccc             <  2:409945/151‑1 (MQ=31)
                                                                   gtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    >  1:255754/1‑126 (MQ=2)
                                                                   gtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    <  2:255769/126‑1 (MQ=2)
                                                                      cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaatagg                                      >  1:841957/1‑117 (MQ=2)
                                                                      cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaatagg                                      <  2:842046/121‑5 (MQ=2)
                                                                      cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAc        >  1:222902/1‑151 (MQ=2)
                                                                      cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAc        >  1:408887/1‑151 (MQ=2)
                                                                      cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAACGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAc        >  2:301172/1‑151 (MQ=2)
                                                                      cAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGaaa         <  2:331664/151‑1 (MQ=31)
                                                                       aaGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    <  1:23775/122‑1 (MQ=2)
                                                                       aaGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    >  2:23777/1‑122 (MQ=2)
                                                                       aaGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAcc       >  1:420590/1‑151 (MQ=2)
                                                                        aGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAAcc       <  2:564397/151‑1 (MQ=31)
                                                                           aGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAAccc      <  1:277342/149‑1 (MQ=31)
                                                                             tGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTg  >  1:949943/1‑151 (MQ=31)
                                                                                                                                                    |                                                                               
AGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTG  >  NC_000913/226215‑226442

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: ‑