Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	380,012	1	.	G	46.9%	101.8 / 143.6	64	intergenic (‑131/+57)	frmR/yaiO	regulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (17/17);  new base G (24/6);  total (41/23)
Fisher's exact test for biased strand distribution p-value = 1.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-02
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAC  >  NC_000913/379880‑380145                                                                                                                                       |                                                                                                                                            aTTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGGGGTCGGAAAACAAAGAGGAA‑A‑‑gggggg                                                                                                                                      <  2:1076639/139‑1 (MQ=255)   ttCGCACCTCATCATGTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑gggggggg                                                                                                                                    >  2:1008836/1‑139 (MQ=255)                           agataCTATAGGGGGGTATTATATATGTCAAGGCATACCCCCCTATAGTATATAAAACTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGgccg‑c                                                                                                         <  2:383296/137‑1 (MQ=255)                                   taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATGGGCAGGGAAGGCCC‑CCCCGGAtg                                                                                                 >  1:143972/1‑138 (MQ=255)                                        gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑gggct                                                                                           >  1:383296/1‑135 (MQ=255)                                        gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑gcgggc                                                                                          >  2:1213771/1‑139 (MQ=255)                                         ggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑gcgggc                                                                                          >  2:1412095/1‑139 (MQ=255)                                                     gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGgc                                                                                                             >  2:436054/1‑110 (MQ=255)                                                     gTCAATGCATACCCCCCTATAGTATATAACACGGGAGAATAAAATTTATCCGGTGAATGGGGTGGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGgc                                                                                                             <  1:436054/110‑1 (MQ=255)                                                          tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGa                                                                          <  2:1143927/138‑1 (MQ=255)                                                          tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGa                                                                          >  1:1143927/1‑138 (MQ=255)                                                          tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAAt                                                                        >  2:869382/1‑139 (MQ=255)                                                            cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCa                                                                            >  1:796576/1‑133 (MQ=255)                                                            cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCa                                                                            <  2:796576/133‑1 (MQ=255)                                                                    ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGAATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                                                               >  1:1159905/1‑139 (MQ=255)                                                                    ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CCAAACCGGCGGGAAAGCCGCCCCCGGAAG‑GCGGGCGGAAGAGGGAAACAGAATTTCCCggcc                                                               >  2:768631/1‑137 (MQ=255)                                                                       gaaGTATAAAAGACTCGAGAATAAAATTTATCCGGTGAAGTTGGTCGGAAAACAAAGAGGAACG‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CGCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                                                            <  1:1008836/137‑1 (MQ=255)                                                                           tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGTTCGGAAAGCAAAGAGGAA‑A‑‑GGGGGGGG‑G‑GTCTAATCCGCCGGCAAGGGCG‑CCCCGGATA‑GCGGGGGGCAGAAGGAAGCAGAATTTCGGGGTCAGCCggg                                                        >  2:572714/1‑139 (MQ=255)                                                                             tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AGGGGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc                                                       >  1:217318/1‑139 (MQ=255)                                                                             tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                                      >  2:281631/1‑139 (MQ=255)                                                                             tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAgg                                                                                                              >  2:803433/1‑84 (MQ=255)                                                                             tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg                                                             <  2:1111677/131‑1 (MQ=255)                                                                             tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg                                                             >  1:1111677/1‑131 (MQ=255)                                                                             tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                                                     >  2:377979/1‑139 (MQ=255)                                                                             tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                                                     >  2:803215/1‑139 (MQ=255)                                                                             tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                                                     >  2:316273/1‑139 (MQ=255)                                                                             tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAAACGGCAGGGAAGGCCG‑CCCCGGATA‑GGGGGGGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                                                     >  1:792267/1‑139 (MQ=255)                                                                             cataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAgg                                                                                                              <  1:803433/83‑1 (MQ=255)                                                                               taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                                                   >  2:291394/1‑139 (MQ=255)                                                                                aaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                                                   >  2:698463/1‑139 (MQ=255)                                                                                  cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AGGGGGGGGGG‑G‑G‑CTAATGGGAAGGGAAGGCC‑‑CCCCGGGTAGGCGGGGGGAAGAGGGAACCAGAATTTCCGGGTAAGACGGGcgccag                                                  >  1:673949/1‑134 (MQ=255)                                                                                  cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AGGGGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCaa                                                  >  2:23189/1‑139 (MQ=255)                                                                                  cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                 >  1:598721/1‑139 (MQ=255)                                                                                  cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                 >  1:1105706/1‑139 (MQ=255)                                                                                  cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                 >  2:23299/1‑139 (MQ=255)                                                                                  cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                                >  2:329778/1‑139 (MQ=255)                                                                                     tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAAttt                                                                      >  2:1423571/1‑115 (MQ=255)                                                                                     tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                    >  2:1497009/1‑133 (MQ=255)                                                                                     tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                              >  2:418221/1‑139 (MQ=255)                                                                                     tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                              >  1:1451982/1‑139 (MQ=255)                                                                                     tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑GTG‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                    <  1:1497009/133‑1 (MQ=255)                                                                                     tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGaa                                                                               >  1:310168/1‑105 (MQ=255)                                                                                     tGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAAttt                                                                      <  1:1423571/115‑1 (MQ=255)                                                                                     tGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGaa                                                                               <  2:310168/105‑1 (MQ=255)                                                                                      ggAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                            <  1:803215/139‑1 (MQ=255)                                                                                       gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                                                            >  1:1274481/1‑123 (MQ=255)                                                                                       gagaATAAAATTTATCCGGTGAAGGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                                                            <  2:1274481/123‑1 (MQ=255)                                                                                        aaaaTAAAATTTATCCGGGGAATGGGGAGGAAAAACAAAGAGGAA‑G‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                          <  2:792267/137‑1 (MQ=255)                                                                                            tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCg                                       >  1:103396/1‑139 (MQ=255)                                                                                            tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCg                                       >  2:1289308/1‑139 (MQ=255)                                                                                                aTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaa                                  <  1:329778/139‑1 (MQ=255)                                                                                                    aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatca                              >  1:573127/1‑139 (MQ=255)                                                                                                     tCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcat                             <  1:869382/139‑1 (MQ=255)                                                                                                         gTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGt                         <  1:291394/139‑1 (MQ=255)                                                                                                             aTGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                     <  2:573127/139‑1 (MQ=255)                                                                                                               gtgGTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCCGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg                    <  1:1412095/139‑1 (MQ=255)                                                                                                                tgGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtg                  <  1:316273/139‑1 (MQ=255)                                                                                                                  gTCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtgt                 >  2:1103902/1‑139 (MQ=255)                                                                                                                  gTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                                      <  2:149643/101‑1 (MQ=255)                                                                                                                  gTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                                      >  1:149643/1‑101 (MQ=255)                                                                                                                   tCGGAAAACAAAGAGGAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAAACAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGCTGGGGTGTc                >  2:1494320/1‑139 (MQ=255)                                                                                                                   tCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAgg                                                                                                              <  2:1443874/46‑1 (MQ=255)                                                                                                                   tCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAgg                                                                                                              >  1:1443874/1‑46 (MQ=255)                                                                                                                   acGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCg               <  1:377979/138‑1 (MQ=255)                                                                                                                    cGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGt              >  2:69603/1‑139 (MQ=255)                                                                                                                       aaaaCAAAGAGGAA‑A‑‑GGGGGGGGAG‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc            <  1:1289308/139‑1 (MQ=255)                                                                                                                                 ggAA‑AG‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCC‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAc  >  2:321308/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                            ATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑A‑‑GGGGGGGG‑G‑G‑CTAATCGGCAGGGAAGGCCG‑CCCCGGATA‑GCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAC  >  NC_000913/379880‑380145

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