Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,256,069	0	G	C	34.0%	92.3 / 15.9	50	P133R (CCG→CGG)	plsB	glycerol‑3‑phosphate 1‑O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (5/28);  new base C (17/0);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 2.67e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTC  >  NC_000913/4255823‑4256305                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                tCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGGTAACTTTGCCTTTTT‑TACACCCC‑Cg‑cgc                                                                                                                                                                                                                                            >  1:568729/1‑251 (MQ=255)      aCGAGCCGGCAAACGTGGGCCTCGGGCGCCCAGACGTTGCCGGCAAAAGTCCATCCCCGCCACGCGCGCCGTTTTCGGAGCGATGGTTTTATCGGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGCCGGCGAGAAACGCACAAACCGGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACGTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGAc                                                                                                                                                                                                                                       <  2:583946/251‑1 (MQ=255)        gAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCCTTTT‑TACGCCCC‑CG‑CGCGCGAccc                                                                                                                                                                                                                                     >  1:689255/1‑250 (MQ=255)        gAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGGTTACCTCGCCCTTTT‑TACTCCCC‑Cg‑cgcccgcccc                                                                                                                                                                                                                                     >  1:249532/1‑244 (MQ=255)         ggCAGGCAAAGGGGGCCCCACGGCAGCCGGGCGTGGACGGGCAAAGTGCATCCGCCCCACGCGCGCCAGTTTGGGGGCGATAGTTTTAGCCGGGCCGTTTTCATCCGCCATACGGCGCAGCGAAACGGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCaa                                                                                                                                                                                                                                    <  2:41822/250‑1 (MQ=255)             ggCAAACGTGGGCTTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCGGAGCGATAGTTTAATCCGTGCCGTGTCCATCCGCCATACGGCGCAGCGAAACGGACGGCGAGAAACGCACAAAGCTGCCGCGACTGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACAt                                                                                                                                                                                                                                <  1:680335/251‑1 (MQ=255)                aaaCGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGGTCCCTTTCCCTTTTT‑TTCGCCCC‑CG‑CGCGCGAcccaaaaaccc                                                                                                                                                                                                                             >  1:718298/1‑242 (MQ=255)                  aCGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGGTTACTTCTCCTTTTTTCAC‑CCCC‑CG‑CGCGCGCCCAAACaaccccc                                                                                                                                                                                                                           >  1:152413/1‑245 (MQ=255)                    ggggcgctccGGCAGCCAGCCGTGGACGGGCAAGGGGCAGACGCGCCCCGCGCGCCAGTTTCTGAGCGTAATTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACGGACGGCGAGAAACGCACAAAACGGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGac                                                                                                                                                                                                                         <  2:536519/242‑1 (MQ=255)                    gTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAATGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGAAGCGGCGGGTTGACGTCGCCGTTTT‑CAAGCCCC‑CG‑CGCCCGACCCACCGTCACCGac                                                                                                                                                                                                                         >  1:37019/1‑251 (MQ=255)                            aCGGCAGCCAGACGTTGACGGGCAAAGTGCATCCGCGCCCCGCGCGCCAGTTTCGGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGc                                                                                                                                                                                                                   <  1:147880/249‑1 (MQ=255)                                       aCGTTGACGGGCAAAGTGCATCCGCGCCACGCGCCCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑t                                                                                                                                                                                                     <  1:13421/251‑1 (MQ=255)                                                       gccaagcgcgccgcgcGCGCCAGTTTCGGGGGGATAGTTTTATCCGGCCCGTGTCCACCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCAAAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCt                                                                                                                                                                                    <  2:79540/245‑1 (MQ=255)                                                                cccccccccgcCGGTTTCTGAGCAGTAGTTTTACCCGGGCCGTGTTCACCCGCCAAACGGCGCAGCAAAACTGACGGCGAGAAACGCACAAAACTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCaa                                                                                                                                                                           <  2:234104/243‑1 (MQ=255)                                                                   cggcgcgccgtttgcggAGCAATAGTTTTAGCCGAGCCGTGTACATCCGCAAGACGGCGCAGCGAACCTGACGGCGAGAAACGCCCAAAACTGTCGCGACCGAGCCACAGCACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTC‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATc                                                                                                                                                                        <  1:113928/235‑1 (MQ=255)                                                                          ccAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCCCTTCTCCTTTTT‑TACGCCCC‑CG‑CGCGCGACCAAACATCACCCACACCGGCCCCATCTCCACC‑TCCCGATTTTGGGTT‑CTACGCTGCCAATCGAGAt                                                                                                                                                                   >  1:43788/1‑249 (MQ=255)                                                                            aGTTTCTGGGCGATAGTTTTATCCGTGCCGTTTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTc                                                                                                                                                               <  1:585057/251‑1 (MQ=255)                                                                            aGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGGTTACTTCGCCTTTTT‑CACCCCCC‑CG‑CGCGCGACCAAACAACACCGACACCGGCCCCAACTGCCCA‑CCCCGATTT‑GGGGTTCTACGGGGCAAAACgaggaaagc                                                                                                                                                               >  1:382224/1‑245 (MQ=255)                                                                             gTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCCTTTT‑CACGCCCC‑CG‑CGCGCGCCCCAACAACACCCCCCCCGGCACCACCTCCACA‑ACCCGATTT‑TGGGTGCTACCGGGCAAATccggaaaagcg                                                                                                                                                              >  1:273367/1‑241 (MQ=255)                                                                                 cggAGCGATAGTTTTACCCGTGCCGTTTTCATCCGCCATACGGCGCAGAGAAACTGAGGGCGAAAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTGCGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACAGCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTgg                                                                                                                                                           <  2:622191/248‑1 (MQ=255)                                                                                   gAGCGTTAGTTTTCTCCTTGCAGTGTTCACCCGCCATCCGGCGCAGCGAACCTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGa                                                                                                                                                          <  1:134461/249‑1 (MQ=255)                                                                                   gAGCGATAGTTTTATCCGTGCCTTGTTCATCCGCCATCCGGCGCAGCGAAACTGACGGCGGGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGa                                                                                                                                                          <  1:312364/249‑1 (MQ=255)                                                                                    aGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGAACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGATCACCTCTCCTTTTT‑CACACCCC‑CC‑CGCGAGGCCAAACACCACCGAAACTGGGACCATCTCCACAACCCAGATTG‑GGGTT‑CTACGGGTCAAAACAAGATAATCGTgggac                                                                                                                                                        >  1:127339/1‑247 (MQ=255)                                                                                        atatttttATCGGTGCCTTGTTCATCCGCCATCCGGCGCCGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCCCTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCtt                                                                                                                                                   <  2:552609/247‑1 (MQ=255)                                                                                        aTAGTTTTACCCGTGCCGTGTTCAGCCGCCATCCGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCtt                                                                                                                                                   <  1:380863/251‑1 (MQ=255)                                                                                            ttttACCCGTGCCGGGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGGGCCCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAAt                                                                                                                                                <  1:653526/250‑1 (MQ=255)                                                                                             tttacccgggcggtttactaccGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCAAAAACCTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CCCGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATa                                                                                                                                               <  2:549860/230‑1 (MQ=255)                                                                                                 tCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTTGCCTTTTT‑TACCCCCC‑CG‑CGCGCGACCAAACaac                                                                                                                                                                                                                               >  1:577860/1‑166 (MQ=255)                                                                                                      gcctgtttaacccGCCATCCGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTcttctt                                                                                                                                     <  1:438966/240‑1 (MQ=255)                                                                                                               aTCCGCCATACGGCGCAGCGAAACTGACGGCGATAAACGCACAAAACGGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTGCGCCGTTAAGCATACGCAGCGGCGGGTTCCCTTCGCCTTTTG‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCgtgt                                                                                                                            <  1:84808/251‑1 (MQ=255)                                                                                                                 ccgccaaccGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACGGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACGTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTaa                                                                                                                          <  2:265935/243‑1 (MQ=255)                                                                                                                      atgcGGCGCAGCGAACCTGACGGCGAAAACCNCACAACACGGTCGCGACGGAGCCACAGTACAGCGAAAAATTGCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGt                                                                                                                     <  2:338619/248‑1 (MQ=255)                                                                                                                                      cTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCCCCTTTTT‑CTCGCCCC‑CG‑CCCCCGGCCCAACATCAACCCCACTCGCACCATCTTCCCC‑TCCACATTT‑TGGGTGCTTCCGGGCACATCTAGACAGTCGTGGGACAGCTTAATAGACCCTTATTTCTGCGGGTAATAGGTGGacacacaccgcccc                                                                                                     >  2:442740/1‑243 (MQ=255)                                                                                                                                        gACGGCGAGAAACTCACAAAACGGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTGCGGCGTGTAATAGGTGAACACACGCGGCccgcc                                                                                                   <  2:249888/251‑1 (MQ=255)                                                                                                                                              gagaAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGa                                                                                                                                                                                                <  2:352689/153‑1 (MQ=255)                                                                                                                                              gagaAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTT‑TACGCCCC‑CG‑CGCGCGCCCAAACAACACCGCCACTGGCACCACCTGCCCA‑TCCAGa                                                                                                                                                                                                >  1:352689/1‑153 (MQ=255)                                                                                                                                               agaAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGGTTACCTTGCCTTTTT‑CACGCCCC‑CG‑CGCCCGACCAAAAATAACCGACACCGGCACCCTCTGCACA‑TACAGATTT‑GGGGTGCTACGGGGCAAATCGGGGAAGTCGTGGGACAGCTTAAAAGACTCTTCTTTCTGCGCGTAATAAGTGAACACACGCGGCccgccccgaa                                                                                              >  1:274442/1‑244 (MQ=255)                                                                                                                                                       aaaaacctgcggcgcccGAGCCACAGTAAAG‑GAAAAATTTTTTGTCCCCCGTAAAGCATCCGCGGCGGCGGGTTCACTCCGCCTTTTC‑CCCGCCCG‑GG‑CGCGCGACCAAACACCCCCGACACTGGCCCCATCTGCACA‑TCCAGTTTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTAAACCCACGCGGCCCGCCGTGAATGAACACATa                                                                                    <  1:362497/236‑1 (MQ=255)                                                                                                                                                       aaaaaacgcgcgagcccGACCCAAAGTCCAGCAAAAAATTTCT‑GTCCGCCGTAAAGCATACGCAGCGGTGGTTTCACTTCGCTTTTTC‑CCCCCCCC‑GG‑CGCCAGCCAAAACATCCCCGCCCCGGGCACCATCTGCACC‑ACCAGTTTG‑GGGTTGCTCCGGGGCAAATCGAGATAGTCGTGGAACAGTTTAATAGTTTTTTCTTCCGGCGTGTAATAGGTGACCACACGCGGCCCGCCGTGAATGAACACATa                                                                                    <  1:6218/236‑1 (MQ=255)                                                                                                                                                                  gcgACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACCTCGCCTTTTT‑TACGCCCC‑CG‑CGCGCGGCCAAACAACACCCACCATGGCACCATCCCCACA‑TCCAGATTT‑GGGGTGCTACCGTGCAAAACGAGATAGTCGTGGGACAGCTTAAATGACTCTTCTTTCTGCGTGGAATAGGTGAACACACGCGGCCCGCCCGGAATGGACACATaaggcgggagt                                                                          >  1:399194/1‑240 (MQ=255)                                                                                                                                                                             acaGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCTCCTTTTT‑TCCCCCCC‑GG‑CGGGCGGCCCAACAACACCCCCACTTGCACCCTCTGCCCA‑TCCAGATTT‑TGGGTGCTACgg                                                                                                                                                                                >  2:178755/1‑137 (MQ=255)                                                                                                                                                                             acaGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACgg                                                                                                                                                                                <  1:178755/137‑1 (MQ=255)                                                                                                                                                                                                   gTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCTCCTTTTT‑TACGCCCC‑CC‑CGCGCGACCAAACAACACCCACCCCGGCCCCCCCTGCCCC‑ACCCGATTT‑TGGTTTCTACGGGGCAAAACGAGATAAGCGTGGGACAGCTTAATAAACTCTTCTTTCGGCGTGTAAAAGGTGAACACACGCCGCCCGCCGTGAAAGAAAACATAACGCGGCAGGAGCGTGCCGGCGactttcaaccggcctcaacgg                                         >  1:96885/1‑231 (MQ=255)                                                                                                                                                                                                       gCCGTTAAGCATACGCAGCGGCGGGTTCCCTTCGCCTTTTT‑CAAGACCC‑CG‑CGCGCGAACAAAAATCACCGGCACTTGGACCCTCTGCACA‑TCCAGATTT‑GGGGTGTTCCCGTGCAAATCGAGATAGGTGAGGGAAAGCTTAATAGACTCTTCTTTCGGGGTGTAATAGGTGAACACACGCCGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGGTTTTCAGCGGCTCTAACggggc                                       >  1:155038/1‑247 (MQ=255)                                                                                                                                                                                                        ccgtaaacaaaccGCAGCGGCGGTCTCCCTTC‑CCTTTTT‑CCCCCCCCCGC‑CGCGCGACAAACCACCCCCGACCCTGGCACCATCAGCAGA‑CCCAGATTG‑GGTTTGCTACGGTGCAAATCGAGATAGTCGTGGACCAGCATAATGGACTTTTCTTGCGGCGTGTAATCGGTGAACCCACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCatg                                    <  1:160281/239‑3 (MQ=255)                                                                                                                                                                                                                 aTACGCGGCGGCGGTTTCTCT‑CCCCTTTTT‑CCCCCCCC‑GGGCGCGCGCCAAAACATCCCCGACACTGGCCCCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGAAACAGCTTAATAGTCTCTTCTTTCGGCGTGTAATGGGTGAACCCCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCAtg                           <  2:382224/251‑1 (MQ=255)                                                                                                                                                                                                                              gggtCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGAGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTgg              >  1:120688/4‑251 (MQ=255)                                                                                                                                                                                                                                  cACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCCGGCAAGTCCTGTGCCAGACACTGGgcgc          >  2:264427/1‑251 (MQ=255)                                                                                                                                                                                                                                    cTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGgcgcgc        >  1:131875/1‑251 (MQ=255)                                                                                                                                                                                                                                           ttttt‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTc  >  1:578252/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                TCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCCAGTTTCTGAGCGATAGTTTTATCCGTGCCGTGTTCATCCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCT‑GTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCC‑GG‑CGCGCGACCAAACATCACCGACACTGGCACCATCTGCACA‑TCCAGATTT‑GGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTC  >  NC_000913/4255823‑4256305

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