Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913226,3570AG70.5% 21.7 / 107.8 129noncoding (599/2904 nt)rrlH23S ribosomal RNA of rrnH operon
Reads supporting (aligned to +/- strand):  ref base A (22/16);  new base G (58/33);  total (80/49)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 129 total aligned reads displayed.

GTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATCTAGCC  >  NC_000913/226208‑226450
                                                                                                                                                     |                                                                                             
gTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGGGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                             >  1:63943/1‑151 (MQ=11)
gTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAGGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAg                                                                                              <  2:642151/151‑1 (MQ=14)
  gAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCag                                                                                           <  1:376862/151‑2 (MQ=11)
  gAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCa                                                                                            <  1:4085/150‑1 (MQ=11)
    aCCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAgc                                                                                         <  2:262461/151‑2 (MQ=11)
      cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGt                                                                                        >  2:367410/1‑151 (MQ=31)
      cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCtt                                                                                       >  1:150063/1‑151 (MQ=11)
          aCCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAAc                                                                                    >  1:113550/1‑151 (MQ=31)
          aCCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAAc                                                                                    <  2:150047/150‑1 (MQ=11)
           ccGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAAcc                                                                                   >  1:737206/1‑151 (MQ=31)
            cGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCg                                                                                  <  2:767279/151‑1 (MQ=31)
               gAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCgtac                                                                               <  1:243642/150‑4 (MQ=11)
                 gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAgg                                                                            >  2:292261/1‑151 (MQ=11)
                 gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTG‑AACCGTGTGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAggg                                                                           >  1:172674/1‑151 (MQ=11)
                  ggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAggg                                                                           <  1:292298/151‑1 (MQ=11)
                    aaaGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAgggg                                                                          >  1:2642/1‑151 (MQ=31)
                    aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAgggt                                                                          >  1:160722/1‑149 (MQ=11)
                       ggCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGc                                                                       >  1:714751/1‑151 (MQ=31)
                       ggCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTAt                                                                               <  2:352273/142‑1 (MQ=11)
                       ggCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTAt                                                                               >  1:352315/1‑142 (MQ=11)
                       ggCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAggc                                                                      >  1:167546/1‑149 (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)
                         cGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCg                                                                     >  2:692472/1‑151 (MQ=31)
                         cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCTTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgt                                                                    >  1:390936/1‑150 (MQ=11)
                         cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                             <  2:471682/126‑1 (MQ=11)
                         cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGc                                                                                             >  1:471735/1‑126 (MQ=11)
                         cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgt                                                                    >  1:507041/1‑150 (MQ=11)
                         cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCg                                                                     >  2:68788/1‑150 (MQ=11)
                          gAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGTGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgt                                                                    >  1:79361/1‑149 (MQ=11)
                           aaaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAg                                                                  <  2:148238/151‑1 (MQ=11)
                            aaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAgc                                                                 >  2:190774/1‑150 (MQ=11)
                            aaaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAgc                                                                 >  1:340748/1‑150 (MQ=11)
                             aaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAg                                                                                          >  1:146745/1‑125 (MQ=11)
                             aaGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAg                                                                                          <  2:146729/125‑1 (MQ=11)
                                aaCCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGa                                                               >  2:324002/1‑150 (MQ=31)
                                 aCCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGaaa                                                             <  1:14136/151‑1 (MQ=31)
                                   cccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAAcc                                                           <  1:271954/151‑1 (MQ=31)
                                   cccGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGTCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGaaagc                                                           >  1:536041/1‑148 (MQ=2)
                                    ccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTCTATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                          >  2:621094/1‑151 (MQ=14)
                                    ccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCg                                                          >  1:84132/1‑151 (MQ=31)
                                     cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGTCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGaaagcaag                                                        <  2:535979/151‑6 (MQ=2)
                                     cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGcg                                                                                            >  2:195872/1‑114 (MQ=11)
                                       gCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGt                                                       <  2:433274/151‑1 (MQ=31)
                                          aGGGGAGTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTGCCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCtt                                                    >  1:16818/1‑150 (MQ=2)
                                           ggggAGTGAAATAGAACCTG‑AACCGTGTGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAAc                                                 <  2:172656/151‑1 (MQ=2)
                                              gAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACt                                                >  2:728845/1‑151 (MQ=31)
                                                gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCa                                                                                            >  1:203001/1‑104 (MQ=11)
                                                gTGAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTGGCGAAAGCGAGTCTTaataggg                                             >  2:4140/1‑146 (MQ=2)
                                                  gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTGCCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                            <  2:16817/150‑1 (MQ=2)
                                                  gAAATAGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTGGCGAAAGCGAGTCTTAATAGGGc                                            <  1:4140/150‑1 (MQ=2)
                                                   aaataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAAcc                                                                                   >  1:189238/1‑110 (MQ=11)
                                                   aaataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAAcc                                                                                   <  2:189219/110‑1 (MQ=11)
                                                   aaataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGCCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGt                                          >  2:371887/1‑151 (MQ=2)
                                                    actaGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGtt                                         >  1:244621/4‑151 (MQ=2)
                                                    aataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGtt                                         >  1:19675/4‑151 (MQ=2)
                                                    aataGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGtt                                         >  1:630479/4‑151 (MQ=2)
                                                    aaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGt                                          <  2:354954/151‑1 (MQ=31)
                                                    aaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGt                                          <  1:653809/151‑1 (MQ=31)
                                                       aGAACCTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGttt                                        >  2:173935/1‑148 (MQ=2)
                                                         aaCCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGt                                     >  2:328758/1‑151 (MQ=14)
                                                            cTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTg                                   >  2:596397/1‑150 (MQ=31)
                                                            cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGct                                 >  1:601118/1‑150 (MQ=2)
                                                            cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGct                                 >  2:251763/1‑150 (MQ=2)
                                                            cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGct                                 >  1:195971/1‑150 (MQ=2)
                                                            cTGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGc                                  >  2:520307/1‑150 (MQ=2)
                                                             tGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataga                                              >  2:250167/1‑133 (MQ=2)
                                                             tGAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataga                                              <  1:250201/137‑5 (MQ=2)
                                                              gAAACCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTgg                               >  2:282889/1‑151 (MQ=2)
                                                               aaaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAgg                               >  2:534496/1‑151 (MQ=31)
                                                               aaaCCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGctg                                >  2:15521/1‑147 (MQ=2)
                                                               aaaCCGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTggg                              >  2:197708/1‑151 (MQ=2)
                                                                  ccGTATGCGTACAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                           >  1:207174/1‑151 (MQ=2)
                                                                   cgtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                            <  1:284079/127‑1 (MQ=2)
                                                                   cgtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGc                                            >  2:284042/7‑133 (MQ=2)
                                                                   cgtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                           >  2:220167/7‑150 (MQ=2)
                                                                   cgtatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCAACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCCGGGtat                           >  1:110912/7‑150 (MQ=2)
                                                                   cGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                           <  2:43797/151‑1 (MQ=31)
                                                                    gtgtACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtata                          <  2:660402/151‑1 (MQ=31)
                                                                    gtgtACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGtat                           >  2:680368/1‑150 (MQ=31)
                                                                     tgtACGTGCAAGCAGTGGGGGCATGCTTAGGCGTGTGACTGCGTCCCTTTTGTATAATGGTTCAGCGACTTATATTCTGTATCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGTtata                          <  2:368872/150‑1 (MQ=11)
                                                                     tatgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGtat                           >  2:132883/5‑148 (MQ=2)
                                                                      atgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                             <  1:343034/126‑6 (MQ=2)
                                                                      atgcgtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTaataggg                                             >  2:342994/4‑124 (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)
                                                                          gtacAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCg                    >  1:728775/1‑150 (MQ=2)
                                                                            acAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGaa                  <  1:778766/151‑1 (MQ=31)
                                                                             cAAGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGTTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAc                >  1:172924/1‑151 (MQ=2)
                                                                              aaGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTTTATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAACAGGGc                                            >  1:384885/1‑122 (MQ=2)
                                                                              aaGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgta                                                                   <  2:142879/99‑3 (MQ=2)
                                                                              aaGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCgta                                                                   >  1:142893/1‑97 (MQ=2)
                                                                              aaGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTTTAGTCGCTGGGTATAGACCCGAAAcc               >  2:739791/1‑151 (MQ=2)
                                                                              aaGCAGTGGGAGCCTACTT‑GTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCAAGCTTAACCGTATAGGGTAGGCGTAGCGAAAGCGAGTCTTAACAGGGc                                            <  2:384842/122‑1 (MQ=2)
                                                                                gCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAAccc              <  2:468243/151‑1 (MQ=31)
                                                                                 cAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCg             >  1:653536/1‑151 (MQ=31)
                                                                                  aGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCgg            >  1:382834/1‑151 (MQ=31)
                                                                                            tGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATCTAGcc  >  2:54024/1‑151 (MQ=11)
                                                                                                                                                     |                                                                                             
GTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCGAATAGGGGAGCCGAAGGGAAACCGAGTCTTAACTGGGCGTTAAGTTGCAGGGTATAGACCCGAAACCCGGTGATCTAGCC  >  NC_000913/226208‑226450

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