Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,799,223	A→G	9.0%	intergenic (‑194/+22)	PP_2452 ← / ← ansB	hypothetical protein/glutaminase‑asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,799,223	0	A	G	9.0%	217.1 / 9.2	78	intergenic (‑194/+22)	PP_2452/ansB	hypothetical protein/glutaminase‑asparaginase
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base G (4/3);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

TAAGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTG  >  NC_002947/2799075‑2799365                                                                                                                                                     |                                                                                                                                               taagtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCAc                                                                                                                                               >  2:937548/1‑150 (MQ=255) taagtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCAc                                                                                                                                               >  2:1278936/1‑150 (MQ=255) taagtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCAc                                                                                                                                               >  2:714692/1‑150 (MQ=255)      aaaTTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGcc                                                                                                                                          <  2:1208661/150‑1 (MQ=255)           cGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAg                                                                                                                                     <  2:932859/150‑1 (MQ=255)               aCAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAggggg                                                                                                                                 <  2:796759/150‑1 (MQ=255)                cAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGa                                                                                                                                <  1:937548/150‑1 (MQ=255)                 aGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGaa                                                                                                                               >  2:948558/1‑150 (MQ=255)                  gCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAg                                                                                                                              <  1:492981/150‑1 (MQ=255)                       tggggAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAATCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTc                                                                                                                         >  2:217662/4‑150 (MQ=255)                         aGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCTGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCTGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAg                                                                                                                       >  2:1029024/1‑150 (MQ=255)                          ggAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAgg                                                                                                                                    <  2:652552/136‑1 (MQ=255)                          ggAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAgg                                                                                                                                    >  1:652552/1‑136 (MQ=255)                           gagaAATCCAGCGTTGGCCTGTAAACCATTAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTa                                                                                                                     <  1:602159/150‑1 (MQ=255)                             gaAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACt                                                                                                                   <  1:984339/150‑1 (MQ=255)                              aaaTCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTc                                                                                                                  <  1:714692/150‑1 (MQ=255)                               aaTCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTcc                                                                                                                 <  1:803267/150‑1 (MQ=255)                                    aGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCagaa                                                                                                            <  2:127141/150‑1 (MQ=255)                                          ggCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCg                                                                                                      >  1:819730/1‑150 (MQ=255)                                              tGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgg                                                                                                                                    >  1:765162/1‑116 (MQ=255)                                              tGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgg                                                                                                                                    <  2:765162/116‑1 (MQ=255)                                               gTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCa                                                                                                 >  2:1217282/1‑150 (MQ=255)                                                tAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAg                                                                                                >  2:164615/1‑150 (MQ=255)                                                        gAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCa                                                                                                               >  1:330723/1‑127 (MQ=255)                                                        gAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCa                                                                                                               <  2:330723/127‑1 (MQ=255)                                                           gCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGt                                                                                     <  2:547679/150‑1 (MQ=255)                                                                  aGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATcc                                                                                                       <  2:1228442/125‑1 (MQ=255)                                                                  aGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGCGGAAGCACTCAGTACTCCCAGAAGATcc                                                                                                       >  1:1228442/1‑125 (MQ=255)                                                                   gAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                         <  2:1010366/138‑1 (MQ=255)                                                                   gAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                         >  1:1010366/1‑138 (MQ=255)                                                                   gAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTc                                                                             >  1:677835/1‑150 (MQ=255)                                                                    aaaGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAgg                                                                                                                                    <  2:1195340/94‑1 (MQ=38)                                                                    aaaGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAgg                                                                                                                                    >  1:1195340/1‑94 (MQ=38)                                                                    aaaGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTc                                                                             >  2:508921/1‑149 (MQ=255)                                                                    aaaGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTc                                                                             <  1:508921/149‑1 (MQ=255)                                                                    aaaGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCt                                                                            >  1:1077105/1‑150 (MQ=255)                                                                      aGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCTGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTg                                                                          <  1:1029024/150‑1 (MQ=255)                                                                      aGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGGAGGGGGGGGGGGGGGGGGGGGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTg                                                                          <  2:105064/150‑1 (MQ=255)                                                                        gTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGt                                                                        >  2:20534/1‑150 (MQ=255)                                                                           aaCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATcc                                                                                                       >  2:490172/1‑116 (MQ=255)                                                                           aaCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATcc                                                                                                       <  1:490172/116‑1 (MQ=255)                                                                             cGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCg                                                                   >  2:135073/1‑150 (MQ=255)                                                                              gATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGc                                                                  >  2:697895/1‑150 (MQ=255)                                                                                    gCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTTGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAACAAACGCAGAAGATCCGCTGCAGTGCTTTTCTGTCCTGGGTCTTGGTTATCGCTACCAt                                                            >  1:978911/1‑150 (MQ=255)                                                                                    gCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTg                                                                                                                                       >  2:1185312/1‑75 (MQ=255)                                                                                    gCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTg                                                                                                                                       <  1:1185312/75‑1 (MQ=255)                                                                                          ttGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAg                                                      >  1:567255/1‑150 (MQ=255)                                                                                                  tctGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCAt                                                                     <  1:370301/127‑1 (MQ=255)                                                                                                  tctGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCAt                                                                     >  2:370301/1‑127 (MQ=255)                                                                                                   ctGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCgcgc                                             >  1:1166919/1‑150 (MQ=255)                                                                                                       gtCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGGACTCAGTa                                                                                                                     >  2:921735/1‑74 (MQ=255)                                                                                                       gtCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTa                                                                                                                     <  1:921735/74‑1 (MQ=255)                                                                                                        tCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACt                                                                                                                          <  2:117108/68‑1 (MQ=255)                                                                                                        tCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACt                                                                                                                          >  1:117108/1‑68 (MQ=255)                                                                                                          tACGCGGTTGTGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAACCACTCAGTACGCCCAGAAGATCCGCTGCAGCTGTTTGCTGTCCTGGGACTTGGTCATCGGTACCATCGCCAGGATCCGCGCTTTTTgt                                      >  2:687928/1‑149 (MQ=255)                                                                                                            cgcgGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCgg                                    <  1:948558/150‑1 (MQ=255)                                                                                                            cgcgGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCgg                                    <  2:677835/150‑1 (MQ=255)                                                                                                                gTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTc                                >  2:1033392/1‑150 (MQ=255)                                                                                                                 ttGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCa                               <  2:1077105/150‑1 (MQ=255)                                                                                                                  tGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAg                              <  1:135073/150‑1 (MQ=255)                                                                                                                   gAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCt                                                                 <  1:1023652/114‑1 (MQ=255)                                                                                                                   gAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCt                                                                 >  2:1023652/1‑114 (MQ=255)                                                                                                                   gAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAgg                             >  1:559427/1‑150 (MQ=255)                                                                                                                     ggAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTc                           >  1:645444/1‑150 (MQ=255)                                                                                                                       aaTTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGt                         <  1:697895/150‑1 (MQ=255)                                                                                                                       aaTTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGt                         >  2:985407/1‑150 (MQ=255)                                                                                                                         ttGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                       <  2:819730/150‑1 (MQ=255)                                                                                                                            gCAACTAGGCGTTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGcc                    <  1:217662/150‑1 (MQ=255)                                                                                                                                 ttcccccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAc               >  2:284389/5‑150 (MQ=255)                                                                                                                                    gCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCa            >  2:239826/1‑150 (MQ=255)                                                                                                                                      ccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGt          >  2:289841/1‑150 (MQ=255)                                                                                                                                      ccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                                >  2:1208311/1‑80 (MQ=255)                                                                                                                                      ccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                                <  1:1208311/80‑1 (MQ=255)                                                                                                                                      ccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCa            <  2:498942/148‑1 (MQ=255)                                                                                                                                      ccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCa            >  1:498942/1‑148 (MQ=255)                                                                                                                                      ccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATACGCGATTGTTGCGGGTTCAGGTCGTGAGCCACGACCCAGt          >  2:621897/1‑150 (MQ=255)                                                                                                                                              ccGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTg  <  2:567255/150‑1 (MQ=255)                                                                                                                                              ccGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTg  <  1:1278936/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               TAAGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTG  >  NC_002947/2799075‑2799365

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