Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,535,979	C→A	6.8%	A332A (GCC→GCA)	pncB →	nicotinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,535,979	0	C	A	6.8%	176.4 / 3.0	59	A332A (GCC→GCA)	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base A (2/2);  total (26/33)
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.01e-01

CCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGAC  >  NC_002947/5535835‑5536125                                                                                                                                                  |                                                                                                                                                   ccTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgc                                                                                                                                                <  4:283968/149‑1 (MQ=255) ccTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgc                                                                                                                                                <  6:22662/149‑1 (MQ=255)  cTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgcg                                                                                                                                               >  4:253329/1‑149 (MQ=255)      gcCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGa                                                                                                                                           <  8:248629/149‑1 (MQ=255)       ccATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGcc                                                                                                                                                    <  4:143343/139‑1 (MQ=255)       ccATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGcc                                                                                                                                                    >  3:143343/1‑139 (MQ=255)       ccATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGAc                                                                                                                                          >  2:122665/1‑149 (MQ=255)        cATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCAGGGCGCGGaaa                                                                                                                                         <  4:351198/149‑3 (MQ=255)        cATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCAGGGCGCGGaa                                                                                                                                          >  3:70646/1‑147 (MQ=255)            aTTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCAGGGCGCGGaaaatt                                                                                                                                      <  4:70646/148‑6 (MQ=255)               cGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCGATTGGCGATTACCAGATACTGGGTATCGACCCGATGACCAA‑GCCGCTGGGGTTTTCCGGCGGCGTGAACCTGACCCGCTCGCTGGAGATGTTCCGCGCCCTGCGCGGCCGGATCCAc                                                                                                                                  >  4:58454/1‑149 (MQ=255)                  gCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTc                                                                                                                               <  4:83012/149‑1 (MQ=255)                      gCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCt                                                                                                                           <  6:124793/149‑1 (MQ=255)                        cGGTAGCCTGGGCGGAGAAGGCCAATTCCCATTAGCAGAAACTGGGTATCGGCCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCttt                                                                                                                         <  4:85461/149‑1 (MQ=255)                         ggTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTg                                                                                                                        >  7:205227/1‑149 (MQ=255)                         ggTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTg                                                                                                                        <  6:2969/149‑1 (MQ=255)                               cTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCg                                                                                                                  <  6:243879/149‑1 (MQ=255)                                  ggCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCgg                                                                                                                 <  3:53577/147‑1 (MQ=255)                                  ggCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCgg                                                                                                                 <  1:292086/147‑1 (MQ=255)                                       agaAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGGACCCTGGTGTTTTCCGACGGCCACAACCTGACCCGCACGCAGGAAATTTACCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCaa                                                                                                           >  1:283874/1‑149 (MQ=255)                                       agaAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAac                                                                                                          >  1:111533/1‑149 (MQ=255)                                       agaAGGCCATTGCCCATTACCACAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAac                                                                                                          <  7:247422/149‑1 (MQ=255)                                         aaGGCCATTGCACATTACCAGAAACTGGGTAGCGACCCGATGACCAA‑GACGCTGGGGTGTTCCGACGGCCTCCACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacct                                                                                                        <  7:333276/149‑1 (MQ=255)                                           ggCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctga                                                                                                      >  6:37542/1‑149 (MQ=255)                                               aTTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgacctg                                                                                                  <  6:87404/149‑1 (MQ=255)                                               aTTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgacctg                                                                                                  <  8:203893/149‑1 (MQ=255)                                                   cccATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGcaa                                                                                               <  2:111533/148‑3 (MQ=255)                                                           cAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAgg                                                                                      >  5:298078/1‑149 (MQ=255)                                                           cAGAAACTGGGTATCGAACCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTCAACAGCTCGCTCGAAATTTTCCGCGCCCTGAGCGGACGCATCAACGTCAGATTTGGCATCGGCACCAAACTCACCTGCGACATCCCAgg                                                                                      >  4:331668/1‑149 (MQ=255)                                                            aGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGgt                                                                                     <  7:124713/149‑1 (MQ=255)                                                            aGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGgt                                                                                     >  5:89826/1‑149 (MQ=255)                                                               aaCTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGgtgt                                                                                   <  2:316453/148‑1 (MQ=255)                                                                    ggTAACGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCCCAGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCATAGGCAACGGCACCAACCAGACCTGCGACATCCCAGGTGTGcccac                                                                             >  5:170433/1‑147 (MQ=255)                                                                           aCCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAAc                                                                      <  1:346450/149‑1 (MQ=255)                                                                             ccGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACa                                                                     >  1:125216/1‑148 (MQ=255)                                                                               gATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGAGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCg                                                                  >  8:22091/1‑149 (MQ=255)                                                                                    ccAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCt                                                              >  1:203433/1‑148 (MQ=255)                                                                                       a‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTaaa                                                          <  3:58454/149‑1 (MQ=255)                                                                                         gACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTTGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTaaaa                                                         <  2:110683/149‑1 (MQ=255)                                                                                            gCTGGTGTTTTCCGACGGCCTCAACCTGTCCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCTTCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCGCATGAACATCGTGCTTAAAATGa                                                      <  2:283874/149‑1 (MQ=255)                                                                                              tGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGGGATGTCGCAGgtcagg                                                                                 >  6:68339/1‑116 (MQ=255)                                                                                                  gtTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCg                                                                                                <  6:9606/101‑1 (MQ=255)                                                                                                  gtTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCg                                                                                                >  5:9606/1‑101 (MQ=255)                                                                                                   ttttCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTg                                               >  5:308015/1‑149 (MQ=255)                                                                                                        cGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAAc                                           <  6:298078/148‑1 (MQ=255)                                                                                                               cTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATccc                                                                                         >  2:123825/1‑95 (MQ=255)                                                                                                               cTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACg                                          <  4:26268/142‑1 (MQ=255)                                                                                                               cTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACg                                          >  3:26268/1‑142 (MQ=255)                                                                                                                tCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCg                                  >  5:168773/1‑149 (MQ=255)                                                                                                                  aaCCTGACCGGCTCTGTGGTATTTTTTCGCGACCTGGGCTGACGCATTATGTTCAGCTTTGTCATCTGTACTTAACTGACCTGCTAAATCCATTGTGTGTGCCCGATGAACATCGTGTTTTACCTGACCGACTGCAACGGCGCGCCGGt                                <  6:380137/149‑1 (MQ=255)                                                                                                                   aCCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTg                               <  7:22091/149‑1 (MQ=255)                                                                                                                     cTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGc                             <  8:205227/149‑1 (MQ=255)                                                                                                                        aCCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGcc                            >  7:146280/1‑147 (MQ=255)                                                                                                                          ccGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGCCGCATCACCGTCAGCTTTGGCATCGGCACCCACCTGACCTGCGCCATCCCAGGTGTGGCCCCGATGACCTGCGTGCGGAAAGTGCCCGCCTGCAACGGCGCGCCGGGGGcc                            >  7:317578/1‑145 (MQ=255)                                                                                                                          ccGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCAAGAAAATCGTGCTTAAAATGACCGACTCCAACGGCGCGCCGGTGGCCaac                         >  8:296890/1‑147 (MQ=255)                                                                                                                                gctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATctct                   <  5:37542/148‑1 (MQ=255)                                                                                                                                gctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATctct                   <  2:125216/148‑1 (MQ=255)                                                                                                                                gctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTg                  <  6:89826/149‑1 (MQ=255)                                                                                                                                  tGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGAc                >  4:89293/1‑149 (MQ=255)                                                                                                                                          gtCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCCTGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgc        >  5:329071/2‑149 (MQ=255)                                                                                                                                                gcCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGTCACCAACCTGACCTGCGACATCCCAGGTGTGGACCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGAc  <  3:331668/149‑1 (MQ=255)                                                                                                                                                  ccTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGgcg                                       <  7:4783/110‑1 (MQ=255)                                                                                                                                                  ccTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGgcg                                       >  8:4783/1‑110 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGAC  >  NC_002947/5535835‑5536125

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