Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,3070GT68.8% 41.4 / 250.1 215noncoding (549/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base G (39/28);  new base T (83/65);  total (122/93)
Fisher's exact test for biased strand distribution p-value = 8.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

GAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTG  >  NC_000913/226158‑226442
                                                                                                                                                     |                                                                                                                                       
   tatGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCgtg                                                                                                                                     >  1:152308/1‑151 (MQ=31)
           gACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCg                                                                                                                             >  1:953425/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)
                     cAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGt                                                                                                                   <  2:166615/151‑1 (MQ=31)
                          cTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAAt                                                                                                              >  2:759384/1‑151 (MQ=31)
                            aaaTACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGt                                                                                                          >  1:168936/1‑151 (MQ=11)
                             aaTACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATggg                                                                                                           >  1:98831/1‑151 (MQ=31)
                                 ctactgactgaCCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGc                                                                                                      >  2:64417/4‑151 (MQ=11)
                                 cTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                                                                                         >  1:441227/1‑149 (MQ=31)
                                      gactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     >  2:46592/1‑147 (MQ=11)
                                       actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACtt                                                                                                 >  1:79909/1‑151 (MQ=31)
                                            ccGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATAtt                                                                                            >  1:414532/1‑151 (MQ=31)
                                             cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACGCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGTCAGCGACTTATATTcag                                                                                         <  2:442553/151‑3 (MQ=11)
                                             cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGtat                                                                                                                 >  2:198284/1‑128 (MQ=11)
                                             cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGtat                                                                                                                 <  1:198270/128‑1 (MQ=11)
                                             cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                                                         >  2:307809/1‑136 (MQ=11)
                                              ggtAGTGGACCAGTACCCTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         <  2:347696/149‑3 (MQ=11)
                                               aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAATCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        <  1:210493/151‑1 (MQ=11)
                                                  gTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAATCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGtgg                                                                                      <  1:517567/150‑3 (MQ=11)
                                                  gTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATAttaagtggc                                                                                     >  1:33312/1‑144 (MQ=11)
                                                     aaCCAGTACAGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGcg                                                                                    <  2:177088/149‑2 (MQ=11)
                                                      aCCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAg                                                                                  <  2:119309/151‑1 (MQ=31)
                                                       ggaGTACAGTGAGGGAAAGGAGAAAAGAAACCAGGAGAGGGGAGTGAAATAGAACCTGAAACAGTATGAGTACAAGCAGTGGGAGCCTAGTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                     >  1:102381/3‑146 (MQ=2)
                                                       ccAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        <  2:158727/143‑1 (MQ=11)
                                                       ccAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCa                                                                                    <  2:102389/147‑1 (MQ=11)
                                                        cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGt                                                                                >  2:544075/1‑151 (MQ=31)
                                                        cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        <  2:232575/142‑1 (MQ=11)
                                                        cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        >  1:232556/1‑142 (MQ=11)
                                                          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                                                                         >  1:445041/1‑150 (MQ=31)
                                                                 gAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGaa                                                                        <  2:572705/150‑1 (MQ=31)
                                                                  gaggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                     <  2:375677/133‑1 (MQ=11)
                                                                    ggAAAGGCTAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                                                         <  2:489120/113‑1 (MQ=11)
                                                                    ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         <  2:294616/129‑3 (MQ=11)
                                                                    ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGtg                                                                                       >  1:75039/1‑130 (MQ=11)
                                                                      aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGAc                                                                                                   <  1:256556/117‑1 (MQ=11)
                                                                      aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGAc                                                                                                   >  2:256579/1‑117 (MQ=11)
                                                                         ggCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAggc                                                              >  2:52429/1‑149 (MQ=11)
                                                                          gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATgggc                                                                                                          <  1:493849/106‑2 (MQ=11)
                                                                          gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATgggc                                                                                                          >  2:493902/1‑105 (MQ=11)
                                                                          gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCa                                                                                    >  2:318831/1‑128 (MQ=11)
                                                                           cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATg                                                                                                             <  1:36884/102‑1 (MQ=11)
                                                                           cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgt                                                            >  2:285529/1‑150 (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)
                                                                                  aaCCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGaaa                                                     <  1:217945/151‑1 (MQ=2)
                                                                                       cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAATCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     <  1:209476/98‑1 (MQ=2)
                                                                                       cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAATCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     >  2:209492/1‑98 (MQ=2)
                                                                                       cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         <  1:351640/110‑3 (MQ=11)
                                                                                       cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCg                                                  >  1:304536/1‑149 (MQ=2)
                                                                                       cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAg                                                <  1:381132/151‑1 (MQ=2)
                                                                                        ggCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                  >  2:915459/1‑149 (MQ=31)
                                                                                        ggCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                  >  1:375737/1‑149 (MQ=31)
                                                                                        ggCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAg                                                <  1:290307/151‑1 (MQ=31)
                                                                                           gAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCt                                             >  1:443250/1‑151 (MQ=31)
                                                                                             ggggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTa                                           >  2:291358/1‑150 (MQ=2)
                                                                                              gggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        <  1:433250/104‑1 (MQ=11)
                                                                                              gggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                        >  2:433298/1‑104 (MQ=11)
                                                                                               ggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                         >  1:178067/1‑100 (MQ=11)
                                                                                               ggAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaat                                         >  1:843082/1‑149 (MQ=2)
                                                                                                  gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                     <  2:207179/87‑1 (MQ=11)
                                                                                                  gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     <  2:356917/151‑6 (MQ=2)
                                                                                                  gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTAAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                     <  2:33314/151‑6 (MQ=2)
                                                                                                    gAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaat                                         >  2:437995/1‑144 (MQ=2)
                                                                                                    gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCg                                   <  2:843174/151‑1 (MQ=2)
                                                                                                    gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAACAGGGc                                    <  1:64412/150‑1 (MQ=2)
                                                                                                       ataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTGAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAg                                                >  1:575716/3‑135 (MQ=2)
                                                                                                       ataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTACAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGttt                                >  2:278231/3‑150 (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)
                                                                                                               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)
                                                                                                                    ccGTATGCGTACAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGTGAGTCTTAATAGGGc                                    >  2:183415/1‑134 (MQ=2)
                                                                                                                     cgtacgcgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    >  2:398396/7‑133 (MQ=2)
                                                                                                                     cgtacgcgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                    <  1:398350/127‑1 (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:314180/1‑151 (MQ=2)
                                                                                                                           cgtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAACAGGGc                                    >  1:276193/1‑127 (MQ=2)
                                                                                                                            gtacAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg            >  1:406657/1‑150 (MQ=2)
                                                                                                                             tacACGCAGTGGGAGCATGCGTAGTCGTGTGACTGCGTACCTTTTGTGTAATGGGTCAGTGACTTATATTCTGTAGCAAGGTTAATCGGATGGGGGATGGGATGGGAATCCGAGTCTTAACTGGGCGTTACGTTGCAGGGtatattc               <  1:694563/147‑4 (MQ=11)
                                                                                                                               cAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGACGCTGGGTATAGACCCGAAAc        >  2:596312/1‑151 (MQ=2)
                                                                                                                               cAAGCAGTGGGAGCCTACTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTAAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAc        >  1:649415/1‑151 (MQ=2)
                                                                                                                                    aGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCgg    >  1:32192/1‑151 (MQ=31)
                                                                                                                                     gTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGt   <  1:464481/151‑1 (MQ=31)
                                                                                                                                                     |                                                                                                                                       
GAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTG  >  NC_000913/226158‑226442

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