Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,3100TG68.7% 39.2 / 231.9 214noncoding (552/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base T (38/29);  new base G (83/64);  total (121/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

TATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTG  >  NC_000913/226161‑226442
                                                                                                                                                     |                                                                                                                                    
tatGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCgtg                                                                                                                                     >  1:152308/1‑151 (MQ=31)
         aCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGt                                                                                                                            >  2:699077/1‑151 (MQ=31)
               ctccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAATCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCtttt                                                                                                                     >  2:210509/1‑151 (MQ=11)
                tccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTg                                                                                                                    >  1:145477/1‑151 (MQ=11)
                 ccAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTg                                                                                                                    <  2:111520/151‑1 (MQ=31)
                  cAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGt                                                                                                                   <  2:166615/151‑1 (MQ=31)
                   aaGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGt                                                                                                                   >  2:350913/1‑150 (MQ=31)
                     ggCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGtat                                                                                                                 >  2:170304/1‑150 (MQ=31)
                       cTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAAt                                                                                                              >  2:759384/1‑151 (MQ=31)
                         aaaTACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGt                                                                                                          >  1:168936/1‑151 (MQ=11)
                              ctactgactgaCCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGc                                                                                                      >  2:64417/4‑151 (MQ=11)
                                   gactgacCGATAGTGAACTAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGtata                                                                                                                >  2:27888/1‑136 (MQ=11)
                                   gactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     >  2:46592/1‑147 (MQ=11)
                                    actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACtt                                                                                                 >  1:79909/1‑151 (MQ=31)
                                    actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACtt                                                                                                 >  1:274192/1‑151 (MQ=31)
                                     ctgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACGCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGTCAGCGACTtat                                                                                               >  1:442505/1‑150 (MQ=11)
                                         ccGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATAtt                                                                                            >  2:443621/1‑151 (MQ=31)
                                         ccGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATAtt                                                                                            >  1:414532/1‑151 (MQ=31)
                                          cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTc                                                                                           >  1:632530/1‑151 (MQ=31)
                                          cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGtat                                                                                                                 <  1:198270/128‑1 (MQ=11)
                                          cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                                                         >  2:307809/1‑136 (MQ=11)
                                            aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAATCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        <  1:210493/151‑1 (MQ=11)
                                                    ggaGTACAGTGAGGGAAAGGAGAAAAGAAACCAGGAGAGGGGAGTGAAATAGAACCTGAAACAGTATGAGTACAAGCAGTGGGAGCCTAGTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                     >  1:102381/3‑146 (MQ=2)
                                                    ccAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        <  2:158727/143‑1 (MQ=11)
                                                     cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGt                                                                                >  2:544075/1‑151 (MQ=31)
                                                     cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        >  1:232556/1‑142 (MQ=11)
                                                      aGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGtt                                                                               <  2:195333/151‑1 (MQ=31)
                                                       gTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                     >  1:472846/1‑143 (MQ=11)
                                                          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)
                                                             tGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGa                                                                         >  1:445041/1‑150 (MQ=31)
                                                              gAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGaa                                                                        <  2:572705/150‑1 (MQ=31)
                                                              gAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAAt                                                                       >  1:50583/1‑151 (MQ=31)
                                                               gaggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                     <  2:375677/133‑1 (MQ=11)
                                                                gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTca                                                                                          <  1:996/129‑2 (MQ=11)
                                                                 gggaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGAc                                                                                                   <  1:655929/116‑1 (MQ=11)
                                                                 ggAAAGGCTAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                                                         <  2:489120/113‑1 (MQ=11)
                                                                 ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGtg                                                                                       >  1:75039/1‑130 (MQ=11)
                                                                   aaaGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAgggg                                                                  >  1:647575/1‑151 (MQ=31)
                                                                   aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAg                                                                                                       >  2:244893/1‑113 (MQ=11)
                                                                   aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGAc                                                                                                   <  1:256556/117‑1 (MQ=11)
                                                                       gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATgggc                                                                                                          >  2:493902/1‑105 (MQ=11)
                                                                       gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCa                                                                                    >  2:318831/1‑128 (MQ=11)
                                                                        cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATg                                                                                                             <  1:36884/102‑1 (MQ=11)
                                                                         gAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCg                                                             >  1:412868/1‑150 (MQ=31)
                                                                          aaaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTTGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     >  1:453750/1‑108 (MQ=11)
                                                                          aaaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTTGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     <  2:453801/108‑1 (MQ=11)
                                                                          aaaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         <  1:432063/120‑3 (MQ=11)
                                                                               aaCCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGaa                                                      <  1:389138/151‑1 (MQ=31)
                                                                               aaCCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGaaa                                                     <  1:217945/151‑1 (MQ=2)
                                                                                 ccccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAAc                                                    >  1:380350/1‑151 (MQ=31)
                                                                                   ccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGaaa                                                     <  1:894013/148‑1 (MQ=31)
                                                                                    cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAATCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     >  2:209492/1‑98 (MQ=2)
                                                                                    cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAATCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     <  1:209476/98‑1 (MQ=2)
                                                                                    cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         <  1:351640/110‑3 (MQ=11)
                                                                                    cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCg                                                  >  1:304536/1‑149 (MQ=2)
                                                                                     ggCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                  >  2:915459/1‑149 (MQ=31)
                                                                                     ggCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                  >  1:375737/1‑149 (MQ=31)
                                                                                     ggCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAg                                                <  1:290307/151‑1 (MQ=31)
                                                                                           gggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        <  1:433250/104‑1 (MQ=11)
                                                                                           gggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        >  2:433298/1‑104 (MQ=11)
                                                                                           gggAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTa                                           >  1:971769/1‑150 (MQ=31)
                                                                                            ggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         <  2:178081/102‑3 (MQ=11)
                                                                                            ggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         >  1:178067/1‑100 (MQ=11)
                                                                                            ggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaat                                         >  1:843082/1‑149 (MQ=2)
                                                                                               gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     <  2:207179/87‑1 (MQ=11)
                                                                                               gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        >  1:490386/1‑100 (MQ=11)
                                                                                               gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     <  2:356917/151‑6 (MQ=2)
                                                                                               gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTAAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     <  2:33314/151‑6 (MQ=2)
                                                                                                tGAAAAAGAACCTGAAACCGTATACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTggg                                     >  2:229587/1‑151 (MQ=14)
                                                                                                 gAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaat                                         <  1:437947/145‑2 (MQ=2)
                                                                                                 gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    <  2:304572/150‑1 (MQ=2)
                                                                                                 gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                   <  2:843174/151‑1 (MQ=2)
                                                                                                   aataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        >  1:327852/4‑96 (MQ=11)
                                                                                                    ataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTGAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAg                                                >  1:575716/3‑135 (MQ=2)
                                                                                                    ataGAACCTGAAACCGTATGCGTACAAGCAGTGGAAGCCTACTTGTTAGGTGACTGCGTACCTTTTGCATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAgagatcggaag                                      >  2:79487/3‑135 (MQ=2)
                                                                                                     taGAACCTGAAACCGTATGCGTACAAGCAGTGGTAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         >  1:346374/2‑91 (MQ=2)
                                                                                                      aGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGttt                                >  2:422283/1‑148 (MQ=2)
                                                                                                          ccTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTg                           <  1:536750/151‑1 (MQ=31)
                                                                                                           cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGct                         >  2:294245/1‑150 (MQ=2)
                                                                                                            tGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGctg                        >  1:360207/1‑149 (MQ=2)
                                                                                                             gAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGc                          >  1:604545/1‑149 (MQ=31)
                                                                                                             gAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTgg                       <  1:294211/151‑1 (MQ=2)
                                                                                                                aCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGt                     <  2:50588/151‑1 (MQ=31)
                                                                                                                  cgtacgcgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    >  2:398396/7‑133 (MQ=2)
                                                                                                                  cGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                   <  2:749835/151‑1 (MQ=31)
                                                                                                                   gtatgcgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                   >  1:243001/6‑149 (MQ=2)
                                                                                                                     gtacgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                   >  2:668462/1‑148 (MQ=31)
                                                                                                                     gtacgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGGATAg                 <  2:343914/150‑1 (MQ=14)
                                                                                                                     atgcgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCTTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    <  1:126140/127‑1 (MQ=2)
                                                                                                                     atgcgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCTTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    >  2:126150/4‑130 (MQ=2)
                                                                                                                        cgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            >  2:440056/1‑151 (MQ=2)
                                                                                                                        cgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            >  2:314180/1‑151 (MQ=2)
                                                                                                                        cgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAACAGGGc                                    <  2:276219/127‑1 (MQ=2)
                                                                                                                         gtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            >  1:406657/1‑150 (MQ=2)
                                                                                                                            cAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGACGCTGGGTATAGACCCGAAAc        >  2:596312/1‑151 (MQ=2)
                                                                                                                                 aGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCgg    >  1:32192/1‑151 (MQ=31)
                                                                                                                                  gTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGt   <  1:464481/151‑1 (MQ=31)
                                                                                                                                                     |                                                                                                                                    
TATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTG  >  NC_000913/226161‑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: