Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,740,308	A→T	6.5%	I279F (ATC→TTC) ‡	gltA →	citrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,740,308	0	A	T	6.5%	230.9 / 4.9	77	I279F (ATC→TTC) ‡	gltA	citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base T (2/3);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

CACGAGCAGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTG  >  NC_002947/4740161‑4740454                                                                                                                                                      |                                                                                                                                                    cACGAGCAGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAAt                                                                                                                                                    >  1:290653/1‑149 (MQ=255)    gAGCAGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCgg                                                                                                                                                 >  6:15440/1‑149 (MQ=255)    gAGCAGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCgg                                                                                                                                                 >  8:44408/1‑149 (MQ=255)     aGCAGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGc                                                                                                                                                <  1:317426/149‑1 (MQ=255)      gCAGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTCCGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGc                                                                                                                                                <  1:342623/149‑1 (MQ=255)      gCAGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCg                                                                                                                                               <  7:112164/149‑1 (MQ=255)       cAGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGa                                                                                                                                              <  2:98535/149‑1 (MQ=255)        aGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGAt                                                                                                                                             <  7:44408/149‑1 (MQ=255)        aGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGAt                                                                                                                                             >  3:390568/1‑149 (MQ=255)         gAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATg                                                                                                                                            >  2:71911/1‑149 (MQ=255)          aaCGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGt                                                                                                                                           >  8:391829/1‑149 (MQ=255)               ctccacctccacCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGa                                                                                                                                      <  7:255218/149‑1 (MQ=255)                   acctccaccGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACAt                                                                                                                                 <  3:67454/149‑1 (MQ=255)                   acctccaccGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACAt                                                                                                                                 <  4:100269/149‑1 (MQ=255)                    cctccaccGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATc                                                                                                                                >  1:265743/1‑149 (MQ=255)                        caccGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGAc                                                                                                                             >  5:355087/1‑148 (MQ=255)                        caccGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACa                                                                                                                            <  3:54510/149‑1 (MQ=255)                                gTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ATATCGACAAATTCATc                                                                                                                    <  4:282674/149‑1 (MQ=255)                                 gcttgCAGGCTCGTCGGGTGCCACCCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑TGGGCCGGT‑GCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATCAAATCGGCTATGTCTCGA‑ACATCGACAAATTCATCg                                                                                                                   <  8:296866/145‑1 (MQ=255)                                    ggCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCgcc                                                                                                                 >  6:247082/1‑148 (MQ=255)                                     gCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGT‑TCGAGACATCGCCGATTTCATcca                                                                                                                  <  4:348323/146‑3 (MQ=37)                                     gCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGc                                                                                                                                                >  4:337219/1‑117 (MQ=255)                                     gCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGc                                                                                                                                                <  3:337219/117‑1 (MQ=255)                                        ggCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCgccaaggc                                                                                                            <  2:133240/149‑1 (MQ=255)                                          cTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGT‑TCGAGACATCGCCGATTTCATcca                                                                                                                  >  7:345695/1‑139 (MQ=37)                                          cTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGT‑TCGAGACATCGCCGATTTCATcca                                                                                                                  <  8:345695/141‑3 (MQ=37)                                           tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGT‑TCGAGACATCGCCGATTTCATccagcatg                                                                                                             <  2:148152/145‑8 (MQ=25)                                           tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGT‑TCGAGACATCGCCGATTTCATccagcatg                                                                                                             >  1:148152/1‑138 (MQ=25)                                           tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAAt                                                                                                                         >  4:418845/1‑133 (MQ=255)                                           tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAAt                                                                                                                         <  3:418845/133‑1 (MQ=255)                                           tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAAGTCATCgccaaggcc                                                                                                           >  8:295858/1‑147 (MQ=255)                                             gtcgTGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCgccaaggccaag                                                                                                        >  3:231877/1‑148 (MQ=255)                                              tcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGa                                                                                                      >  2:1205/1‑149 (MQ=255)                                                     ccAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATc                                                                                                                    >  7:282576/1‑128 (MQ=255)                                                     ccAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATc                                                                                                                    <  8:282576/128‑1 (MQ=255)                                                        aCCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGAc                                                                                            <  2:265743/149‑1 (MQ=255)                                                                gCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTc                                                                                      >  1:229661/1‑147 (MQ=255)                                                                gCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCaa                                                                                    <  4:390568/149‑1 (MQ=255)                                                                    gTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACaaa                                                                                                                          >  5:132170/1‑107 (MQ=255)                                                                    gTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACaaa                                                                                                                          <  6:132170/107‑1 (MQ=255)                                                                    gTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCgc                                                                                                                  <  6:427988/115‑1 (MQ=255)                                                                    gTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCgc                                                                                                                  >  5:427988/1‑115 (MQ=255)                                                                     tATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGa                                                                               <  7:371307/149‑1 (MQ=255)                                                                     tATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGa                                                                               <  7:391829/149‑1 (MQ=255)                                                                       tCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCgccaaggccaa                                                                                                         <  8:6564/121‑1 (MQ=255)                                                                       tCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCgccaaggccaa                                                                                                         >  7:6564/1‑121 (MQ=255)                                                                       tCGCAGACGTTATCGAAGAACTGT‑GGGTCCAGG‑CCAAATTCAGTGCGAACGAAGACGTACTGACCATGCTGGATGTAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGACAAGGCAAAGGACAAGAACGACCCCTTCTAGCTGAtc                                                                             >  4:70823/1‑148 (MQ=255)                                                                        cgccgcCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGACAAGGCCAAGGACAAGAACGACCCGTTCAAGCGGATgg                                                                            >  8:290221/1‑149 (MQ=255)                                                                           cgccgGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCCAGCTGATGGGCt                                                                         >  6:202496/1‑149 (MQ=255)                                                                            gccgGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCtt                                                                        <  1:98589/149‑1 (MQ=255)                                                                              cgGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCg                                                                      <  2:290653/149‑1 (MQ=255)                                                                                gTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGt                                                                    >  1:317769/1‑149 (MQ=255)                                                                                 tATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGt                                                                    >  1:420081/1‑148 (MQ=255)                                                                                   tCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCa                                                                  <  2:342488/148‑1 (MQ=255)                                                                                      ccgcACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCGAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACcgc                                                              >  7:141271/1‑149 (MQ=255)                                                                                        gcACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCgt                                                            <  2:229661/149‑1 (MQ=255)                                                                                            gccgtgggccgggccCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCgtg                                                           >  6:431244/6‑148 (MQ=255)                                                                                                ggggCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTGAAGCTGATGGGCTTCGGTCACCGCGTGTAc                                                        >  1:116012/1‑146 (MQ=255)                                                                                                ggggCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAg                                                     <  5:15440/149‑1 (MQ=255)                                                                                                  ggcccgg‑cccACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCCAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGaa                                                   >  7:102704/1‑149 (MQ=255)                                                                                                   gcccgg‑cccACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAAc                                                  >  6:431162/1‑149 (MQ=255)                                                                                                            cACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGAccc                                          >  8:197359/1‑149 (MQ=255)                                                                                                                      cgaacgaaGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCAcc                                >  6:36195/1‑149 (MQ=255)                                                                                                                          cgaaGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGtga                            <  1:71911/149‑1 (MQ=255)                                                                                                                           gaaGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGtgat                           >  6:376579/1‑149 (MQ=255)                                                                                                                            aaGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGtgatg                          <  5:247082/149‑1 (MQ=255)                                                                                                                               ccGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAg                       <  7:197359/149‑1 (MQ=255)                                                                                                                                   aCTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGa                   >  8:364174/1‑149 (MQ=255)                                                                                                                                    cTGACCATGCTCGTTGAAATCCGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTTTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGAc                  <  4:238617/149‑1 (MQ=255)                                                                                                                                       aCCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCt                <  5:431162/148‑1 (MQ=255)                                                                                                                                        ccATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGc              <  5:431244/149‑1 (MQ=255)                                                                                                                                         cATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGcg             >  3:33868/1‑149 (MQ=255)                                                                                                                                          aTGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGcga            <  7:364174/149‑1 (MQ=255)                                                                                                                                             cTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGcgacg          <  2:317769/148‑1 (MQ=255)                                                                                                                                              tGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGaa        <  2:116012/149‑1 (MQ=255)                                                                                                                                                gATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCCACGaa        <  7:392771/147‑1 (MQ=255)                                                                                                                                                   gAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCt   >  3:213457/1‑149 (MQ=255)                                                                                                                                                    aaaTCGGCGTTGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTg  >  6:411107/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CACGAGCAGAACGCCTCCACCTCCACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGT‑GGGGCCCGG‑CCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGA‑ACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTG  >  NC_002947/4740161‑4740454

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