Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 226,307 4 bp→TTAG noncoding (549‑552/2904 nt) rrlH → 23S ribosomal RNA of rrnH operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,3070GT81.5% 158.4 / 175.7 294noncoding (549/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base T (132/106);  total (162/132)
Fisher's exact test for biased strand distribution p-value = 8.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913226,3080CT81.6% 165.0 / 144.8 294noncoding (550/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base C (29/25);  new base T (133/107);  total (162/132)
Fisher's exact test for biased strand distribution p-value = 8.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913226,3090GA81.9% 167.0 / 143.4 293noncoding (551/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base A (133/107);  total (161/132)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913226,3100TG82.2% 166.7 / 159.8 293noncoding (552/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base G (134/107);  total (161/132)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 294 total aligned reads displayed.

ATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATCTAGCC  >  NC_000913/226162‑226450
                                                                                                                                                 ||||                                                                                                                                            
         ccATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTAc                                                                                                                                  >  2:404688/1‑151 (MQ=1)
              atccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCtttt                                                                                                                             >  2:158577/2‑151 (MQ=11)
                    ggCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGtat                                                                                                                         >  2:454513/1‑150 (MQ=31)
                     gCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAGCCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAAt                                                                                                                      <  1:365822/151‑1 (MQ=11)
                        aaaTACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATgg                                                                                                                    <  1:346015/150‑1 (MQ=11)
                              tactgactgaCCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             >  1:218288/3‑151 (MQ=11)
                              tactgactgaCCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             >  1:249728/3‑151 (MQ=11)
                               actgactgaCCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             <  2:57110/149‑1 (MQ=11)
                                 tgactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTTACTGCGTACCTTTTGTATAAt                                                                                                                      >  1:336031/1‑139 (MQ=11)
                                 tgactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTTACTGCGTACCTTTTGTATAAt                                                                                                                      <  2:335991/139‑1 (MQ=11)
                                  gactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGa                                                                                                            <  1:769123/149‑1 (MQ=31)
                                      gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTATCTTTTGTATAATGGGTCAgtcagtagt                                                                                                       <  2:85576/149‑9 (MQ=11)
                                      gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                                                                     <  2:24601/151‑1 (MQ=11)
                                      gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                                                                     >  2:639753/1‑151 (MQ=11)
                                      gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCCGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                                                                     <  2:118412/151‑1 (MQ=11)
                                        ccGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATgg                                                                                                                    <  2:38650/134‑1 (MQ=11)
                                        ccGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATgg                                                                                                                    >  1:38654/1‑134 (MQ=11)
                                        ccGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTc                                                                                                   <  1:404733/151‑1 (MQ=11)
                                        ccGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GATAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTc                                                                                                   >  1:107626/1‑151 (MQ=11)
                                         cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTcg                                                                                                                >  1:283851/1‑136 (MQ=11)
                                         cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             >  1:231598/1‑140 (MQ=11)
                                         cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             <  2:231568/140‑1 (MQ=11)
                                         cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTc                                                                                                   >  1:262495/1‑150 (MQ=11)
                                          gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCt                                                                                                  <  2:1989/151‑1 (MQ=31)
                                          gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             <  1:125051/139‑1 (MQ=11)
                                             aGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGtg                                                                                               <  1:291822/150‑2 (MQ=11)
                                              gTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGGGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                             >  1:63943/1‑151 (MQ=11)
                                                gAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCag                                                                                           <  1:376862/151‑2 (MQ=11)
                                                  aCCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAgc                                                                                         <  2:262461/151‑2 (MQ=11)
                                                    cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGt                                                                                        >  2:367410/1‑151 (MQ=31)
                                                      gTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             <  1:323296/127‑1 (MQ=11)
                                                      gTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             >  2:323257/1‑127 (MQ=11)
                                                        aCCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAAc                                                                                    >  1:113550/1‑151 (MQ=31)
                                                         ccGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAAcc                                                                                   >  1:737206/1‑151 (MQ=31)
                                                           gTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGt                                                                                                                           <  1:384158/108‑1 (MQ=11)
                                                           gTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGt                                                                                                                           >  2:384115/1‑108 (MQ=11)
                                                           gTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGAGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGAc                                                                                                           <  1:258076/124‑1 (MQ=11)
                                                               gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGc                                                                                                              <  1:200331/117‑1 (MQ=11)
                                                               gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAgg                                                                            >  2:292261/1‑151 (MQ=11)
                                                               gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTG‑AACCGTGTGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAggg                                                                           >  1:172674/1‑151 (MQ=11)
                                                                ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGt                                                                                                                  <  1:192014/112‑1 (MQ=11)
                                                                ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                                 <  1:181748/129‑3 (MQ=11)
                                                                ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                                 <  2:326214/129‑3 (MQ=11)
                                                                ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                                 >  2:181730/1‑127 (MQ=11)
                                                                ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                                >  1:331272/1‑130 (MQ=11)
                                                                ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                                <  2:331232/130‑1 (MQ=11)
                                                                  aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             <  2:310686/115‑1 (MQ=11)
                                                                  aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                                 >  1:303164/1‑125 (MQ=11)
                                                                  aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAgggt                                                                          >  1:160722/1‑149 (MQ=11)
                                                                    aGGCGAAAAGAACCCCGGAGAGTGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGGCTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTcag                                                                                                 <  1:149445/125‑3 (MQ=11)
                                                                     ggCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             >  2:137510/1‑112 (MQ=11)
                                                                      gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCg                                                                                                             >  2:389959/1‑111 (MQ=11)
                                                                      gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCg                                                                     >  2:369086/1‑151 (MQ=11)
                                                                       cGAAAGGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgt                                                                    >  1:156087/1‑150 (MQ=11)
                                                                       cGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCg                                                                     <  2:494139/151‑1 (MQ=31)
                                                                       cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCTTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgt                                                                    >  1:390936/1‑150 (MQ=11)
                                                                       cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                             >  1:471735/1‑126 (MQ=11)
                                                                       cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgt                                                                    >  2:286449/1‑150 (MQ=11)
                                                                       cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCg                                                                     >  2:68788/1‑150 (MQ=11)
                                                                          aaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAgc                                                                 >  2:190774/1‑150 (MQ=11)
                                                                          aaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACAGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                                >  2:149431/1‑120 (MQ=11)
                                                                           aaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAg                                                                                          <  2:146729/125‑1 (MQ=11)
                                                                              aaCCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACCTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGa                                                               >  2:670955/1‑150 (MQ=14)
                                                                              aaCCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTca                                                                                                  >  2:237538/1‑113 (MQ=11)
                                                                              aaCCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTca                                                                                                  <  1:237568/114‑2 (MQ=11)
                                                                               aCCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGaaa                                                             <  1:14136/151‑1 (MQ=31)
                                                                                 cccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAAcc                                                           <  1:271954/151‑1 (MQ=31)
                                                                                 cccGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGTCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGaaagc                                                           >  1:536041/1‑148 (MQ=2)
                                                                                  ccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                          >  1:84132/1‑151 (MQ=31)
                                                                                  ccGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                                                                <  2:212484/112‑1 (MQ=11)
                                                                                   cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGcg                                                                                            >  2:195872/1‑114 (MQ=11)
                                                                                        aGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTGCCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCtt                                                    >  1:16818/1‑150 (MQ=2)
                                                                                            gAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACt                                                >  2:728845/1‑151 (MQ=31)
                                                                                              gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTGGCGAAAGCGAGTCTTaataggg                                             >  2:4140/1‑146 (MQ=2)
                                                                                                gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTGCCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                            <  2:16817/150‑1 (MQ=2)
                                                                                                gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTGGCGAAAGCGAGTCTTAATAGGGc                                            <  1:4140/150‑1 (MQ=2)
                                                                                                  aataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTtata                                                                                                      <  1:150823/87‑1 (MQ=11)
                                                                                                  aaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGt                                          <  2:354954/151‑1 (MQ=31)
                                                                                                   ataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGttt                                        >  2:345726/3‑150 (MQ=2)
                                                                                                     aGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGttt                                        >  2:173935/1‑148 (MQ=2)
                                                                                                      gAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaat                                                 <  2:252586/139‑2 (MQ=2)
                                                                                                          cTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTg                                   >  2:596397/1‑150 (MQ=31)
                                                                                                          cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGct                                 >  1:195971/1‑150 (MQ=2)
                                                                                                          cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGct                                 >  1:601118/1‑150 (MQ=2)
                                                                                                            gAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                                                                     <  1:31098/81‑1 (MQ=11)
                                                                                                            gAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTgg                               >  2:282889/1‑151 (MQ=2)
                                                                                                             aaaCCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGctg                                >  2:15521/1‑147 (MQ=2)
                                                                                                                ccGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                           >  1:207174/1‑151 (MQ=2)
                                                                                                                 cgtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                           >  2:220167/7‑150 (MQ=2)
                                                                                                                 cgtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                           >  1:49923/7‑150 (MQ=2)
                                                                                                                   tatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                           >  2:132883/5‑148 (MQ=2)
                                                                                                                    atgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                             <  1:343034/126‑6 (MQ=2)
                                                                                                                    atgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                            <  1:9782/127‑1 (MQ=2)
                                                                                                                        gtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaatagg                                              <  2:144047/124‑5 (MQ=2)
                                                                                                                        gtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGtt                                         <  1:144270/129‑1 (MQ=2)
                                                                                                                        gtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCg                    <  2:439420/151‑1 (MQ=31)
                                                                                                                              gCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAAccc              <  2:468243/151‑1 (MQ=31)
                                                                                                                                   gggAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTg          >  2:527005/1‑150 (MQ=31)
                                                                                                                                    ggAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGAt        <  1:692708/151‑1 (MQ=31)
                                                                                                                                          tGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATCTAGcc  >  2:54024/1‑151 (MQ=11)
                                                                                                                                                 ||||                                                                                                                                            
ATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATCTAGCC  >  NC_000913/226162‑226450

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