Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2940TC41.0% 135.3 / 25.1 95intergenic (+100/+129)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (16/40);  new base C (39/0);  total (55/40)
Fisher's exact test for biased strand distribution p-value = 6.08e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCT  >  NC_000913/3423145‑3423444
                                                                                                                                                      |                                                                                                                                                          
acccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAAACCATTAATTGCATGCTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCCCCCCCCTc                                                                                                                                                          >  2:1169389/3‑151 (MQ=255)
     aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGc                                                                                                                                                     <  2:1333715/151‑1 (MQ=255)
     aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCcgg                                                                                                                                                      >  1:725151/1‑150 (MQ=255)
      aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCGGGcc                                                                                                                                                    >  1:266711/1‑151 (MQ=255)
       aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCCCCCGGcc                                                                                                                                                    >  1:849782/1‑150 (MQ=255)
          agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGCATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCGCCACCCTCCggg                                                                                                                                                      >  2:709585/1‑145 (MQ=255)
          agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                <  2:1020349/150‑1 (MQ=255)
           gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGCCCCGGCCGCCCCCCCCCCGGcccgtg                                                                                                                                                >  2:904631/1‑146 (MQ=255)
              cGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCGCCCGCCACCCCCCcgg                                                                                                                                                      >  2:742935/1‑141 (MQ=255)
               gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGGTCGCCCCCCAACCCCCTGGCCGt                                                                                                                                                 >  2:1299687/1‑144 (MQ=37)
                  cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                      <  2:78504/151‑1 (MQ=255)
                    aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                    <  1:104597/151‑1 (MQ=255)
                     gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                                   <  2:250902/151‑1 (MQ=255)
                        cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTCGCCACAtt                                                                                                                                  >  1:1035337/1‑149 (MQ=25)
                        cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCCCCCCCCCGGCCGCTGGtt                                                                                                                                            >  2:797947/1‑140 (MQ=37)
                         aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                              <  1:1047441/151‑1 (MQ=255)
                         aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTATGGCAGTTTATGGCGGGCGGTCCTCCCTCCACCCCCCCGTCCGGGTCTTCGcccagctta                                                                                                                                 >  2:1288614/1‑142 (MQ=255)
                             ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGTTCAAATcc                                                                                                                          >  1:1346088/1‑151 (MQ=25)
                              tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCCCCCCCCCGGCCGGTGCTTCGCAACGTTCAAATccc                                                                                                                         >  1:416553/1‑150 (MQ=255)
                                  gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCCCCCCCCCCCCGGCCCTTTCTTCGCAACATTCAAAACCGCTcc                                                                                                                     >  1:1135709/1‑151 (MQ=255)
                                    aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCGCTcccc                                                                                                                   >  1:1388850/1‑150 (MQ=255)
                                     ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAAAcccccccc                                                                                                                    >  1:766768/1‑149 (MQ=21)
                                           ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                             <  1:1466829/150‑1 (MQ=255)
                                            cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCccgcggggtt                                                                                                           >  1:1272993/1‑147 (MQ=18)
                                            cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGTTTAAATCCCCTCCccgcggggtt                                                                                                           >  1:1311088/1‑147 (MQ=255)
                                             aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCCGCCCGCCCCCCCCCGCGCCGTGGGTTCCCAACGCTTCCATCCGCCCCCCGCCGGttt                                                                                                           >  1:397191/1‑150 (MQ=255)
                                             aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                          <  1:323719/151‑1 (MQ=255)
                                               tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCCGTCCGCCACCCCCCCGGCCGTTTCTTTGCTACGTTCAAAACCCGCCCCCGCCGGttt                                                                                                            >  2:539827/1‑148 (MQ=18)
                                                  ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCCCCCGGCCGTTGCTTCGCAACGGTCAAAACCGCCCCCCGCGGAttttttc                                                                                                       >  1:1455465/1‑145 (MQ=18)
                                                   tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCCGCCCCCCACCCCCCCGGCCGGTGCTGCGCAACGGTCACAACCCCTCCCCCCGCATTTGTTCttc                                                                                                    >  2:148815/1‑149 (MQ=18)
                                                    gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCACCCCCCCGGCCGTTTCTTCCCAAAGCTCAAAATCGCCCCCCGCGGGTTTGTTCTAc                                                                                                    >  2:1246775/1‑150 (MQ=255)
                                                        gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCgg                                                                                                <  1:175471/150‑1 (MQ=255)
                                                         ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGTCCCCCGCCCTCCCGGGCGTTGCTTTGAAACCTTCAAAACACCTCCCCGGGGGATTGTCCTACTc                                                                                                  >  2:863213/1‑147 (MQ=17)
                                                          tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGGTAGTTCGCAACCATCCAATTCCCACCCCGCCGCTTTGTCCCACTCcgg                                                                                               >  2:1512243/1‑149 (MQ=17)
                                                               cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtgt                                                                                        <  1:477220/151‑1 (MQ=255)
                                                               cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                         <  2:1324881/150‑1 (MQ=255)
                                                                  aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGGCTGCCCGCCCCCCCCCCGGCCGTTGCTTTGCAACCTTTAAATCCCCTCCCGGCGGAGTTGTTCTACTCcgggg                                                                                             >  1:732653/1‑141 (MQ=17)
                                                                       aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTCCCCCCCCCCCCCCCGGCCGGTTTGTTGCAACCTTCACATTCGCCCCCCGGGGGTTTGTCCCAccctgg                                                                                               >  2:845058/1‑131 (MQ=255)
                                                                       aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGAGGCCTGCCCGCCACCCCCCCGGCCGGTGCTTCCCCAAGTTTACAACTGCCCCCCGCGGGTTTTTTCTACTCTCGAGGGGGTTCTCCGac                                                                                >  2:1194788/1‑151 (MQ=255)
                                                                        aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGGCCTCCCCCCCCCCTCCCGGCCGGGTCTTTGCAACCATCAAATACGCCCCCCGCGGGGTTGTTCTTCTCGGGAGGGGGttttcc                                                                                   >  2:1106600/1‑143 (MQ=255)
                                                                        aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCCGGCCTGCCCCCCCCCCTCCCGGGCCGCTCTTCGCAAAATTTAAATACCCCCCCCGCGGGGTTGTCCTACTCCGGGGAgggtgctccc                                                                                  >  2:1259320/1‑142 (MQ=255)
                                                                          aatgCCATCCTTTCGGCTGTCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCGCCACCCTCCGGGTCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGTATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                             <  2:237985/148‑1 (MQ=17)
                                                                          aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                             <  1:1173829/151‑1 (MQ=255)
                                                                            ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                           <  2:109735/151‑1 (MQ=255)
                                                                            ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                           <  2:884933/151‑1 (MQ=255)
                                                                            ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                           <  2:1424295/151‑1 (MQ=255)
                                                                              ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCTCCCGGCCGTTGCTTCCCAACGTTCAAATCCCCTCCCCGCGGGTTTGTCCTACTCCGGGGAGTGGTCACCGccaaacaa                                                                          >  1:1573775/1‑150 (MQ=17)
                                                                               cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                        <  1:323475/151‑1 (MQ=255)
                                                                               cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCGGGCCGTTGCTTCGCCAAGATCAAATCCCCCCCCCGCGGGTTTGTCCTACTCCGGGGAGTGTTTACCCacaaacaa                                                                          >  1:954444/1‑149 (MQ=17)
                                                                               cATCCTTTCGGATGGCCTTTCGCTCGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                        <  1:1283453/151‑1 (MQ=255)
                                                                                  ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGGCCACCAACCCCCCAGCCGTTGCCTACCAAACGTCCAATCCCCTCCCCGCGGGTTTTTTCCACTCCGGGCAGAGGTCACCCacaaagaaaaaa                                                                      >  1:963655/1‑146 (MQ=9)
                                                                                   cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGGTTTGCCACCTTCAAATACCCCCCCCGCGGGTTTTTTCTTCTCCGGGGAGGGTTTACCCacaaaaaa                                                                          >  1:468817/1‑145 (MQ=255)
                                                                                    tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTGCTTCGCAACGTTCAAAACCCCCCCCCGCGGGTTTTTTCTACTCGGGGGAGTGTTCACCCAAAAAAAACAGATa                                                                    >  1:6686/1‑150 (MQ=14)
                                                                                     ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTTCTTCCCAACGGTCAAATCCCCTCCCCGCGGGTTTTTCCCCCTCGGGGGAGAGGTCACCCACAAACAAACGATaa                                                                   >  1:917828/1‑150 (MQ=14)
                                                                                      tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCTCCCGGCCGGTTGTTTCCAAAGTTTAAAACCCCCCCCCGCGGATTTGTCCTACTCCGGAGAGGGTTTACCCACCAACAAAAGAaaaa                                                                  >  1:1456553/1‑150 (MQ=9)
                                                                                       cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCCCTCCCCGCGGGTTTTTCCTACTCGGGGGAGTGGTTACCGACAAACAACAGATaaaa                                                                 >  1:28093/1‑150 (MQ=14)
                                                                                        ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACa                                                               <  1:1146701/151‑1 (MQ=255)
                                                                                        ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCCGCGGGTTTTTCCTACTCGGGAGAGGGTTCACCCACAAACAACAGATaaaa                                                                 >  1:603749/1‑149 (MQ=14)
                                                                                          aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaa                                                             <  2:1153615/151‑1 (MQ=255)
                                                                                          aTGGCCTTTCCCTTGATTTGATGTCTGGCAGTTTATGGCGGTCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaa                                                             <  2:814750/151‑1 (MQ=255)
                                                                                               cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTTAAAACCCCTCCCCGCGGGTTTTTCCTACTCCGGGGAGGGTTTACCCACAAACAAAAGATAaaaaaaaa                                                            >  1:1148344/1‑147 (MQ=14)
                                                                                                tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCCCCCCGCGGGTTTTTCCTACTCCGGAGAGTGTTCACCGACAAACAACAGATAAAAAAAAAAGccc                                                       >  1:1282376/1‑151 (MQ=14)
                                                                                                tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCCGCCGTTGCTTCCCAACGTTTCAAACCCCCCCCCGCGGGTTTTTTCTCCTCCGGGGGGGGTTTGCCCCCCAACAACAGATAaaaaaaaaagac                                                        >  2:1454432/1‑146 (MQ=9)
                                                                                                 ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGGTGCTTCTCAACCTTCAAAAACGCCCCCCGCGGATCTGTCCTACTCGGGGGGGGGGTCACCGACAAACAACCGATAaaaaaaaa                                                            >  2:1575074/1‑145 (MQ=14)
                                                                                                 ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACGGATAAAACAAAAGGcc                                                        <  1:802482/149‑1 (MQ=255)
                                                                                                  tCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCa                                                      <  1:1030182/150‑1 (MQ=255)
                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                     <  2:101309/150‑1 (MQ=255)
                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGt                                                    <  1:1300633/151‑1 (MQ=40)
                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTCCCCCCCACCCCCCCGGGCGGTGCTTCGCAAAGGTCAAATTCGCCCCCCGCGGAGTTTTCCCACCCGGGAGAGTGTTGACCCGCAAACAACAGATAAAAAAAAAAGcc                                                        >  1:118277/1‑147 (MQ=255)
                                                                                                      ttgatttgatGTCTGGCAGTTTATGGCGGGCGGCCTCCCCGCCCCCCTCCCGGGCGTTGCTTCGCCACGTTCTAAACCGCCCCCCGCGGGTTTTTCCTCCTCCGGGGGGGGGTTTCCCTCAAACAACAGAGAAAACaaaaaacc                                                        >  2:747326/1‑140 (MQ=2)
                                                                                                       tgatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCt                                                  <  2:1270203/149‑1 (MQ=40)
                                                                                                         atttgatGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGTTGCTTCGCAACGTTCAAAACCCCTCCCCGCGGGTTTTTCCTTCTCCGGGGAGAGTTTACCGACAAACAAAAGATAAAAAAAAAAGGCCCGTCTTcc                                               >  1:458888/1‑150 (MQ=2)
                                                                                                         atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                              <  1:516253/151‑1 (MQ=38)
                                                                                                         atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCGGGCCGTTTCTTCGCAAACGTTCAATCCCCCCCCGGCGGGATTTTCCTCCTCAGGAGAGTGTTTACCCACAAACAAACAATAAAAAAAAAAGACCAGTCTTcc                                               >  1:1618754/1‑150 (MQ=2)
                                                                                                            tgatGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCACCCCCCCGGGCGGGGCTTCGCAACGTTCAAAACCGCCCCCGCCGGGGTTTTCCTATTCGGGGGAGGGTTCACCCACCAACAAAAGCTAAAAAAAAAAGccccagcttccccc                                            >  2:1517153/1‑139 (MQ=2)
                                                                                                                   tGGCAGTTTATGGCGGGCGTGCTGCCCCCCCCCCCCCCGGCCGGTTCTTTTCCACGTTCAAATTCGCCCCCCGCCGGATTGTTCTTCTTGGGGGAGGGGTTACCGACCAACAACAACAAAAAAAAAAAGCACCGCCTTCCGCCTGCGCCtt                                    >  2:1090978/1‑151 (MQ=2)
                                                                                                                      cAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                                 <  2:36922/151‑1 (MQ=37)
                                                                                                                           tATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGTATTTGTCCTACTCGGGAGAGTGTGCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTTATCCTTTTGtttt                             <  2:557555/150‑1 (MQ=255)
                                                                                                                           tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                             <  1:1388401/150‑1 (MQ=38)
                                                                                                                               gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                        <  2:1400080/151‑1 (MQ=38)
                                                                                                                                tgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                        >  2:950476/2‑150 (MQ=38)
                                                                                                                                 gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                      <  1:803324/151‑1 (MQ=38)
                                                                                                                                      gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc                   <  1:1346141/150‑1 (MQ=40)
                                                                                                                                      gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt                  >  2:323276/1‑151 (MQ=40)
                                                                                                                                       tCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                 >  2:846037/1‑151 (MQ=40)
                                                                                                                                        ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                <  2:624464/151‑1 (MQ=40)
                                                                                                                                           gcacGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGACCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTGGATGTCggg                <  2:334136/145‑1 (MQ=255)
                                                                                                                                           gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc               <  1:1483683/149‑1 (MQ=40)
                                                                                                                                           gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa              <  2:724860/150‑1 (MQ=40)
                                                                                                                                               gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc         <  1:1155050/151‑1 (MQ=40)
                                                                                                                                               gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc         >  1:1117816/1‑151 (MQ=40)
                                                                                                                                                ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGCGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGCCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc        >  2:1375753/1‑151 (MQ=25)
                                                                                                                                                  aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt       >  1:1591921/1‑150 (MQ=40)
                                                                                                                                                   cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt       <  2:1392542/149‑1 (MQ=40)
                                                                                                                                                   cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTa      <  2:159768/150‑1 (MQ=40)
                                                                                                                                                      tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct  >  2:279597/1‑151 (MQ=40)
                                                                                                                                                      tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGCCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct  >  2:817153/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                          
ATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCT  >  NC_000913/3423145‑3423444

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