Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,673,177	A→C	6.6%	D164A (GAT→GCT)	PP_2340 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,673,177	0	A	C	6.6%	184.3 / 4.2	61	D164A (GAT→GCT)	PP_2340	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base C (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

TCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAT  >  NC_002947/2673030‑2673320                                                                                                                                                      |                                                                                                                                                tCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTAACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGa                                                                                                                                                 <  1:49447/149‑1 (MQ=255)  ccGTGCCCTGGCGGCCTAGGGCATGGTGTTGCTGGCGGAGGTCCGGGACTCGGAACTGCCCCAATTGCGCGTTGCGATCG‑ATCGCGAGTTCTTCCAGCATTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCAGCGCGAt                                                                                                                                                <  6:197130/148‑1 (MQ=255)  ccGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATg                                                                                                                                               >  5:300232/1‑149 (MQ=255)  ccGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATg                                                                                                                                               >  7:205859/1‑149 (MQ=255)  ccGTGACCTGGAGGCCTATGGCATGGTGTAGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCAAGCTTCCTGCGCGATg                                                                                                                                               >  5:305889/1‑149 (MQ=255)            gAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCa                                                                                                                                      <  1:278234/148‑1 (MQ=255)                    tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAgg                                                                                                                             >  8:108045/1‑149 (MQ=255)                    tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCa                                                                                                                               <  8:19848/147‑1 (MQ=255)                      gCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCCCGCAGGGc                                                                                                                           >  7:136734/1‑149 (MQ=255)                           gtgtTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTg                                                                                                                      <  2:81646/149‑1 (MQ=255)                           gtgtCGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGt                                                                                                                       <  7:166648/148‑1 (MQ=255)                                gctgcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCg                                                                                                                 >  8:113614/1‑149 (MQ=255)                                   gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                              <  8:205859/149‑1 (MQ=255)                                    ccGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                              <  6:284291/148‑1 (MQ=255)                                       gAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAAt                                                                                                          <  6:300232/149‑1 (MQ=255)                                         gTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCa                                                                                                        >  4:67645/1‑149 (MQ=255)                                             cGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCg                                                                                                    <  3:67645/149‑1 (MQ=255)                                                gACTGGGAACTGCCCCAATTGCGCGATGAGATCG‑ATCGCGAGTTATTCCTTCTTTAAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTACTGCGCGATG‑CATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCgg                                                                                                >  1:23512/1‑149 (MQ=255)                                                 aCTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAgg                                                                                                                             >  3:154832/1‑120 (MQ=255)                                                 aCTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAgg                                                                                                                             <  4:154832/120‑1 (MQ=255)                                                     ggAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTa                                                                                            <  3:335425/149‑1 (MQ=255)                                                       aaCTGCCCCAATTGCGCGATGCGATCG‑AGCGCGAGGTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAtt                                                                                          <  3:244660/149‑1 (MQ=255)                                                        aCTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTTTATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTc                                                                                         <  3:293412/149‑1 (MQ=255)                                                            ccccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCg                                                                                                 <  2:69601/137‑1 (MQ=255)                                                            ccccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCg                                                                                                 >  1:69601/1‑137 (MQ=255)                                                             cccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAAc                                                                                    >  8:319293/1‑149 (MQ=255)                                                              ccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAAcc                                                                                   >  3:23867/1‑149 (MQ=255)                                                              ccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAAcc                                                                                   >  8:22324/1‑149 (MQ=255)                                                               cAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCg                                                                                  >  6:98406/1‑149 (MQ=255)                                                                  ttGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                               <  1:309264/149‑1 (MQ=255)                                                                    tcgcgAT‑CGATCGCATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGt                                                                             >  7:200731/2‑149 (MQ=255)                                                                    gcgcGATGCGATCA‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCTATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTgg                                                                              <  7:124947/148‑1 (MQ=255)                                                                     cgcgATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGt                                                                             <  7:319293/148‑1 (MQ=255)                                                                      gcgatgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCCGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTg                                                                            >  4:219147/1‑148 (MQ=255)                                                                        catgctatCG‑ATCGCGAGTTCTTACTGCTTTACCACGCGCACT‑ACCAGCGCTTTATCGCCCCCAGCTTGCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAg                                                                         <  8:136734/148‑1 (MQ=255)                                                                         atgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTCCGCGATTCCATTGGCAGCCCGCAGGGCCAGTGACCCGGGCGAATCAGTAGCCGTTGTATTCGAAACCCCTGGTGCAgt                                                                        >  4:194086/1‑148 (MQ=255)                                                                               tcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGcccat                                                                                                                                          <  7:125905/77‑4 (MQ=255)                                                                               tcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGcccat                                                                                                                                          >  8:125905/1‑74 (MQ=255)                                                                                    tcgCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                              >  5:298048/1‑149 (MQ=255)                                                                                      gcgAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGc                                                            >  6:129393/1‑149 (MQ=255)                                                                                          gTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAAc                                                        <  5:301516/149‑1 (MQ=255)                                                                                               tcCTGCTTTCAGCGGCGCACT‑ACCCGCGCTATAGCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                    <  4:23867/148‑1 (MQ=255)                                                                                                         agagGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCg                                          <  7:22324/148‑1 (MQ=255)                                                                                                             gcgcACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATG‑CATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCa                                    <  2:23512/149‑1 (MQ=255)                                                                                                                 aCT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCg                                 <  7:76532/149‑1 (MQ=255)                                                                                                                        aGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGa                           >  8:338688/1‑149 (MQ=255)                                                                                                                         gcgcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                  <  6:177602/109‑1 (MQ=255)                                                                                                                         gcgcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                  >  5:177602/1‑109 (MQ=255)                                                                                                                           gcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCt                        <  8:200731/149‑1 (MQ=255)                                                                                                                           gcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCt                        >  2:19066/1‑149 (MQ=255)                                                                                                                             tataTCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTgg                      >  5:34402/1‑149 (MQ=255)                                                                                                                                atCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                   >  4:34389/1‑149 (MQ=255)                                                                                                                                 tCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTccc                  <  1:19066/149‑1 (MQ=255)                                                                                                                                 tCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                   >  5:340398/1‑148 (MQ=255)                                                                                                                                    cccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGAc                                                 <  1:171492/115‑1 (MQ=255)                                                                                                                                    cccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGAc                                                 >  2:171492/1‑115 (MQ=255)                                                                                                                                      cccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGAt              >  3:341492/1‑148 (MQ=255)                                                                                                                                              ccTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGc     >  6:59034/1‑149 (MQ=255)                                                                                                                                              ccTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGc     >  2:132678/1‑149 (MQ=255)                                                                                                                                              ccTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGc     >  8:239141/1‑149 (MQ=255)                                                                                                                                                 gcgcgcTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt  >  4:272914/1‑149 (MQ=255)                                                                                                                                                 gcgcGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt  >  3:261361/1‑149 (MQ=255)                                                                                                                                                 gcgcGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt  >  5:281894/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                TCCGTGACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACT‑ACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAT  >  NC_002947/2673030‑2673320

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