Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	855,668	G→A	10.4%	P173P (CCG→CCA)	murI →	glutamate racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	855,668	0	G	A	10.4%	180.4 / 8.9	67	P173P (CCG→CCA)	murI	glutamate racemase
Reads supporting (aligned to +/- strand):  ref base G (26/34);  new base A (2/5);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

GTGCCAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTG  >  NC_002947/855523‑855816                                                                                                                                                  |                                                                                                                                                     gTGCCAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTa                                                                                                                                                   <  7:64034/149‑1 (MQ=255)    ccAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACATGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTg                                                                                                                                                >  5:322268/1‑149 (MQ=255)           cGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCAGTAGCGGCTGCa                                                                                                                                          <  2:96184/148‑1 (MQ=255)            gCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAgg                                                                                                                                        >  6:15074/1‑149 (MQ=255)            gCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAg                                                                                                                                         <  3:186309/148‑1 (MQ=255)            gCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCAGTAGCGGCTGCac                                                                                                                                         <  3:144731/148‑2 (MQ=255)               gCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCAGTAGCGGCTGCacataa                                                                                                                                     <  2:93097/149‑6 (MQ=255)                 cTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAgg                                                                                                                                        >  6:48021/1‑144 (MQ=255)                 cTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAgg                                                                                                                                        <  5:48021/144‑1 (MQ=255)                     cTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATa                                                                                                                                <  6:101506/148‑1 (MQ=255)                      tCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCTGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACg                                                                                                                              >  3:87204/1‑149 (MQ=255)                      tCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCAGTAGCGGCTGCacataacc                                                                                                                                   >  8:234970/1‑137 (MQ=255)                      tCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCAGTAGCGGCTGCacataacc                                                                                                                                   <  7:234970/144‑8 (MQ=255)                       cGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGc                                                                                                                             <  7:120937/149‑1 (MQ=255)                        gATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGc                                                                                                                             <  5:15074/148‑1 (MQ=255)                            gCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGAt                                                                                                                         <  6:273694/148‑1 (MQ=255)                                        aCGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCAc                                                                                                             <  6:123466/148‑1 (MQ=255)                                             cAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTa                                                                                                       >  5:92290/1‑149 (MQ=255)                                                gTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTAccc                                                                                                    <  1:34387/149‑1 (MQ=255)                                                 ttattaCCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCAt                                                                                                         >  2:164549/1‑143 (MQ=255)                                                 ttattaCCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCAt                                                                                                         <  1:164549/143‑1 (MQ=255)                                                 ttattaCCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCt                                                                                                   >  5:130364/1‑149 (MQ=255)                                                   tttaggcAGCCTTGCCCGTGGCTGGTCGAGCTGATCGATTCCTGTGACGTGTGAAGCCCGGGGGTGCGCCAGCTGCTGCAAGGGTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCtt                                                                                                  <  5:174184/142‑1 (MQ=255)                                                       cccAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCtt                                                                                             >  2:253500/1‑149 (MQ=255)                                                       cccAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCtt                                                                                             >  2:134600/1‑149 (MQ=255)                                                          aGCCTTGCCCGGGGCTGGTCGAGCTGATCGTAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCATGGGTTTGTGCAGCCGCTACTGGCTGGAGGTTGCGATACGCTGATCCTCGGCTGCACCCTTTACCCTTTCCTTCGc                                                                                          <  7:81515/149‑1 (MQ=255)                                                                    ggggCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTgg                                                                                                                                               <  5:59970/86‑1 (MQ=255)                                                                    ggggCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTgg                                                                                                                                               >  6:59970/1‑86 (MQ=255)                                                                    ggggCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTgg                                                                                <  1:157484/149‑1 (MQ=255)                                                                     gtgCTGGTCGAGCTGAGCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGATGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTAACCTTTCCTTCGCCCATTGCTGGc                                                                               <  6:322268/147‑1 (MQ=255)                                                                       gCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCg                                                                             <  6:92290/149‑1 (MQ=255)                                                                        cTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCgg                                                                            >  3:241866/1‑149 (MQ=255)                                                                         tGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCgg                                                                            >  8:184624/1‑148 (MQ=255)                                                                              gAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATgg                                                                       <  2:89607/148‑1 (MQ=255)                                                                                 cTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTccc                                                                   >  5:183608/1‑149 (MQ=255)                                                                                    aTCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGa                                                                >  4:204020/1‑149 (MQ=255)                                                                                    aTCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGa                                                                >  8:109715/1‑149 (MQ=255)                                                                                    aTCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGa                                                                <  5:41079/149‑1 (MQ=255)                                                                                       gAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCt                                                                                                                 >  4:260779/1‑97 (MQ=255)                                                                                       gAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCt                                                                                                                 <  3:260779/97‑1 (MQ=255)                                                                                       gAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATg                                                              <  6:183608/148‑1 (MQ=255)                                                                                           ccGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTg                                                         <  1:253500/149‑1 (MQ=255)                                                                                           ccGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTg                                                         <  4:241866/149‑1 (MQ=255)                                                                                            cGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTgg                                                        <  7:184624/149‑1 (MQ=255)                                                                                              gTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGc                                                       >  6:162853/1‑148 (MQ=255)                                                                                                  ccTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCcatca                                                  >  6:135197/1‑149 (MQ=255)                                                                                                   cTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCgg                                                                                                                   >  7:303610/1‑83 (MQ=255)                                                                                                   cTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCgg                                                                                                                   <  8:303610/83‑1 (MQ=255)                                                                                                    tGGTAAGCCCGGCGGTGCGTCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGACCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCatcatc                                                >  2:6218/1‑149 (MQ=255)                                                                                                    tGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGa                                                                                                                          >  8:261942/1‑75 (MQ=255)                                                                                                    tGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGa                                                                                                                          <  7:261942/75‑1 (MQ=255)                                                                                                      gTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCatcatc                                                <  6:130364/147‑1 (MQ=255)                                                                                                        aaGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGaca                                            >  5:262216/1‑149 (MQ=255)                                                                                                                cggTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGcgcg                                    <  4:87204/149‑1 (MQ=255)                                                                                                                        cAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGgc                            <  1:134600/149‑1 (MQ=255)                                                                                                                          gctgctGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTgg                                                        <  2:210722/119‑1 (MQ=255)                                                                                                                          gctgctGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTgg                                                        >  1:210722/1‑119 (MQ=255)                                                                                                                                  aaGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGACCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTg                  <  1:6218/149‑1 (MQ=255)                                                                                                                                     gTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCac               >  8:34346/1‑148 (MQ=255)                                                                                                                                          ctgCAGCCAGTAGCGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGtt                   >  2:209033/2‑140 (MQ=255)                                                                                                                                          ctgCAGCCAGTAGCGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGtt                   <  1:209033/139‑1 (MQ=255)                                                                                                                                            gCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTg                              >  4:132886/1‑127 (MQ=255)                                                                                                                                            gCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTg                              <  3:132886/127‑1 (MQ=255)                                                                                                                                              aGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTg      <  5:162853/149‑1 (MQ=255)                                                                                                                                                ccGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTGCCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCca                                                     >  7:102999/1‑100 (MQ=255)                                                                                                                                                ccGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCca                                                     <  8:102999/100‑1 (MQ=255)                                                                                                                                                 cGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGt   <  7:109715/149‑1 (MQ=255)                                                                                                                                                  gCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTg  >  4:176394/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GTGCCAAGTTCGCAGCCTTGCTCGATCGCTTCGCAAATGACGTACAGGTTATTACCCAGCCTTGCCCGGGGCTGGTCGAGCTGATCGAAGCCGGTGACCTGGTAAGCCCGGCGGTGCGCCAGCTGCTGCAAGGTTATGTGCAGCCGCTACTGGCTGCAGGTTGCGATACGCTGATCCTCGGCTGCACCCATTACCCTTTCCTTCGCCCATTGCTGGCCGGCATGGTCCCGGATGACGTGGCCATCATCGACACCGGCGCGGCCGTGGCGCGTCAGTTGCAACGGCTGTTGGGTG  >  NC_002947/855523‑855816

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