Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,345,974	G→A	7.3%	G458S (GGC→AGC)	acsA‑II →	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,345,974	0	G	A	7.3%	198.2 / 2.0	67	G458S (GGC→AGC)	acsA‑II	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base A (2/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGGCAGACCGAAACCGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTGCCCGTCGCGATGCCGACGGC  >  NC_002947/5345827‑5346114                                                                                                                                                    |                                                                                                                                             tggtggCAGACCGAAACCGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCccggcgccgg                                                                                                                                             >  3:117956/1‑149 (MQ=255)    tggCAGACCGAAACCGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGt                                                                                                                                          <  1:383183/149‑1 (MQ=255)            cGAAACCGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAAACCGGCTGCGCCACCCAGACGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCAGCACGCCGGCgc                                                                                                                                  >  6:28231/1‑148 (MQ=255)              aaaCCGGCGGCATCATGCTCACCCCTCTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATc                                                                                                                                <  1:400319/149‑1 (MQ=255)              aaaCCGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATc                                                                                                                                <  8:264138/149‑1 (MQ=255)               aaCCGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCAGCACGCCGGCgccgg                                                                                                                               <  5:28231/149‑5 (MQ=255)                aCCGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCAGCACGCCGGCgccggg                                                                                                                              <  6:96075/149‑6 (MQ=255)                 ccGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCaaa                                                                                                                             >  3:379469/1‑149 (MQ=255)                 ccGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCaaa                                                                                                                             >  3:287852/1‑149 (MQ=255)                 ccGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCaaa                                                                                                                             <  6:97961/149‑1 (MQ=255)                 ccGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCa                                                                                                                               >  5:97961/1‑147 (MQ=255)                     cggcATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCa                                                                                                                         <  2:71385/149‑1 (MQ=255)                        catcatGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCt                                                                                                                      <  4:287852/149‑1 (MQ=255)                        catcatGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCt                                                                                                                      >  5:287787/1‑149 (MQ=255)                        catcatGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCt                                                                                                                      >  8:192915/1‑149 (MQ=255)                            atGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGc                                                                                                                   <  6:26702/148‑1 (MQ=255)                                  aCCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCaaa                                                                                                                             >  6:181871/1‑132 (MQ=255)                                  aCCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCaaa                                                                                                                             <  5:181871/132‑1 (MQ=255)                                  aCCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGc                                                                                                             <  2:339763/148‑1 (MQ=255)                                      cACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATc                                                                                                        >  1:367081/1‑149 (MQ=255)                                           cccGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGc                                                                                                     <  6:287787/147‑1 (MQ=255)                                            ccGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGc                                                                                                  <  7:192915/149‑1 (MQ=255)                                              ggCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGt                                                                                                >  3:161749/1‑149 (MQ=255)                                                cacacaGTCACTCAAACCCGTCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCt                                                                                              <  5:410424/149‑1 (MQ=255)                                                         aCTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGac                                                                                      >  3:261625/1‑148 (MQ=255)                                                           tCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCAGCACGCCGGCgccg                                                                                                                                >  7:163713/1‑101 (MQ=255)                                                           tCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCAGCACGCCGGCgccg                                                                                                                                <  8:163713/104‑4 (MQ=255)                                                             aaaCCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGAccacca                                                                                 >  2:402597/1‑149 (MQ=255)                                                                           cACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGac                                                                    <  2:182738/148‑1 (MQ=255)                                                                           cACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGac                                                                    >  1:182738/1‑148 (MQ=255)                                                                            aCCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGacac                                                                  >  8:11036/1‑149 (MQ=255)                                                                             cccAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACAcc                                                                 <  7:309857/149‑1 (MQ=255)                                                                                   cGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCGCGGGCAGATCCGCATATTCTACTGCGACGACCAGCGGATGGTCGACACCTACTTc                                                           <  4:379469/149‑1 (MQ=255)                                                                                          cGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGcccc                                                    >  7:223689/1‑149 (MQ=255)                                                                                          cGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGcccc                                                    >  5:62051/1‑149 (MQ=255)                                                                                          cGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTCGCGATCAACGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGcccc                                                    >  6:413978/1‑149 (MQ=255)                                                                                            gCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGc                                                                                                                                    >  5:275658/1‑67 (MQ=255)                                                                                            gCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGc                                                                                                                                    <  6:275658/67‑1 (MQ=255)                                                                                              gTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTAc                                                                                            >  1:44961/1‑105 (MQ=255)                                                                                              gTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTAc                                                                                            <  2:44961/105‑1 (MQ=255)                                                                                                gCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGacc                                                                                     <  4:292106/110‑1 (MQ=255)                                                                                                gCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGacc                                                                                     >  3:292106/1‑110 (MQ=255)                                                                                                  aaCCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGc                                                                                         <  5:273237/104‑1 (MQ=255)                                                                                                  aaCCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGc                                                                                         >  6:273237/1‑104 (MQ=255)                                                                                                   aCCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTc                                                                       >  5:11928/1‑121 (MQ=255)                                                                                                   aCCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTc                                                                       <  6:11928/121‑1 (MQ=255)                                                                                                     cGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGctac                                         <  6:62051/149‑1 (MQ=255)                                                                                                     cGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGctac                                         <  1:402597/149‑1 (MQ=255)                                                                                                             ctgGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCAc                                 >  5:24088/1‑149 (MQ=255)                                                                                                              tgGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACt                                >  1:174420/1‑149 (MQ=255)                                                                                                               ggACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTg                               >  7:247328/1‑149 (MQ=255)                                                                                                                gACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTgg                              <  8:223689/149‑1 (MQ=255)                                                                                                                   gAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGa                           >  2:276155/1‑149 (MQ=255)                                                                                                                   gAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGa                           >  4:137093/1‑149 (MQ=255)                                                                                                                     aaaaGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTa                                                               <  2:4507/111‑1 (MQ=255)                                                                                                                     aaaaGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTa                                                               >  1:4507/1‑111 (MQ=255)                                                                                                                      aaaGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACgg                        >  5:403254/1‑149 (MQ=255)                                                                                                                        aGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTg                      <  2:174420/149‑1 (MQ=255)                                                                                                                         ggCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCa                                  >  2:391109/1‑136 (MQ=255)                                                                                                                         ggCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCa                                  <  1:391109/136‑1 (MQ=255)                                                                                                                         ggCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTGc                     <  6:403254/149‑1 (MQ=255)                                                                                                                           cAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACCGCGACCACCAGCCGATGGTCGACACCTTCTACAAGCCCCTGCCCCGCAACTACTTCCCTGGCGAcaggcacc                   >  5:122199/1‑142 (MQ=255)                                                                                                                                  aTCGAAGGCCCCGGCGCCGGCGGGCGGGCGAGCATCGCGAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGCTGGGCTACACCTCCTTCACGCCCCTGCCCGGCTACTCCTTCACTGGCGACGGTGCCCGTCGGGa            <  8:95447/149‑1 (MQ=255)                                                                                                                                   tCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTGCCCGTcg               >  7:217923/1‑145 (MQ=255)                                                                                                                                   tCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTGCCCGTcg               <  8:217923/145‑1 (MQ=255)                                                                                                                                   tCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTGCCCGTCGCGAt           >  4:312807/1‑149 (MQ=255)                                                                                                                                   tCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTGCCCGTCGCGAt           <  4:168935/149‑1 (MQ=255)                                                                                                                                     gAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTGCCCGTCGCGATGc         >  5:18236/1‑149 (MQ=255)                                                                                                                                            ccggcgccggcgTGCTGGCGATCAAAGCCAGCTGGCCCGTGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTGCCCGTCGCGATGCCGACGGc  <  1:242406/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                             TGGTGGCAGACCGAAACCGGCGGCATCATGCTCACCCCACTGCCCGGCACACAGTCACTCAAACCCGGCTGCGCCACCCAGCCGATGTTCGGCGTGCAACCGGTGCTGCTGGACGAAAAAGGCAAACTGATCGAAGGCCCCGGCGCCGGCGTGCTGGCGATCAAAGCCAGCTGGCCCGGGCAGATCCGCAGCGTCTACGGCGACCACCAGCGGATGGTCGACACCTACTTCAAGCCCCTGCCCGGCTACTACTTCACTGGCGACGGTGCCCGTCGCGATGCCGACGGC  >  NC_002947/5345827‑5346114

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