Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,095,232	G→C	9.8%	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	9.8%	158.2 / 8.8	51	intergenic (‑240/+77)	htpX/alaA	protease HtpX/glutamate‑pyruvate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (22/24);  new base C (3/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

AACTGTGTACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAG  >  NC_002947/2095083‑2095381                                                                                                                                                      |                                                                                                                                                       aaCTGTGTACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCcgcg                                                                                                                                                        >  1:22662/1‑150 (MQ=255)    tgtgTACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgg                                                                                                                                                     <  1:1634109/150‑1 (MQ=255)     gtgtACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcggg                                                                                                                                                    <  1:1241826/150‑1 (MQ=255)      tgtACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggc                                                                                                                                                   <  2:1646618/150‑1 (MQ=255)       gtACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggcg                                                                                                                                                  <  2:722937/150‑1 (MQ=255)            gcggAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                             >  2:1509020/1‑150 (MQ=255)            gcggAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                             <  2:584090/150‑1 (MQ=255)                  tATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCt                                                                                                                                       >  1:970497/1‑150 (MQ=255)                   aTGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTTAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTg                                                                                                                                      <  1:544434/150‑1 (MQ=255)                        agagGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGccc                                                                                                                                 >  2:770763/1‑150 (MQ=255)                        agagGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGccc                                                                                                                                 <  1:144599/150‑1 (MQ=255)                            gATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCccacc                                                                                                                             >  2:724581/1‑150 (MQ=255)                              tGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                         <  2:538342/136‑1 (MQ=255)                              tGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGAGCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                         >  1:538342/1‑136 (MQ=255)                              tGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCCAAGTGGCTTGGTAAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCcccccac                                                                                                                           >  2:964515/1‑150 (MQ=38)                                              ggATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTgcgc                                                                                                                                   <  2:1626021/126‑1 (MQ=255)                                              ggATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTgcgc                                                                                                                                   >  1:1626021/1‑126 (MQ=255)                                              ggATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAAc                                                                                                           >  2:1637918/1‑150 (MQ=255)                                               gATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCCCGTCGCAcc                                                                                                          >  2:935445/1‑150 (MQ=255)                                                 tAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACctc                                                                                                        >  2:1468828/1‑150 (MQ=255)                                                 tAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACctc                                                                                                        <  2:960031/150‑1 (MQ=255)                                                 tAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACctc                                                                                                        >  1:960031/1‑150 (MQ=255)                                                 tAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTATTGCTTGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACctc                                                                                                        <  1:919787/150‑1 (MQ=39)                                                   gagCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCtctc                                                                                                      <  1:1361193/150‑1 (MQ=255)                                                   gagCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCtctc                                                                                                      <  1:1465089/150‑1 (MQ=255)                                                   gagCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCtctc                                                                                                      >  2:1465089/1‑150 (MQ=255)                                                      cGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTa                                                                                                   <  1:1318403/150‑1 (MQ=255)                                                        ttGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑c                                                                                                <  1:296618/150‑1 (MQ=255)                                                                 gTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGACGCAAGTGGTGTGGTGCCTGTGAGATGGAGCGCCGCGCTGGCGGCGTTCGATCTGCGACCCACCACCAACCCCCGGTCCGAACTCTCGTAGACCCCCgggc                                                                                        >  1:707364/1‑148 (MQ=255)                                                                   aaGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGccc                                                                                     <  1:1468828/150‑1 (MQ=255)                                                                     gAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGccccc                                                                                   >  2:424297/1‑150 (MQ=255)                                                                        aCATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTc                                                                                <  1:724581/150‑1 (MQ=255)                                                                          aTTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGa                                                                              >  2:866605/1‑150 (MQ=255)                                                                           ttCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCCCCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGaa                                                                             <  2:1550952/150‑1 (MQ=255)                                                                             cGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCCCGTAG‑CCCCCGGCGCCACAGTCGAAcc                                                                           <  2:22662/150‑1 (MQ=255)                                                                                       gCTGCGTCAAGGCCGCTGGGGTTTTGCCGCAAGTGGTGTGGTGCCTGTGAGCTCGAGCGCCGCGCGGGCGGCGCTCGATCTGTGCCCCACCCCCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGaa                                                                 <  1:902355/150‑1 (MQ=39)                                                                                        cTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAAt                                                                <  1:1509020/150‑1 (MQ=255)                                                                                             tCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAAccc                                                           <  1:866605/150‑1 (MQ=255)                                                                                               aaGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGccc                                                                                     >  1:1086500/1‑122 (MQ=255)                                                                                               aaGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGccc                                                                                     <  2:1086500/122‑1 (MQ=255)                                                                                                        gggTGTTTGCCGCAAGTGGTGGGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAAtt                                                <  2:305867/150‑1 (MQ=255)                                                                                                         ggTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTg                                               <  1:1637918/150‑1 (MQ=255)                                                                                                                gCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATcc                                        <  1:1269404/150‑1 (MQ=255)                                                                                                                 ccGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCCCCCCCGCGTCGAACCTCTCGTAG‑CCCCCCGCGCCCCCGTCGAACCCTACTGCGAATAACCCCTCAGGAAATTGCCGATccc                                       >  1:917316/1‑149 (MQ=255)                                                                                                                 ccGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCCGCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCg                                       >  1:589436/1‑150 (MQ=255)                                                                                                                  cGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTa                                                                        >  2:65894/1‑116 (MQ=255)                                                                                                                  cGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTa                                                                        <  1:65894/116‑1 (MQ=255)                                                                                                                   gCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGc                                                                    <  1:926300/119‑1 (MQ=255)                                                                                                                   gCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGc                                                                    >  2:926300/1‑119 (MQ=255)                                                                                                                        tggtgtggtgCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAAtt                                <  1:1438894/150‑1 (MQ=255)                                                                                                                                  ccTGTGAGATCGAGCGCCGCGCGCGCGGCGCTCGTTCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGATCTCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCa                      <  2:333484/150‑1 (MQ=255)                                                                                                                                        agaTCGAGCGCCGCCCGCGCGGCGCTCGATCTGCGCCCCCCCACCACCCCCCCGTCGAAACTCTCGTAA‑CGCCCCGCCCCCCCGTCGAACCCTACTGCGAATAAACCTTCAGGAAAATTCCGGTCCCGCCAATTTCCCTCTTCAGGTCAt                >  1:742244/1‑150 (MQ=255)                                                                                                                                         gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCCGCGCCCCCGGCGAAACCCCCTGCGAATAAACCTTCAGGAAATTGCCGCTCCGGCCAATTGCCGCTTCCCgggtccc               >  1:1006640/1‑145 (MQ=255)                                                                                                                                          aTCGAGCGCCGCCCGCGCGGCGCTCGATCTGCGCCCCACCACCACACCCGCGTAGGACCTTGCGTAG‑CCCCCGGCGCCCCCGGCGAACCCTACTGCGAATAACCCTTAAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATcc              >  1:200884/1‑150 (MQ=255)                                                                                                                                                 gccgcGCGGGCGGCGCTCGATCTGCGCCCCCCCACCACCCCCGCGTCGAACCTCCCGTAG‑CCCCCCGCGCCCCCCCCGAAACCTACTGCGAATAACCCCTCACGAAAATGCCCATCCGGCCCATTGCCCCCTCCACGTCCTccccccccc       >  2:1449875/1‑147 (MQ=255)                                                                                                                                                 gccgcGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGGCGAACCCCTCGTAG‑CCCCCCGCGCCCCCGTCGAACCCTTCTTCGAAAAAACCTTCAGGAAATTGTCGGTCCGGCCAATTTCCCCTTCCAAGTCATCCACCcccg       >  2:1558276/1‑150 (MQ=255)                                                                                                                                                   cgcgcgGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCCAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCCTCCCGGTCATCCCCCcgccgc     >  2:808911/1‑150 (MQ=255)                                                                                                                                                     cgcgGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCaa   <  1:424297/150‑1 (MQ=255)                                                                                                                                                      gcgggcggCGCTCGATCTGCGCCCTACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAg  <  1:1310008/150‑1 (MQ=255)                                                                                                                                                      gcgggcggCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAg  <  2:707364/150‑1 (MQ=255)                                                                                                                                                      gcgggcggCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAg  <  1:1558276/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       AACTGTGTACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAG‑CCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAG  >  NC_002947/2095083‑2095381

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: ‑