Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,095,232	G→C	11.2%	intergenic (‑240/+77)	htpX ← / ← alaA	protease HtpX/glutamate‑pyruvate aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,095,232	0	G	C	11.2%	183.1 / 13.5	63	intergenic (‑240/+77)	htpX/alaA	protease HtpX/glutamate‑pyruvate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base C (3/4);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

TGTGTACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAA  >  NC_002947/2095086‑2095380                                                                                                                                                   |                                                                                                                                                         tgtgTACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTACGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgg                                                                                                                                                       <  2:1255965/150‑1 (MQ=255)    gtACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGTGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggcg                                                                                                                                                    <  2:699411/150‑1 (MQ=255)         gcggAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                               <  1:1754451/150‑1 (MQ=255)           ggAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                             <  1:1560916/150‑1 (MQ=255)            gAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                            <  1:2272058/150‑1 (MQ=255)                   gCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑gc‑‑gc                                                                                                                                  >  1:1658342/1‑150 (MQ=255)                      gagGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑Gc‑ccc                                                                                                                              >  2:545766/1‑150 (MQ=255)                      gagGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑Gc‑ccc                                                                                                                              >  2:2021283/1‑150 (MQ=255)                      gagGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑Gc‑ccc                                                                                                                              <  1:2021283/150‑1 (MQ=255)                       agGATGGCGTGGGGGCTGTGGGATAGACCGTTGCAGGATGTAAGACACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑Ccca                                                                                                                             <  1:918941/150‑1 (MQ=255)                                  ggggCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCAccccc                                                                                                                  >  1:2293818/1‑150 (MQ=255)                                       tgtgGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTc                                                                                                             <  2:2090088/150‑1 (MQ=255)                                       tgtgGGATAGAGCGTGACAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTATTTGCCGCAAGTGGTGTGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTc                                                                                                             <  2:46275/150‑1 (MQ=255)                                         tgGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGgcggcg                                                                                                                                                 >  2:307139/1‑116 (MQ=255)                                         tgGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGgcggcg                                                                                                                                                 <  1:307139/116‑1 (MQ=255)                                            gATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCCAAcc                                                                                                        >  2:231037/1‑150 (MQ=255)                                                gagCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                                >  2:2077780/1‑110 (MQ=255)                                                gagCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                                <  1:2077780/110‑1 (MQ=255)                                                 agCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCg                                                                                                   <  2:2283880/150‑1 (MQ=255)                                                 agCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCg                                                                                                   <  2:681194/150‑1 (MQ=255)                                                 agCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCg                                                                                                   <  2:1904911/150‑1 (MQ=255)                                                 agCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCCCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCg                                                                                                   <  2:1562472/150‑1 (MQ=255)                                                                 aGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCCCCACCACCCCCCCGTCGCACCTCCCGCAGCCCCCCCCGcccc                                                                                   >  2:507324/1‑150 (MQ=38)                                                                  gAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGccccc                                                                                  >  1:214614/1‑150 (MQ=255)                                                                   aaaCATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCg                                                                                 >  2:1181363/1‑150 (MQ=255)                                                                       aTTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGAGCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCCACACCACCCCCGCGCCCGACCTCTCGTAGCCCCCGGCGCCCCCCTCgg                                                                             >  1:2258303/1‑149 (MQ=38)                                                                          cGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAAcc                                                                          <  2:1417265/150‑1 (MQ=255)                                                                                    gCTGCGTCAAGGCCGCTGGGTGTTTGCCGCACGTGGGGTGGTGCCTGTGAGATCGCGCCGCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGTGCCCCCGTCGAACCCTACTGCGaa                                                                <  1:636216/150‑1 (MQ=255)                                                                                      tGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAAc                                                                                                         <  1:1492275/107‑1 (MQ=255)                                                                                      tGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAAc                                                                                                         >  2:1492275/1‑107 (MQ=255)                                                                                      tGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCtctc                                                                                                    >  1:1725399/1‑112 (MQ=255)                                                                                      tGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCtctc                                                                                                    <  2:1725399/112‑1 (MQ=255)                                                                                           cAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCt                                                         <  2:836873/150‑1 (MQ=255)                                                                                              ggCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGCGCGTCACTCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                      <  1:149362/150‑1 (MQ=255)                                                                                              ggCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                      <  2:2293818/150‑1 (MQ=255)                                                                                              ggCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                      <  2:2302829/150‑1 (MQ=255)                                                                                              ggCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCa                                                      <  1:1061409/150‑1 (MQ=255)                                                                                                     gggTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAAcc                                                                          >  1:2158023/1‑123 (MQ=255)                                                                                                     gggTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAAcc                                                                          <  2:2158023/123‑1 (MQ=255)                                                                                                      ggTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCCAACCTCTCGTAACCCCCGGCGCCCCCGTCGAACCCCACTGCGAATAAACCTTCAGGAAATTg                                              >  2:1850174/1‑150 (MQ=255)                                                                                                       gtgtTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGc                                             <  2:1658342/150‑1 (MQ=255)                                                                                                          tttGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGa                                          <  2:551747/150‑1 (MQ=255)                                                                                                              ccGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCCCCTCGAACCCCTCGGAGCCCCCCGCCCCCCCCCCCAACCCCACTGCCAAAAAACCCTCAGGAAAATGCCGATCCg                                      >  2:2038344/1‑150 (MQ=25)                                                                                                               cGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCCTCGAACCTCTCGTAGCCCCCGGCGCCCCCCCCCAACCCCACTGCGAATAAACCCTCAGGAAAATGCCCATcccg                                     >  2:1657831/1‑148 (MQ=38)                                                                                                                     tggtgtggtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCTTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAAtt                               <  2:214614/150‑1 (MQ=255)                                                                                                                       gtgtggtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCCCCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGCCGCAACCCACTGCGAATAACCCCTCCCGAAAATGCCGATCCGGCCCAtttc                             >  2:445736/1‑148 (MQ=37)                                                                                                                        tgtggtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTa                                                                                                 >  1:10970/1‑81 (MQ=255)                                                                                                                        tgtggtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTa                                                                                                 <  2:10970/81‑1 (MQ=255)                                                                                                                            gtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTc                                                       <  1:116044/119‑1 (MQ=255)                                                                                                                            gtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTc                                                       >  2:116044/1‑119 (MQ=255)                                                                                                                                  gtgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGt                  <  1:1181363/150‑1 (MQ=255)                                                                                                                                  cgcagaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGC‑ACCA‑CACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAg                    <  2:2204934/147‑1 (MQ=25)                                                                                                                                   tgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACt                                                                     <  1:1032856/98‑1 (MQ=255)                                                                                                                                   tgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACt                                                                     >  2:1032856/1‑98 (MQ=255)                                                                                                                                    gagaTCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCCCCACCCCCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGGACCCTACTTCGGAAAACCCTTCAGGAAAATGCCGCTCCGGCCAATTTCCGCTTCCCGGTcc                >  1:1681506/1‑149 (MQ=255)                                                                                                                                       aTCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCCCCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTTAGGAAATTGCCGATCCGGCCCATTTCCGCTTCCAGGTTATcc             >  1:987510/1‑150 (MQ=255)                                                                                                                                       aTCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAACCCCCGGCGCCCCCGTCCAACCCTACTGCGAATAACCCTTCAAGAAATTGCCGATCCGGCCACTTGCCCCTTCCAGGtcctcc             >  2:160401/1‑150 (MQ=255)                                                                                                                                        tCGAGCGCCGCCCGCGCGGCGCTCGATCT‑‑C‑‑ACAGGCACCACACCACTTGCGTCAAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAACTTGCCGATCCGGCCAATGGCCGCTTCCAGGTCATccc            >  2:1411633/1‑149 (MQ=18)                                                                                                                                          gAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCCCCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTc                                                                               >  2:1508319/1‑81 (MQ=255)                                                                                                                                          gAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTc                                                                               <  1:1508319/81‑1 (MQ=255)                                                                                                                                          gAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCCCGCCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGCACCCTACTTCGAATAACCCTTCCGGGAAACTCCGATCCGGCCCATTGCCGCTTCCAGGTCATCcccc          >  1:1833757/1‑150 (MQ=255)                                                                                                                                          gAGCGCCGCCCGCGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGGACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACCGCGAATAACCCTTCAGGACATTGACCATGCGGCCAATTGCCGCTTCCAGGTCATCCAcc          >  2:192053/1‑150 (MQ=255)                                                                                                                                           aGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCtctc                                                                                                    <  1:1552027/59‑1 (MQ=255)                                                                                                                                           aGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCtctc                                                                                                    >  2:1552027/1‑59 (MQ=255)                                                                                                                                               ccgcGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCgg     >  1:48746/1‑150 (MQ=255)                                                                                                                                                cgcgcgGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGGCGAACCCCTCGTAGCCCCCCGCGCCCCCGGCGAACCCTACCGCGAATAACCCCTCAGGAAAATGCCGATCCGGCCCAATTCCGCTTCCCGGGCATCCCCCCCCGGc    >  2:2384076/1‑150 (MQ=255)                                                                                                                                                  cgcgGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCaa  <  1:545766/150‑1 (MQ=255)                                                                                                                                                  cgcgGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCaa  <  1:2314012/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                         TGTGTACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCT‑GC‑‑GC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAA  >  NC_002947/2095086‑2095380

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