Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	974,449	0	A	G	28.8%	96.3 / 18.7	59	D44G (GAT→GGT)	mukF	chromosome partitioning protein MukF
Reads supporting (aligned to +/- strand):  ref base A (15/27);  new base G (17/0);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 3.57e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATGAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGCGTCT  >  NC_000913/974209‑974692                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        tATGAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGGGCGGCTGGGTGGGgaggtg                                                                                                                                                                                                                                               >  1:375756/1‑248 (MQ=255)    gAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGt                                                                                                                                                                                                                                            <  1:80868/251‑1 (MQ=255)        cATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATCTTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTTC‑CAC‑GCTCAACGGCGAGCGGCTTGGTGGGGAGATggtcga                                                                                                                                                                                                                                          >  2:219903/1‑244 (MQ=255)           tACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAgga                                                                                                                                                                                                                                      <  2:346071/250‑2 (MQ=255)                gAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGGGCGGCTGGGTGGTGGGAAGAGTGGAGGCGGGCtg                                                                                                                                                                                                                                >  1:421045/1‑251 (MQ=255)                    gaGAAACGCGTTATTGCCGGCAGGAACCGTATATTACCCTGGGGCGGTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGTGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGTTGAGATGAGTGAATGCGAGCtggtg                                                                                                                                                                                                                             <  2:469831/249‑1 (MQ=255)                      gaagcgcgTTATTGCCGTCAGGAACCGGATAGTACCCTGGGGCGTGATATTCATGTCACCTCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAGCTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCTTCTCGCTGCCGGTAGACCGGCTCTCTTTTCGGCTGGCGGTTGC‑GAT‑TCTGAACGGCGAGCGTCTGGATGGTGAGATGTGTGAAGGCGAGGGGGTGGAt                                                                                                                                                                                                                          <  2:296000/247‑1 (MQ=255)                       aaaCGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGGGAAGGCGAGCtggggggtg                                                                                                                                                                                                                         >  1:42131/1‑245 (MQ=255)                        atcgcgTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGc                                                                                                                                                                                                                        <  1:444212/249‑1 (MQ=255)                            cgTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCCGGACGGCGAGGGTCTGGGTGGGGAGATGGGGGGAGGCGCGCtggggggggcattc                                                                                                                                                                                                                    >  2:110138/1‑240 (MQ=255)                             gTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTcc                                                                                                                                                                                                                   <  1:110138/251‑1 (MQ=255)                              ttATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGGCTGGATGGGGAGATGAGTGAAGGCGAGCTGGTGGATggctttcc                                                                                                                                                                                                                  >  1:233861/1‑244 (MQ=255)                                   gCCGTCAGGCACCGTATATTAACCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGGTTC‑CACCGC‑GAACGGCGGGCGCCTGGGCGGTTGGATGAGTTAAGGCGAGCTGGAGGATGCATTCCGCCAtt                                                                                                                                                                                                             >  1:367651/1‑250 (MQ=255)                                          gtaaCCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTTAAGGCGAGCTGGTGGATGCATTCCGCCATgtgagtga                                                                                                                                                                                                      <  2:337357/249‑1 (MQ=255)                                                 tataTTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGCGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGG‑CAC‑GCCTAAAGGCGAGCGTCTGGGTTGGGGGGTGAGTGAAGGCGGGCGGGGGGAGGCATTCCGTCATGTGTGTGGTGCG‑ttt                                                                                                                                                                                              >  2:379303/1‑251 (MQ=255)                                                        cccTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGGACGGGGAGCGTCTGGGTGGTGAGAAGGGGGGAGGGGAGGTGGTGGATGCGTTCCGcccggggggggatggggttgggcaa                                                                                                                                                                                         >  2:179252/1‑227 (MQ=255)                                                               gCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGGGAGATGAGGGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCg                                                                                                                                                                                 >  1:14452/1‑250 (MQ=255)                                                                 gTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAACGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGCTGCACAG‑GC‑GAACGGCGAGCGGCTGGGTGGGGGGGAGCGGGAAGGAGAGCCGGAGGATGGACTCCGCCATGTGAGTGCTGCggctgaggcaaacagcgaa                                                                                                                                                                                >  1:234096/1‑232 (MQ=255)                                                                 gTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGGGAGGTGAGTGAAGGCGAGCTGGGGGATGCATTTCGCCAAGGGAGTGATGCG‑TTTGGGCAAACCCGCGaa                                                                                                                                                                               >  1:222304/1‑250 (MQ=255)                                                                         cATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGG‑CAC‑GCTGAACGGCGAGCGCCTGGATGGGGGGATGAGTGAAGGCGGGGTGGGGGGTGCCTTCCGCCATGTGAGTGAGGCG‑GTTTAGCAAAACAGCGGAACCATCggg                                                                                                                                                                      >  1:432222/1‑250 (MQ=255)                                                                          aTGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGGGAGATGAGTGAAGGCGAGCTGGGGGATGCATTCCGCCATGTGAGGGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCg                                                                                                                                                                     >  1:205014/1‑251 (MQ=255)                                                                                ccGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGAGGGCGCCTGGGTGGGGAGATGAGTGAAGGCGCGCTGGGTGATGGATTCCCCCCTGTGAGTGATGCG‑TTTGTGCAAACCAGCGAAACCATCGGCGTgggcg                                                                                                                                                               >  1:395874/1‑251 (MQ=255)                                                                                        aaaCCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGGGAGTGATGCG‑TGTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑Ac                                                                                                                                                     <  2:432222/251‑1 (MQ=255)                                                                                                      ggATAAAGTATGAGTGAATTTTCCCAGAAAGTCCCCGACCTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATTCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACAt                                                                                                                                       <  2:322080/251‑1 (MQ=255)                                                                                                        aTAAAGTATGAGAGCATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGCCTCTCTTTTCTGCTGGCGGTTTC‑CAT‑GCTGAACTGCGAGCGGCTGGATGGAGTGCTGGGTGCAGGAGCGCCGGTTGCAGCATTCCGCCATGTGGGTGATTCGGTTTG‑GCAACACAGCGAAAAAATCGGGGTGCGCGCCCAATACACGCGacacgacaagg                                                                                                                                       >  2:382230/1‑241 (MQ=255)                                                                                                          aaaGTATGAGTGAATTTTCCCAGACAGGCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTGTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAGCCTGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGt                                                                                                                                    <  1:431811/250‑1 (MQ=255)                                                                                                                     gAATTTTCCCAGACAGTCCCCGCACTGTTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCCCCCGGGTGGGGGGATGAGTGAAAGAGAGCTGGGGGCTGCACTCCGCCATGTGAGGGATGTG‑TTTTAGGAAACAACAGAAACCATCGGGGGGCGAGCC‑CCTA‑ACACGATCAAAGACATGGTGC‑gctaaagca                                                                                                                        >  2:293000/1‑242 (MQ=255)                                                                                                                     gAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGGGCTGGGGGATGCATTCCTCCATGGGGGGGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATT‑ACGCGATCAACCACATGGTGC‑GTCAACGTct                                                                                                                       >  1:276496/1‑251 (MQ=255)                                                                                                                         tttCCCAGACAGTCCCCGACCTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGGGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTctgct                                                                                                                    <  2:386853/250‑1 (MQ=255)                                                                                                                             ccAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGGGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑AAGCGATCCAAGACATGGGGA‑GTCAACGGCTGCTGAAcc                                                                                                               >  1:99869/1‑251 (MQ=255)                                                                                                                                            aCTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTTC‑CAC‑GCTGAACCGCGGGCGCCTGGATGGGGAGATGAGTGAAGGCGAGCTGGGGGATGCATTCCGCCCTGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑CATA‑AAGCGATCAACCACATGGTGC‑GGCAACGTCTGCTGGAACGCtttttcagccagg                                                                                                >  1:55727/1‑241 (MQ=255)                                                                                                                                                  tGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCg                                                                                           <  1:100278/250‑1 (MQ=255)                                                                                                                                                       gggCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGTTGC‑GTCAGCGGCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAggg                                                                                     <  1:455783/251‑1 (MQ=255)                                                                                                                                                         gCCAGACAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CCC‑GCTGGACGGCGAGGGTCTGGGTGGTGGGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAAACCGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCCACGACCTGGTGC‑GTCAACGCCTGCTGAACCGCTTTACCAGCCAGCAGGCgggaggga                                                                                    >  1:144208/1‑250 (MQ=255)                                                                                                                                                         gCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGGGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGaa                                                                                   <  1:219903/251‑1 (MQ=255)                                                                                                                                                         gCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGGGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAAAGACATGGTGC‑GTCAACGGCTGCTTAACCGCTTTTCCCACGAGCAGGCgggaggga                                                                                    >  1:51373/1‑250 (MQ=255)                                                                                                                                                            aGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCCACTCTCTTTTCTGCTGGCGGTTGA‑CCC‑GCCGAACGGCGGGCGTCTGGGTGGTGAGATGAGTGAAGGAGAGGTGGTGGGTGCATTCCGCCATGTGTGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGGGCGCGCC‑CAAA‑AAGCGGTCAAAGGAATGGGGG‑GGCGACGACTGGTGGACCGCCTTACCAGCGAGCAGGCGGGAGGGAACGc                                                                                >  1:383674/1‑251 (MQ=255)                                                                                                                                                                  aaTGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGT‑CAC‑GCTGAACGGCGAGCGTCCGGATGGTGAGGTGAGTGAAGGCGAGCTGGTGGGGGCATTCCGCCATGTGAGTGATGCG‑CTTGAGCACACCAGCGAAACCATCGGCGTGCGGGCC‑AAAA‑ACGCGATTTACGACATGGGGCCGACA‑CGTCCGCTGAACCGCTTTACCCGCGCGACGGCGGGAGGGAACGCAAttt                                                                           >  1:268028/1‑250 (MQ=255)                                                                                                                                                                     gACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTAcc                                                                        <  1:459572/250‑1 (MQ=255)                                                                                                                                                                     gACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCGGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTAcc                                                                        <  2:14452/250‑1 (MQ=255)                                                                                                                                                                         tctcCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGGTGGGGAGATGAGTGAAGGCGAGATGGTGGATGCATTCCGCCATGTGGGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGGGC‑GTCAACGTCTGCTGGACCGCTTTACCACCGAGCAGGCGGAAGGGAACGccaattaacctctt                                                                   >  1:274755/1‑238 (MQ=255)                                                                                                                                                                            ccATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGGTGGGGAGATGAGTGAAGGCGAGCTGGGGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑AAGCGATCAACGACATGGTGC‑GTCAACGGCTGCTGGACCGCCTTTCCAACGAGCCGGCGGGAGGGGACGCAATTTACCGGCTGAcg                                                                >  1:244673/1‑251 (MQ=255)                                                                                                                                                                             cATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGCGAAGGCGAGCTGGTGGATGCATTTCGCCCTGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGGGCC‑AATA‑ACGCGATCAACGAAATGGGGC‑GTCAACGTCTTCTGAAACGCttttacag                                                                                                     >  1:150275/1‑208 (MQ=255)                                                                                                                                                                             cATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAg                                                                                                     <  2:150275/213‑1 (MQ=255)                                                                                                                                                                                  cgcggCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCtc                                                          <  2:421045/247‑1 (MQ=255)                                                                                                                                                                                  cGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGAcgcc                                                              <  1:179252/247‑1 (MQ=255)                                                                                                                                                                                  cGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCtc                                                          <  1:383460/251‑1 (MQ=255)                                                                                                                                                                                  cGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCAGGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGAcgcc                                                              <  2:72270/247‑1 (MQ=255)                                                                                                                                                                                     tGCCGGTAGACCGACTCGCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAAc                                                                                                                           <  2:386260/183‑1 (MQ=255)                                                                                                                                                                                     tGCCGGTAGACAGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGGTGGGGAGATGAGTGAAGGCGAGCTGGGGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGGTCCACGCCCAGGTGC‑GTCAAc                                                                                                                           >  1:386260/1‑183 (MQ=255)                                                                                                                                                                                        cGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCtcggcatc                                                    <  1:1472/251‑1 (MQ=255)                                                                                                                                                                                              aCCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGGCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCAtt                                              <  1:291599/251‑1 (MQ=255)                                                                                                                                                                                                           tCTGCTGGCGTGTGC‑CAC‑GCTGAACGGCGGGCGTCTGGAGGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TGTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACtactac                                  <  2:394737/250‑1 (MQ=255)                                                                                                                                                                                                                  gCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCa                          <  2:375756/251‑1 (MQ=255)                                                                                                                                                                                                                            aC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCCACGTCTGCTGAACCGCTTTACCAGCGAGCCGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCCCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGtt                 >  1:297915/1‑251 (MQ=255)                                                                                                                                                                                                                               gCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGtttt               >  1:99715/1‑251 (MQ=255)                                                                                                                                                                                                                                   aaCGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGTTTTCTAc           <  1:379303/251‑1 (MQ=255)                                                                                                                                                                                                                                         gAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGTTCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTCCCAGCGAGCAGGCGGAAGGGAACGGAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGCATAATACATCCGTCAGCGCGAGTGTTTTACGCTGCg     >  2:127980/1‑251 (MQ=255)                                                                                                                                                                                                                                            cGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGCGTCt  >  1:149461/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        TATGAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGC‑CAC‑GCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCG‑TTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATA‑ACGCGATCAACGACATGGTGC‑GTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGCGTCT  >  NC_000913/974209‑974692

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