Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,345,549	C→A	7.1%	P316Q (CCG→CAG)	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,549	0	C	A	7.1%	195.2 / 3.9	65	P316Q (CCG→CAG)	acsA‑II	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (26/34);  new base A (3/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

CCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGTGA  >  NC_002947/5345402‑5345690                                                                                                                                                    |                                                                                                                                                  ccAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCg                                                                                                                                                  >  4:33914/1‑149 (MQ=255) ccAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCg                                                                                                                                                  >  2:352868/1‑149 (MQ=255)    aaGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTg                                                                                                                                               <  6:337459/149‑1 (MQ=255)         gTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCaa                                                                                                                                          <  5:103549/149‑1 (MQ=255)            ctgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATg                                                                                                                                        <  7:304720/148‑1 (MQ=255)            ctgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCAGCGGGCCataga                                                                                                                                       <  8:206121/149‑5 (MQ=255)             tgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGgc                                                                                                                                      >  3:275642/1‑149 (MQ=255)                      ccGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCg                                                                                                                            <  1:352868/149‑1 (MQ=255)                       cggcggCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGc                                                                                                                           <  3:72593/149‑1 (MQ=255)                       cggcggCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGc                                                                                                                           >  2:250569/1‑149 (MQ=255)                        ggcggcTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCt                                                                                                                          <  4:275642/149‑1 (MQ=255)                        ggcggcTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCt                                                                                                                          <  2:42525/149‑1 (MQ=255)                        ggcggcTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCt                                                                                                                          <  2:254506/149‑1 (MQ=255)                                   gctgcAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAg                                                                                                               >  5:354500/1‑149 (MQ=255)                                    ctgcAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGaa                                                                                                                <  3:33914/147‑1 (MQ=255)                                    ctgcAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAg                                                                                                               <  4:370357/148‑1 (MQ=255)                                    ctgcAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCGAGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAg                                                                                                               <  6:358348/148‑1 (MQ=255)                                               gATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAAc                                                                                                    >  3:99812/1‑148 (MQ=255)                                                    ccTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTAcccc                                                                                              <  2:248180/149‑1 (MQ=255)                                                     cTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTAcccc                                                                                              >  5:214361/1‑148 (MQ=255)                                                      ttCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCg                                                                                             >  5:331820/1‑148 (MQ=255)                                                      ttCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGa                                                                                            >  6:13524/1‑149 (MQ=255)                                                      ttCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATATTCGAAGGCGTGCCCAACTACCCCGa                                                                                            >  5:303031/1‑149 (MQ=255)                                                      ttCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCCGATGTTCGAAGGCGTGCCCAACTACCCCGa                                                                                            >  5:69400/1‑149 (MQ=255)                                                        gaaGGTGGTGTTCGTTTATCGTGACGGCGAGGTGTTCTGGTGCAGCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GGAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGa                                                                                            <  5:70496/146‑1 (MQ=255)                                                        cAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCg                                                                                                                            >  2:18593/1‑115 (MQ=255)                                                        cAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGa                                                                                            <  4:159775/147‑1 (MQ=255)                                                         aaGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑Ac                                                                                      <  5:122701/149‑1 (MQ=255)                                                         aaGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑Ac                                                                                      <  7:63827/149‑1 (MQ=255)                                                            gtggtgTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTc                                                                                   <  5:358153/149‑1 (MQ=255)                                                               gtgtTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTcgcg                                                                                >  4:140213/1‑149 (MQ=255)                                                                 gtTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCt                                                                              >  7:143262/1‑149 (MQ=255)                                                                     gATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCt                                                                           <  6:344358/148‑1 (MQ=255)                                                                       ttATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCaa                                                                                                     >  6:156417/1‑123 (MQ=255)                                                                       ttATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCaa                                                                                                     <  5:156417/123‑1 (MQ=255)                                                                       ttATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTgg                                                                         >  1:55426/1‑148 (MQ=255)                                                                            gTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAAAGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTACTGGCAggtg                                                                   >  2:307468/1‑149 (MQ=255)                                                                             tGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAggtgg                                                                  <  5:13524/149‑1 (MQ=255)                                                                             tGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAggtgg                                                                  >  2:241046/1‑149 (MQ=255)                                                                             tGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCAAATGTC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAggtgg                                                                  >  7:238051/1‑149 (MQ=255)                                                                               aCGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAggtggt                                                                 <  1:250569/148‑1 (MQ=255)                                                                                cGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGtggtgg                                                               <  6:217169/149‑1 (MQ=255)                                                                                 ggCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGa                                                              >  8:182756/1‑149 (MQ=255)                                                                                   cGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACa                                                            <  1:241046/149‑1 (MQ=255)                                                                                      ggagtTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGc                                                         <  3:200397/146‑1 (MQ=255)                                                                                        tgtTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCAc                                                       <  4:99812/149‑1 (MQ=255)                                                                                            cTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAgg                                                   >  2:164273/1‑149 (MQ=255)                                                                                               gTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGt                                                  <  6:214361/147‑1 (MQ=255)                                                                                                  cACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGGCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACa                                             >  1:366389/1‑149 (MQ=255)                                                                                                        cgATGTCGGCTGGGGGACCGGCCACATCTACATCCTCTAATGCCCGCTTGCCAATTGCTGGAATCTCGCTGATTTTCGAATGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACAtcttc                                        <  1:308002/149‑1 (MQ=255)                                                                                                                      tGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGAcc                         >  6:328925/1‑149 (MQ=255)                                                                                                                       gACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCg                        >  3:331781/1‑149 (MQ=255)                                                                                                                               acaGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCTTGCCCAACTACCCAGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCATgcgc                <  3:73438/149‑1 (MQ=255)                                                                                                                                caGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGATCTCGCCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTgc                  <  1:307468/149‑1 (MQ=38)                                                                                                                                   cTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGCCCGCCCTGCGCGCa             >  8:141192/1‑148 (MQ=255)                                                                                                                                   cTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCAc            <  6:331820/149‑1 (MQ=255)                                                                                                                                   cTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCAc            <  6:105363/149‑1 (MQ=255)                                                                                                                                     aCATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGTACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACt           >  4:284786/1‑148 (MQ=255)                                                                                                                                         cGTCTATGGACAGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGc      >  4:184333/1‑149 (MQ=255)                                                                                                                                          gccattGGCCAGCGGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTATCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATg       <  1:20264/142‑1 (MQ=255)                                                                                                                                          gccattGGCCAGCGGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTATCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATg       >  2:20264/6‑147 (MQ=255)                                                                                                                                          gTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGc      <  3:140213/148‑1 (MQ=255)                                                                                                                                           tCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGt    <  5:138294/149‑1 (MQ=255)                                                                                                                                            cattGGCCAGCGGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGt    >  8:73503/4‑148 (MQ=255)                                                                                                                                               tGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGTGa  <  7:182756/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGC‑GCAATCTCGCTGATGTTCGAAGGCGTGCCCAACTACCCCGAT‑‑‑ACCTCGCGCTTCTGGCAGGTGGTGGACAAGCACCAGGTGAACATCTTCTACACCGCCCCGACCGCCCTGCGCGCACTGATGCGTGA  >  NC_002947/5345402‑5345690

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