Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,679,138	A→G	10.0%	D255G (GAC→GGC)	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,138	0	A	G	10.0%	162.1 / 7.8	60	D255G (GAC→GGC)	bioA	adenosylmethionine‑‑8‑amino‑7‑oxononanoate aminotransferase BioA
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base G (4/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GACCCTGGCCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTAC  >  NC_002947/5678992‑5679283                                                                                                                                                    |                                                                                                                                                  gACCCTGGCCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGCGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCacgac                                                                                                                                                  >  2:95053/1‑148 (MQ=255)        gCCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCacga                                                                                                                                                   <  4:127822/140‑1 (MQ=255)        gCCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCacga                                                                                                                                                   >  3:127822/1‑140 (MQ=255)        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                                                                                                                                              <  1:240996/136‑1 (MQ=255)                 cacGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACgaga                                                                                                                                              >  2:240996/1‑136 (MQ=255)                     ccTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCAcgacga                                                                                                                                                >  6:89298/1‑130 (MQ=255)                     ccTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCAcgacga                                                                                                                                                <  5:89298/130‑1 (MQ=255)                         gATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCg                                                                                                                           <  6:24878/147‑1 (MQ=255)                          aTCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGGCGATCTCGTCGTggata                                                                                                                                  >  4:39447/1‑136 (MQ=255)                          aTCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGGCGATCTCGTCGTggata                                                                                                                                  <  3:39447/139‑4 (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                                                                                             <  7:282155/123‑1 (MQ=255)                                                                               gTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGc                                                                                             >  8:282155/1‑123 (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                                                                    <  7:183248/130‑1 (MQ=255)                                                                                                  cAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGc                                                                    >  8:183248/1‑130 (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                                                                                                    >  4:139261/1‑80 (MQ=255)                                                                                                                    cTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCg                                                                                                    <  3:139261/80‑1 (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                                       <  2:99910/114‑1 (MQ=255)                                                                                                                                               ccACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCt                                       >  1:99910/1‑114 (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)                                                                                                                                                    |                                                                                                                                                  GACCCTGGCCGAACACCACGCCTCGATCAGCGCCGTGATCGTCGAGCCGCTGATACAGGGTGCCGGTGGCATGCGCATGTACCACCCGGTGT‑ACCTCAAGTTGTTGCGCGAGGCCTGCGACCGCTATGACGTGCACCTTATCCACGACGAGATCGCCGTGGGCTTCGGCCGCACCGGCACGATGTTCGCCTGCGAGCAAGCCGGCATCCGCCCGGATTTCCTGTGCCTGTCCAAGGCCCTGACCGGCGGCTACCTGCCGCTGGCCGCCTGCCTGACCACCGACAAGGTGTAC  >  NC_002947/5678992‑5679283

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