Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,3010GA67.0% 60.9 / 244.8 221noncoding (543/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base G (42/31);  new base A (84/64);  total (126/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 222 total aligned reads displayed.

CCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTG  >  NC_000913/226151‑226442
                                                                                                                                                      |                                                                                                                                             
ccTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATg                                                                                                                                               <  1:445576/151‑1 (MQ=31)
    tCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCtt                                                                                                                                            >  1:663402/1‑150 (MQ=31)
     cTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCtt                                                                                                                                            >  1:162918/1‑149 (MQ=31)
       gAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAgg                                                                                                                                         <  2:289401/150‑1 (MQ=31)
       gAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAg                                                                                                                                          <  1:588457/149‑1 (MQ=31)
          tatGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCgtg                                                                                                                                     >  1:152308/1‑151 (MQ=31)
                  gACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCg                                                                                                                             >  1:953425/1‑151 (MQ=31)
                          tccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTg                                                                                                                    >  1:145477/1‑151 (MQ=11)
                               ggCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGtat                                                                                                                 >  2:305310/1‑150 (MQ=31)
                                 cTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAAt                                                                                                              >  2:759384/1‑151 (MQ=31)
                                 cTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATg                                                                                                             >  2:208914/1‑151 (MQ=11)
                                   aaaTACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATg                                                                                                             <  1:37887/149‑1 (MQ=11)
                                    aaTACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATggg                                                                                                           >  1:98831/1‑151 (MQ=31)
                                        ctactgactgaCCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGc                                                                                                      >  2:64417/4‑151 (MQ=11)
                                             gactgacCGATAGTGAACTAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGtata                                                                                                                >  2:27888/1‑136 (MQ=11)
                                             gactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     >  2:46592/1‑147 (MQ=11)
                                              actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACtt                                                                                                 >  1:274192/1‑151 (MQ=31)
                                               ctgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACGCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAAT‑GGTCAGCGACTtat                                                                                               >  1:442505/1‑150 (MQ=11)
                                                    cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACGCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAAT‑GGTCAGCGACTTATATTcag                                                                                         <  2:442553/151‑3 (MQ=11)
                                                    cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTc                                                                                           >  1:632530/1‑151 (MQ=31)
                                                    cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGtat                                                                                                                 <  1:198270/128‑1 (MQ=11)
                                                     ggtAGTGGACCAGTACCCTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         <  2:347696/149‑3 (MQ=11)
                                                     gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATACTcag                                                                                         <  1:208899/151‑3 (MQ=11)
                                                     gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGA‑GGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        <  2:168950/151‑1 (MQ=11)
                                                         gTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAg                                                                                      <  2:308144/151‑1 (MQ=31)
                                                         gTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAATCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGtgg                                                                                      <  1:517567/150‑3 (MQ=11)
                                                         gTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATAttaagtggc                                                                                     >  1:33312/1‑144 (MQ=11)
                                                            aaCCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAg                                                                                  <  2:145488/151‑1 (MQ=11)
                                                            aaCCAGTACAGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGcg                                                                                    <  2:177088/149‑2 (MQ=11)
                                                             aCCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAg                                                                                  <  2:119309/151‑1 (MQ=31)
                                                              ggaGTACAGTGAGGGAAAGGAGAAAAGAAACCAGGAGAGGGGAGTGAAATAGAACCTGAAACAGTATGAGTACAAGCAGTGGGAGCCTAGTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                     >  1:102381/3‑146 (MQ=2)
                                                              ccAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCa                                                                                    <  2:102389/147‑1 (MQ=11)
                                                               cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGt                                                                                >  2:544075/1‑151 (MQ=31)
                                                               cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        >  1:232556/1‑142 (MQ=11)
                                                                   aCCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTaa                                                                             <  1:788509/150‑1 (MQ=31)
                                                                    ccGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAAcc                                                                           >  1:62081/1‑151 (MQ=31)
                                                                      gTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCg                                                                          >  2:684611/1‑150 (MQ=31)
                                                                      gTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGa                                                                         <  1:312514/151‑1 (MQ=31)
                                                                       tGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGa                                                                         >  2:526108/1‑150 (MQ=31)
                                                                        gAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGaa                                                                        <  2:572705/150‑1 (MQ=31)
                                                                         gaggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                     >  1:375634/3‑135 (MQ=11)
                                                                           ggAAAGGCTAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                                                         >  1:489068/1‑113 (MQ=11)
                                                                           ggAAAGGCTAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                                                         <  2:489120/113‑1 (MQ=11)
                                                                           ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAg                                                                                  <  1:82257/136‑1 (MQ=11)
                                                                             aaaGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAgggg                                                                  >  1:647575/1‑151 (MQ=31)
                                                                             aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAg                                                                                                       >  2:244893/1‑113 (MQ=11)
                                                                             aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGAc                                                                                                   <  1:256556/117‑1 (MQ=11)
                                                                               aGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAgggg                                                                  >  2:691459/1‑149 (MQ=31)
                                                                                 gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATgggc                                                                                                          >  2:493902/1‑105 (MQ=11)
                                                                                 gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATgggc                                                                                                          <  1:493849/106‑2 (MQ=11)
                                                                                  cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATg                                                                                                             >  2:36887/1‑102 (MQ=11)
                                                                                  cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         >  1:217684/1‑120 (MQ=11)
                                                                                  cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         <  2:217700/122‑3 (MQ=11)
                                                                                    aaaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTTGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     >  1:453750/1‑108 (MQ=11)
                                                                                    aaaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTTGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     <  2:453801/108‑1 (MQ=11)
                                                                                             ccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGaaa                                                     <  1:894013/148‑1 (MQ=31)
                                                                                              cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         <  1:351640/110‑3 (MQ=11)
                                                                                              cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         >  2:351678/1‑108 (MQ=11)
                                                                                              cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCg                                                  >  1:304536/1‑149 (MQ=2)
                                                                                              cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAg                                                <  1:381132/151‑1 (MQ=2)
                                                                                               ggCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                  >  2:915459/1‑149 (MQ=31)
                                                                                               ggCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAg                                                <  1:290307/151‑1 (MQ=31)
                                                                                                     gggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        >  2:433298/1‑104 (MQ=11)
                                                                                                      ggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         >  1:178067/1‑100 (MQ=11)
                                                                                                       gAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACt                                        <  1:443573/151‑1 (MQ=31)
                                                                                                           gAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaat                                         >  2:437995/1‑144 (MQ=2)
                                                                                                           gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                   <  2:843174/151‑1 (MQ=2)
                                                                                                           gAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGc                                    <  1:433355/151‑1 (MQ=31)
                                                                                                             aataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        <  2:327888/93‑1 (MQ=11)
                                                                                                             aaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGt                                  >  1:733558/1‑151 (MQ=31)
                                                                                                              ataGAACCTGAAACCGTATGCGTACAAGCAGTGGAAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGCATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAgagatcggaag                                      >  2:79487/3‑135 (MQ=2)
                                                                                                               taGAACCTGAAACCGTATGCGTACAAGCAGTGGTAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         >  1:346374/2‑91 (MQ=2)
                                                                                                               aaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTa                                <  2:162931/151‑1 (MQ=31)
                                                                                                                aGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGttt                                >  2:422283/1‑148 (MQ=2)
                                                                                                                    ccTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTg                           <  1:536750/151‑1 (MQ=31)
                                                                                                                     cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGct                         >  2:294245/1‑150 (MQ=2)
                                                                                                                      tGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGctg                        >  1:360207/1‑149 (MQ=2)
                                                                                                                       gAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGc                          >  1:604545/1‑149 (MQ=31)
                                                                                                                       gAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTgg                       >  1:417148/1‑151 (MQ=2)
                                                                                                                       gAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTgg                       <  1:294211/151‑1 (MQ=2)
                                                                                                                         aaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAggg                      <  2:3181/151‑1 (MQ=31)
                                                                                                                          aCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGt                     <  2:50588/151‑1 (MQ=31)
                                                                                                                           ccGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGTGAGTCTTAATAGGGc                                    >  2:183415/1‑134 (MQ=2)
                                                                                                                            cgtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     <  1:433081/62‑1 (MQ=2)
                                                                                                                            cgtacgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    <  1:398350/127‑1 (MQ=2)
                                                                                                                             gtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                   >  1:243001/6‑149 (MQ=2)
                                                                                                                              tatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaat                                         <  1:209832/122‑2 (MQ=2)
                                                                                                                              tatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaat                                         >  2:209848/5‑125 (MQ=2)
                                                                                                                              tatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                   >  1:413846/5‑148 (MQ=2)
                                                                                                                               gtacgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                   >  2:668462/1‑148 (MQ=31)
                                                                                                                               gtacgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGGATAg                 <  2:343914/150‑1 (MQ=14)
                                                                                                                               atgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCTTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    <  1:126140/127‑1 (MQ=2)
                                                                                                                                 gcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGCCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGAccc             >  2:74060/2‑151 (MQ=2)
                                                                                                                                  cgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAACAGGGc                                    <  2:276219/127‑1 (MQ=2)
                                                                                                                                   gtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            >  1:406657/1‑150 (MQ=2)
                                                                                                                                   gtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            >  2:382346/1‑150 (MQ=2)
                                                                                                                                      cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGACGCTGGGTATAGACCCGAAAc        >  2:596312/1‑151 (MQ=2)
                                                                                                                                           aGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCgg    >  1:32192/1‑151 (MQ=31)
                                                                                                                                            gTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGt   <  1:464481/151‑1 (MQ=31)
                                                                                                                                              gggAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTg  >  2:198619/1‑150 (MQ=31)
                                                                                                                                                      |                                                                                                                                             
CCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTG  >  NC_000913/226151‑226442

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