Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,287	0	G	C	42.9%	109.0 / 38.2	78	intergenic (+93/+136)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (11/33);  new base C (33/0);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 2.11e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTC  >  NC_000913/3423147‑3423436                                                                                                                                              |                                                                                                                                                           ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCAcccccc                                                                                                                                                   >  2:937760/1‑150 (MQ=255)  cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCG‑GTCCTGGCCCCCAccccccc                                                                                                                                                  >  1:181782/1‑149 (MQ=255)   tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCGCCAcccccccg                                                                                                                                                 >  2:871288/1‑148 (MQ=255)     aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGcc                                                                                                                                              <  2:1080344/151‑1 (MQ=255)        cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCTGCCCGCCCCCCTCCggg                                                                                                                                                >  2:297700/1‑146 (MQ=255)         agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                           <  2:861703/150‑1 (MQ=255)         agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGCCCCCCCCCCTCCGGGCCGGTg                                                                                                                                          >  2:97066/1‑151 (MQ=255)           agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                         <  1:302089/150‑1 (MQ=255)            gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCGTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                       <  1:572202/151‑1 (MQ=255)            gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                       <  2:60458/151‑1 (MQ=255)             cGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCCCCCGGGCggt                                                                                                                                           >  1:247501/1‑144 (MQ=37)               tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCCCCCGGCCGGTTCtt                                                                                                                                       >  2:206130/1‑148 (MQ=37)                  gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTTCCCGCCCCCCTCCCGGCCGGTGCTTCG‑cca                                                                                                                                 >  1:530158/1‑148 (MQ=255)                   aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACg                                                                                                                              <  1:850731/151‑1 (MQ=255)                    gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCTCCGGGCCCTTGCTTCG‑C‑aa                                                                                                                                >  2:60307/1‑148 (MQ=255)                      gCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCACCCTCCGGGCCGTTTCTTTG‑C‑AACGtt                                                                                                                            >  2:119806/1‑150 (MQ=255)                                cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCt                                                                                                                  <  1:28945/150‑1 (MQ=255)                                cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCCCCCGGCCGGTGCTTCC‑C‑AACGTTCaaaatcgccc                                                                                                                 >  2:400388/1‑144 (MQ=21)                                cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GCCCCGCCCCCCACCCCCCCGGCCCGTCCTTTG‑C‑AAAATTTCAATCCGcac                                                                                                                 >  1:794361/1‑149 (MQ=255)                                  cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTccc                                                                                                               <  2:326169/151‑1 (MQ=255)                                  cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTTG‑C‑CAAGTTCAAAAcccccccc                                                                                                               >  1:604006/1‑151 (MQ=21)                                    ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCcgg                                                                                                             <  2:427075/151‑1 (MQ=255)                                     gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCcggc                                                                                                            <  2:186566/151‑1 (MQ=255)                                        aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGa                                                                                                        <  2:256021/151‑1 (MQ=255)                                          ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCTGCCCCCCCCCCCCCCGGCCGCTGCTTCG‑C‑AACCTTCAAAAACGCTCCCGGC‑GGGtt                                                                                                      >  2:255497/1‑151 (MQ=255)                                             ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTACTGCCCCCCCCCCCCCGGGCCGGTGCTTCG‑C‑AAAGTTTAAAACCGCCCCCCGC‑GGGttt                                                                                                     >  2:215412/1‑149 (MQ=18)                                              tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCCCCCTCCCCCCCCCCGGGCGTTGGGTCG‑C‑AACGTTCTAATCCCCTCCccgc‑gggctt                                                                                                     >  2:925090/1‑144 (MQ=18)                                               agttGAATGTTAGTTCAAAAAGCAAAAAGGCCATCCTTTCGGAGGGCCTTTGGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTCCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTg                                                                                                    <  1:879049/146‑1 (MQ=255)                                                 ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCTGCCCCCCACCCTCCCGGCCGCTGCTTCC‑C‑AACGGTCAAATCCGCTCCCCGC‑GGGTTTTTcc                                                                                                 >  2:647593/1‑149 (MQ=255)                                                 ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGA‑GGCCCGCCCCCCCCCCCCCCCGCCGGTGCTTCG‑C‑TACCTTTAAAACCCCCCCccgccggct                                                                                                       >  1:170533/1‑139 (MQ=255)                                                  tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTTG‑C‑AACCTTTCAATCCGCCCCCCGC‑GGGTTTGttc                                                                                                 >  1:481821/1‑146 (MQ=255)                                                    aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑c                                                                                          <  2:90991/151‑1 (MQ=255)                                                     aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCCCCCGGCCCGTGCTTCT‑C‑AACCGTCAAATCCGCCCCCCGC‑GGGTTTTTCC‑Ttct‑c                                                                                          >  1:495446/1‑150 (MQ=17)                                                       gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑Cggg                                                                                       <  1:1162742/151‑1 (MQ=255)                                                              cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCTTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAgtgt                                                                                <  1:994215/151‑1 (MQ=255)                                                                aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTCCCCGCCCCCCCCCCGGCCGGTTTTTTG‑C‑AACGCTCAAATTCGCCCCCCCC‑CGCGTTGTCC‑CTCT‑Cgggggggggtt                                                                               >  1:50635/1‑142 (MQ=12)                                                                 aatgCAAAAAGGCCATCCGTCAGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGC‑GTCCTCCCCCCCCCCCCCCCGGCCGGTGCTTCT‑C‑AACCGTCCAATCCACCCCCCGC‑GGGTTTTTTC‑Ttct‑cccggaagggtt                                                                               >  1:322159/4‑138 (MQ=2)                                                                 aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCa                                                                             <  2:939661/151‑1 (MQ=255)                                                                  aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCAc                                                                            <  1:581185/151‑1 (MQ=255)                                                                   aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCC‑AAATTTCAAA‑CCCCTCCCCGC‑GGGGTTTTTC‑TTCT‑CCGGAGAGGGtttcc                                                                            >  2:486432/1‑147 (MQ=9)                                                                      aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCCCCCGGGCCTTGCTTCG‑C‑CACGGTCAAAACCGCCCCCCGC‑GGGTTTGTTC‑TACT‑Ccgggggaaggtcacc                                                                           >  2:286385/1‑136 (MQ=14)                                                                       aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTTG‑C‑AACGTTTAAAACCCCCCCCCGC‑GGGTTTTTCC‑TTCT‑CGGGggaggggttacccac                                                                        >  2:1103629/1‑139 (MQ=17)                                                                           ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGacaaaca                                                                   <  2:47186/151‑1 (MQ=255)                                                                            gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCCGCCCCCCACCCTCCCGGCCGGTGCTTCG‑C‑AACGTTCAAAACCGCCCCCCGC‑GGGTTTTTCC‑TACT‑CCGGAGAGGGTTCACCGacaaa                                                                     >  2:637222/1‑148 (MQ=255)                                                                                tCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGCCCCCCCCCCTCCCGGGCGGTGCTTTG‑C‑CAACTTTAAATACCCCCCCCGC‑GGGATTTTTC‑TACT‑CGGGGGGGGGTTTAccccccaaaaa                                                                  >  1:1134154/1‑138 (MQ=9)                                                                                 ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑CGTCTCCCCCCCCCCCCCCCGGCCGTTTCTTCT‑C‑AACCCTTAAAACCCCTCCCGGCGGGATTTTTCC‑TCCC‑CGGG‑GAGGGGTCCCCCacaaaca                                                                   >  1:1187511/1‑145 (MQ=255)                                                                                    ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCG‑C‑AACGTTTAAATCCGCCCCCCGC‑GGGTTTTTTC‑TACT‑CCGGAGAGTGGTTACCCGCAAACAACAAATaa                                                           >  2:940838/1‑150 (MQ=14)                                                                                     tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATaaa                                                          <  2:77149/150‑1 (MQ=255)                                                                                      cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAAc                                                        <  1:683438/151‑1 (MQ=255)                                                                                       ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCCCCCGGCCCTTGCTTCG‑C‑AAAGTTCAAAACCGCTCCCCGC‑GGGTTTTTTC‑TACT‑CTGGAGAGTGTTTACCCACAAACAAAAGATaaaa                                                         >  2:396685/1‑149 (MQ=14)                                                                                         aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCTGGCCCCCAACCCCCCGGCCGTTGGGTTG‑C‑GACACTCAAAATCCCCCCCCGC‑GGAGTTTTCC‑TCCA‑CGGGAGAGAGGTCACCGACAAAAAAACGATAaaaaaa                                                      >  2:339316/1‑150 (MQ=255)                                                                                         aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCTGCCCCCCACCCTCCCGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCCGC‑GGGTTTTTTC‑TACT‑CTGGAGAGGGGTCACCGACAAACAACAGATAaaaaaa                                                      >  2:928047/1‑150 (MQ=14)                                                                                           ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaaa                                                    <  2:792059/150‑1 (MQ=255)                                                                                            gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCCCCCGGCCGGTGCTTCG‑C‑AACGTTCAAATCCGCCCCCCGC‑GGGTTTTTCC‑TACT‑CCGGAGAGGGTTTTCCCACAAACAACAGAAAAAACaaaa                                                    >  2:58797/1‑149 (MQ=14)                                                                                              cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTGCCGCCCCCCACCCCCCCGGCCGTTTCTTCG‑C‑CACCTTTTAATCCCCTCCCCGC‑GGGTTTGTCC‑CTCT‑CGGGGGAGAGTTTTCCGCAAAACAACAGAaaaaaaaaaaa                                                   >  1:1025093/1‑147 (MQ=9)                                                                                               tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTTCCCCCCACCCTCCCGGCCCTTTCTTTT‑CAAAGGTTCAAA‑CCGCCCCCCGC‑GGGTTTTTCC‑TTCTCCGGGGGGGGGTTTACC‑ACAAACAAAAGATAaaaaaaaaaggcc                                               >  1:910997/1‑146 (MQ=255)                                                                                                ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTTCCCCCCCCCCCCCGGG‑CGCTGCTTTG‑C‑AACGTTCAAAATCGCCCCCCGC‑GGGTTTGTCC‑TCCT‑CGGGAGGGGGGTTTCCGTCCAAAAAAAGAaaaaaaaaaa                                                    >  1:472962/1‑145 (MQ=255)                                                                                                ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCCCCCCCCCGGGCGTTGCTTTT‑C‑CAACCTCAAATCCCCCCCCCGC‑CGGTTTTTCCTTCCC‑CGGGGGGGGGTTC‑CCCACCAACAACAGAAAAAACaaaaaaccc                                               >  1:684497/1‑145 (MQ=9)                                                                                                 tCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCCGCCCCCCCCCCTCCTGGCCGTTGCTTTG‑C‑CACGTTCCAATCCACCCCCGCC‑GGGCTTGTCC‑CACT‑TCCGGGGGGGTTCACCGACCAAAAACACAAAAAACaaaa                                                    >  2:449387/1‑144 (MQ=255)                                                                                                  cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGt                                            <  1:736020/151‑1 (MQ=40)                                                                                                   gCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCGCCACCCTCCGGGGCGGTGGCTCG‑C‑AACCCTCAAAACCCCTCCCCGG‑GGATTTTTTC‑TTCT‑CGGGGGAGGGTTTACCCCCAAAAAAAAAATAaaaaaaaaagaccc                                              >  2:939389/1‑142 (MQ=255)                                                                                                        atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                      <  2:247501/151‑1 (MQ=38)                                                                                                        atttgatGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCCCCCGGCCGCTGCTTCGCC‑AACG‑TCAAATCCGCTCCCCGC‑GGGTTTGTCC‑TACT‑CCGGGGAGGGGTTACCCACCAACAAAAGATAAAACaaaaaaccccgactttc                                       >  1:1147977/1‑137 (MQ=255)                                                                                                           tgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCTCCCGGGCGGTGCTTCT‑C‑CACGTTTAAAATCCCTCCCCGC‑CGGGTTTTCC‑CCCT‑CGGGAGAGGGGTTCCCCACAAACCACCGAAAAAAAAAAAAGCCCACTCTTcc                                       >  1:234842/1‑147 (MQ=2)                                                                                                            gatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCG‑C‑AACGGTCAAATCCGCCCCCGGC‑GGATTTTTCC‑TACT‑CCGGGGAGGGTTTACCGACAAACAACAAATAAAACAAAAAGCCCAGTCTTcc                                       >  2:351224/1‑146 (MQ=11)                                                                                                             atGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTg                                  <  2:693075/150‑1 (MQ=38)                                                                                                                   ggCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                             <  2:762322/149‑1 (MQ=35)                                                                                                                         ttATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGttt                      <  2:1183264/150‑1 (MQ=38)                                                                                                                          tATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                     <  1:1134330/150‑1 (MQ=38)                                                                                                                           aTGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                    <  2:488880/150‑1 (MQ=38)                                                                                                                                gggc‑gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt              <  1:1248/151‑1 (MQ=38)                                                                                                                                     gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt          <  2:768498/151‑1 (MQ=40)                                                                                                                                     gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt          >  2:1163343/1‑151 (MQ=40)                                                                                                                                          gCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc       <  1:874330/149‑1 (MQ=40)                                                                                                                                          gCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa      <  2:1187520/150‑1 (MQ=40)                                                                                                                                           cccGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa      >  2:78397/1‑149 (MQ=40)                                                                                                                                           cccGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa      <  1:78397/149‑1 (MQ=40)                                                                                                                                             cGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc  >  1:523253/1‑151 (MQ=40)                                                                                                                                              |                                                                                                                                                           CCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCG‑C‑AACGTTCAAATCCGCTCCCGGC‑GGATTTGTCC‑TACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTC  >  NC_000913/3423147‑3423436

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