Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,654,726	T→C	5.5%	F402L (TTC→CTC)	aceA →	isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,654,726	0	T	C	5.5%	214.8 / 2.1	73	F402L (TTC→CTC)	aceA	isocitrate lyase
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base C (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

CAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAG  >  NC_002947/4654581‑4654871                                                                                                                                                  |                                                                                                                                                  cAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTc                                                                                                                                                 <  5:411341/148‑1 (MQ=255)    gCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCCACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGACTTCGCTTACGCCCGCAAGGGCTACACCTTCGTgg                                                                                                                                             >  2:267868/1‑149 (MQ=255)      aTCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGgcg                                                                                                                                           <  4:125241/149‑1 (MQ=255)         aGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCAc                                                                                                                                        >  2:89412/1‑149 (MQ=255)             catcaCCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACcag                                                                                                                                     <  7:393682/148‑1 (MQ=255)              atcaCCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACcag                                                                                                                                     >  6:15340/1‑147 (MQ=255)              atcaCCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACcag                                                                                                                                     >  2:24359/1‑147 (MQ=255)                caCCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTc                                                                                                                                                 >  4:31034/1‑133 (MQ=255)                caCCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTc                                                                                                                                                 <  3:31034/133‑1 (MQ=255)                 aCCCTGGCCCGCAACCACCACATGTGGCACGGCATGAACAACCTGGCGCACGACTACGCCCGCACCGACAAGCCCGCCAACGAGACGCTGCAGGAGCAGGACTTCGCTGACGCCCGCAAGGGCTACACCTTCGTGGCGCACCAGCCgg                                                                                                                                 >  4:348398/1‑148 (MQ=255)                  cccTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCACCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAgg                                                                                                                                 >  4:65157/1‑147 (MQ=255)                  cccTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCACCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAgg                                                                                                                                 <  3:65157/147‑1 (MQ=255)                    cTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGt                                                                                                                             <  3:57954/149‑1 (MQ=255)                       gCCGGCATCCACAACTTGTGGCACGGCATGTTCATCCTTGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGGTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTggg                                                                                                                          <  4:86899/149‑1 (MQ=255)                       gCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTg                                                                                                                            <  2:403911/147‑1 (MQ=255)                                 acaacaTGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTAc                                                                                                                >  3:332229/1‑149 (MQ=255)                                  caacaTGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACt                                                                                                               <  8:329385/149‑1 (MQ=255)                                   aacaTGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACtt                                                                                                              <  5:15340/149‑1 (MQ=255)                                      aTGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTcga                                                                                                           >  1:218232/1‑149 (MQ=255)                                       tgtgGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTcgac                                                                                                          <  1:99146/149‑1 (MQ=255)                                        gtgGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTcga                                                                                                           <  5:285189/147‑1 (MQ=255)                                        gtgGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTcga                                                                                                           >  6:285189/1‑147 (MQ=255)                                           gCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACa                                                                                                      >  2:423769/1‑149 (MQ=255)                                              cGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGa                                                                                                   >  5:228214/1‑149 (MQ=255)                                              cGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTCGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGa                                                                                                   <  1:301022/149‑1 (MQ=255)                                                      tCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCg                                                                                                                                                >  2:217857/1‑96 (MQ=255)                                                        aaCCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCATCAAGGTCTACACCTTCGTGGCGCACCATCAGGAAGTGGCCACTGTCTACTTCGACGACATGACCACCGTGa                                                                                          <  1:267868/148‑1 (MQ=255)                                                          ccTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATc                                                                                        >  7:290207/1‑148 (MQ=255)                                                           cTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCa                                                                                      <  8:174425/149‑1 (MQ=255)                                                                  aCGCCTACGCCCGCAACGACATGACCGCCTTCGTGAATCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGAACACCGTGATCCAGggcggc                                                                               <  7:361730/149‑1 (MQ=255)                                                                      cTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCtt                                                                           <  5:441679/149‑1 (MQ=255)                                                                       tACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAACTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTgg                                                                                                                           <  5:30963/100‑1 (MQ=255)                                                                       tACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAAATCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTgg                                                                                                                           >  6:30963/1‑100 (MQ=255)                                                                           cccGCCACGCCAGGCGCGCCTACGTGAAGCTGCAGGCGCAGGAATTCGGGGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGGGGGCACTGGCTACTTCGACGTCAGGCCCACCGTGATCCAGGGCGGCGCTtcgtcg                                                                      <  3:348398/149‑1 (MQ=255)                                                                            ccGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCgg                                                                     >  7:391625/1‑149 (MQ=255)                                                                                  cGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTg                                                               <  8:391625/149‑1 (MQ=255)                                                                                      aTGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCACGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACt                                                           <  6:121719/149‑1 (MQ=255)                                                                                         aCCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGAc                                                        >  1:434603/1‑149 (MQ=255)                                                                                         aCCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGAc                                                        <  4:332229/149‑1 (MQ=255)                                                                                             ccTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACt                                                                <  6:132435/137‑1 (MQ=255)                                                                                             ccTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACt                                                                >  5:132435/1‑137 (MQ=255)                                                                                             ccTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGt                                                    <  1:24359/149‑1 (MQ=255)                                                                                             ccTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGt                                                    >  7:136479/1‑149 (MQ=255)                                                                                              cTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGtt                                                   >  6:142033/1‑149 (MQ=255)                                                                                                    gAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGAcc                                              >  1:147432/1‑148 (MQ=255)                                                                                                         tGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCgaag                                          >  3:329125/1‑147 (MQ=255)                                                                                                          gcaggagcaggaATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCgaag                                          >  3:287588/1‑146 (MQ=255)                                                                                                          ccaggagcagcaATTAGCTGACCCCAGCAACTGCTACACCTTCGTGGCGCACAAGAAGGAACTGATCAATGGCTACTTCGAAGACATGACAACAGTGATACAGTGAGCCGCTTCGTCTGTGACTGCACTGACCGCTTCGACCgaataa                                        >  4:49619/2‑148 (MQ=255)                                                                                                             ggagcaggaATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGc                                    <  4:287588/149‑1 (MQ=255)                                                                                                              gagcaggaATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAagag                                     >  1:212200/1‑147 (MQ=255)                                                                                                               agcaggaATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAg                                  <  2:1406/149‑1 (MQ=255)                                                                                                                  atgaATTCGCTGACGCCAGCAAGGTCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGTCTACTTCGACGACATGACCACCGTGATCCAGGTCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTc                               <  1:167322/147‑1 (MQ=255)                                                                                                                  aggaATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTc                               >  3:427452/1‑149 (MQ=255)                                                                                                                    gaATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTtc                                                                          <  2:43933/104‑1 (MQ=255)                                                                                                                    gaATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTtc                                                                          >  1:43933/1‑104 (MQ=255)                                                                                                                    gaATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCa                             >  1:27710/1‑149 (MQ=255)                                                                                                                      aTTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCAc                            <  2:434603/148‑1 (MQ=255)                                                                                                                          gCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGa                         >  7:81989/1‑147 (MQ=255)                                                                                                                          gCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATa                       <  3:163074/149‑1 (MQ=255)                                                                                                                           cTGACGCCAGCAAGGGCTACACCCTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAg                      >  4:317489/1‑149 (MQ=255)                                                                                                                           cTGACGCCAGCAAGGGCTACACCCTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAg                      >  4:311857/1‑149 (MQ=255)                                                                                                                            tGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGt                     >  7:95118/1‑149 (MQ=255)                                                                                                                            tGACGCCAGCAAGAGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCTGTGCCTGCACTGACCGGTTCGAGCGAAGAAGAGCTGTTGCAGTGTTAGt                     >  7:371879/1‑149 (MQ=255)                                                                                                                              aCGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGt                     <  2:27710/147‑1 (MQ=255)                                                                                                                                gCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGt                 >  8:154108/1‑149 (MQ=255)                                                                                                                                gCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGGGATCCAGGGCGGCGCTTCGTCGGGGACTGCACTGACCGGTTCGCCCGAAGAAGAGCAGTTCGACTGCTAGTgggg                 >  4:336565/1‑148 (MQ=255)                                                                                                                                 ccAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGtt                >  2:62939/1‑149 (MQ=255)                                                                                                                                         ggCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGa        <  3:49619/149‑1 (MQ=255)                                                                                                                                           cTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGaa       >  7:62680/1‑148 (MQ=255)                                                                                                                                              cacCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACaa    <  6:228214/148‑1 (MQ=255)                                                                                                                                               acCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAg  >  7:390577/1‑149 (MQ=255)                                                                                                                                               acCCTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAg  <  3:311857/149‑1 (MQ=255)                                                                                                                                               acCCTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAg  <  3:317489/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCAACCTGGCGCACGACTACGCCCGCAACGACATGACCGCCTACGTGAAGCTGCAGGAGCAGGAATTCGCTGACGCCAGCAAGGGCTACACCTTCGTGGCGCACCAGCAGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAG  >  NC_002947/4654581‑4654871

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