Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,679,143	A→C	9.5%	I257L (ATC→CTC)	bioA →	adenosylmethionine‑‑8‑amino‑7‑oxononanoate aminotransferase BioA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,679,143	0	A	C	9.5%	184.2 / 9.5	63	I257L (ATC→CTC)	bioA	adenosylmethionine‑‑8‑amino‑7‑oxononanoate aminotransferase BioA
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base C (4/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GCCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTAC  >  NC_002947/5678999‑5679283                                                                                                                                                  |                                                                                                                                             gCCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATc                                                                                                                                            <  3:1934/147‑1 (MQ=255)    gAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGt                                                                                                                                       >  2:207723/1‑149 (MQ=255)      acacCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTgg                                                                                                                                     <  7:176087/149‑1 (MQ=255)      acacCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTgg                                                                                                                                     <  1:95053/149‑1 (MQ=255)         ccacGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCCCGGCGATCTCGTCGTggata                                                                                                                                  <  4:208649/149‑4 (MQ=255)          cacGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACgaga                                                                                                                                              >  2:240996/1‑136 (MQ=255)          cacGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACgaga                                                                                                                                              <  1:240996/136‑1 (MQ=255)                  gATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCg                                                                                                                           <  6:24878/147‑1 (MQ=255)                   aTCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGGCGATCTCGTCGTggata                                                                                                                                  <  3:39447/139‑4 (MQ=255)                   aTCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGGCGATCTCGTCGTggata                                                                                                                                  >  4:39447/1‑136 (MQ=255)                           cGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGa                                                                                                                <  5:123348/149‑1 (MQ=255)                              gATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGt                                                                                                             >  4:77805/1‑149 (MQ=255)                               aTCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGt                                                                                                             <  5:237282/148‑1 (MQ=255)                                    cgcgCCGCTGATACCGGGTGCCGGTGGCAGGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGCCGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCt                                                                                                       <  7:131448/146‑1 (MQ=255)                                      aGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGTGACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTg                                                                                                      >  1:233657/1‑149 (MQ=255)                                      aGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGc                                                                                                     <  1:207723/149‑1 (MQ=255)                                                aCAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCg                                                                                           >  1:283641/1‑149 (MQ=255)                                                                     cATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCtgt                                                                      <  2:233657/149‑1 (MQ=255)                                                                      ctGTACCACCCTGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGAttt                                                                           <  6:322100/142‑1 (MQ=255)                                                                      aTGTACCACACGGTGT‑ACCTCAAGTTGTTGCTCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGAGAGCAAGCCGGCATCCGCCCGGAttt                                                                           >  5:322100/1‑143 (MQ=255)                                                                        gTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGc                                                                                             >  8:282155/1‑123 (MQ=255)                                                                        gTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGc                                                                                             <  7:282155/123‑1 (MQ=255)                                                                           ccaccCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCt                                                                        <  2:306637/141‑1 (MQ=255)                                                                           ccaccCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCt                                                                        >  1:306637/1‑141 (MQ=255)                                                                            caccCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTc                                                               >  5:292902/1‑149 (MQ=255)                                                                            caccCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTc                                                               >  5:324228/1‑149 (MQ=255)                                                                              cccGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCtgt                                                                      >  5:180367/1‑140 (MQ=255)                                                                              cccGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCtgt                                                                      <  6:180367/140‑1 (MQ=255)                                                                                       aCCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGccc                                                                                 <  1:139851/121‑1 (MQ=255)                                                                                       aCCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGccc                                                                                 >  2:139851/1‑121 (MQ=255)                                                                                        ccTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGa                                                    >  4:236820/1‑149 (MQ=255)                                                                                         cTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGAc                                                   >  3:255575/1‑149 (MQ=255)                                                                                          tCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGAcc                                                  >  3:160294/1‑149 (MQ=255)                                                                                           cAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGc                                                                    >  8:183248/1‑130 (MQ=255)                                                                                           cAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGc                                                                    <  7:183248/130‑1 (MQ=255)                                                                                               ttgttgCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACcggcgg                                             <  6:324228/149‑1 (MQ=255)                                                                                                 gttgCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCACCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCggcggc                                            <  4:282250/148‑1 (MQ=255)                                                                                                      gcgAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTg                                      <  3:236820/149‑1 (MQ=255)                                                                                                       cgAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCTTCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGc                                     >  5:126462/1‑149 (MQ=255)                                                                                                         aGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCg                                   <  6:292902/149‑1 (MQ=255)                                                                                                             cTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCg                                                                                                    <  3:139261/80‑1 (MQ=255)                                                                                                             cTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCg                                                                                                    >  4:139261/1‑80 (MQ=255)                                                                                                             cTGCGACCGCTATGACGTGCACCATATCCACCACGACATCGCCCTCAGCTTCGGCCGCACCCGCACCAACTTCGCCTGCGAGCAAGCCGGCATCCGCCCCTAAATCCTGTGCCAGTCCGAGGCCCTGACCGGAGGCTACCTGCCGCAgg                               >  6:82840/1‑149 (MQ=255)                                                                                                              tGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGgc                              >  7:88062/1‑149 (MQ=255)                                                                                                                   ccGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCtgt                                                                      >  2:129803/1‑104 (MQ=255)                                                                                                                   ccGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCtgt                                                                      <  1:129803/104‑1 (MQ=255)                                                                                                                    cGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCgcctg                        >  2:206891/1‑149 (MQ=255)                                                                                                                      cTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGACCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCgcctgcc                      >  3:162412/1‑149 (MQ=255)                                                                                                                       tATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCgcctgcct                     >  4:44107/1‑149 (MQ=255)                                                                                                                         tGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGa                   >  7:97910/1‑149 (MQ=255)                                                                                                                             gTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGCCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGaccac               <  4:162412/149‑1 (MQ=255)                                                                                                                               gCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGaccacc              >  6:233095/1‑148 (MQ=255)                                                                                                                                   gaagcccACGGCGATCTCGTCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGAc           >  2:308403/6‑147 (MQ=255)                                                                                                                                      aTCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCGAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAg        <  5:82840/147‑1 (MQ=255)                                                                                                                                        ccACGGCGATCTCGTCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGgtg     >  1:50337/1‑148 (MQ=255)                                                                                                                                        ccACGGCGATCTCGTCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCGGCGAGCAAGCCGGTATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGgtg     >  6:259988/1‑148 (MQ=255)                                                                                                                                        ccACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCt                                       >  1:99910/1‑114 (MQ=255)                                                                                                                                        ccACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCt                                       <  2:99910/114‑1 (MQ=255)                                                                                                                                        ccACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTACCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGgtg     >  1:99971/1‑148 (MQ=255)                                                                                                                                         cACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGgtgt    <  3:77805/148‑1 (MQ=255)                                                                                                                                          acgacgAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCTTCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTAc  <  6:126462/149‑1 (MQ=255)                                                                                                                                          acgacgAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGaccac               <  3:151328/136‑1 (MQ=255)                                                                                                                                          acgacgAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGaccac               >  4:151328/1‑136 (MQ=255)                                                                                                                                            gacgaGATCGCCGTGGGCTTCGGCCGCAGCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAg        >  1:177487/1‑141 (MQ=255)                                                                                                                                            gacgaGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAg        <  2:177487/141‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                             GCCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTAC  >  NC_002947/5678999‑5679283

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