Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	156,383	C→T	7.4%	intergenic (‑8/+35)	ygdQ ← / ← citN	transport protein/citrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	156,383	0	C	T	7.4%	192.8 / 3.5	67	intergenic (‑8/+35)	ygdQ/citN	transport protein/citrate transporter
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base T (2/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

ATGCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCG  >  NC_002947/156236‑156529                                                                                                                                                     |                                                                                                                                                    aTGCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCg                                                                                                                                                    <  5:347968/149‑1 (MQ=255)            cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcg                                                                                                                                         >  2:234955/1‑149 (MQ=255)             agtaTCGAGATCATGATGCTGTTTTCGTTACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCGATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCTGACCGGGGcgcg‑t                                                                                                                                       <  3:67323/146‑1 (MQ=255)              gAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcg‑tc                                                                                                                                      >  7:271972/1‑149 (MQ=255)                aTCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCTGACCGGGGcgcg‑ttgc                                                                                                                                    <  4:117438/149‑1 (MQ=255)                aTCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcg‑tcgc                                                                                                                                    <  1:237382/149‑1 (MQ=255)                 tCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcg‑tcgcg                                                                                                                                   >  2:249780/1‑149 (MQ=255)                 tCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcg‑tcgcg                                                                                                                                   >  7:365943/1‑149 (MQ=255)                      aTCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGt                                                                                                                              <  3:100612/149‑1 (MQ=255)                           gatgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATAc                                                                                                                         >  6:289079/1‑149 (MQ=255)                             tgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACg                                                                                                                        <  6:59313/148‑1 (MQ=255)                               atgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGa                                                                                                                      <  1:234955/148‑1 (MQ=255)                                tgttgtCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCATT‑CGCGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATg                                                                                                                    >  3:351612/1‑149 (MQ=255)                                 gttgtCGATTACGATGACGATCTTAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGAt                                                                                                                     <  1:295724/148‑1 (MQ=255)                                 gttgtCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGACTGGTCAGCCATCCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGt                                                                                                                   >  8:22002/1‑149 (MQ=255)                                  ttgtCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTa                                                                                                                  >  3:342032/1‑149 (MQ=255)                                   tgtCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTa                                                                                                                  <  8:365943/148‑1 (MQ=255)                                      cGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTa                                                                                                               >  1:204381/1‑148 (MQ=255)                                       gaaaCCGTGTACTATCTCAAGCGCGTTAAGGGTCAAAACGCCACCCCTGATTCCCGGGCTGTTCAGCCATTCCATGTATTTTCCTT‑CGTGCGTTAATGTCAACGCGCCCCGGTCAGGGCGCG‑TTGCGTACGTGATACGGATGTATTaaa                                                                                                             <  1:96946/146‑1 (MQ=255)                                       gATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTaaa                                                                                                             <  7:196192/149‑1 (MQ=255)                                          aCCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGac                                                                                                          >  3:121144/1‑149 (MQ=255)                                           ccGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGac                                                                                                          <  8:271972/148‑1 (MQ=255)                                              aGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGactact                                                                                                      <  7:121289/149‑1 (MQ=255)                                                 aCGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGaa                                                                                                   <  6:234090/149‑1 (MQ=255)                                                   gATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGaa                                                                                                   >  6:381548/1‑147 (MQ=255)                                                   gATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACa                                                                                                 >  4:158851/1‑149 (MQ=255)                                                       tcAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTAc                                                                                                                                >  7:364233/1‑114 (MQ=255)                                                       tcAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTAc                                                                                                                                <  8:364233/114‑1 (MQ=255)                                                        cAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCgg                                                                                             >  4:65592/1‑148 (MQ=255)                                                         aaGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCggg                                                                                            <  4:121144/148‑1 (MQ=255)                                                         aaGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTA‑TTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGa                                                                                           >  7:182778/1‑148 (MQ=255)                                                           gcgcGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGaaa                                                                                         <  5:105727/149‑1 (MQ=255)                                                             gcgGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAAt                                                                                       <  3:158851/149‑1 (MQ=255)                                                                 tAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATcccc                                                                                   <  8:275851/149‑1 (MQ=255)                                                                   aGGGTGAAAAAGGCAACCCAGATTTCCGGACTGGTCAGCCATCCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCAt                                                                                 <  7:22002/149‑1 (MQ=255)                                                                    gggTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATc                                                                                <  6:219485/149‑1 (MQ=255)                                                                      gTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTa                                                                                                                                 >  4:45080/1‑98 (MQ=255)                                                                      gTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTa                                                                                                                                 <  3:45080/98‑1 (MQ=255)                                                                             aGGCAACCCAGATTTCCGGGCTCGACAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGAAGTgcg                                                                       >  3:306063/1‑149 (MQ=255)                                                                               gCAACCCAGATTTACGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGACCCGGTCAGGGCGCG‑TCGAGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGcc                                                                     >  8:21886/1‑149 (MQ=255)                                                                                cAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGcc                                                                     >  5:63393/1‑148 (MQ=255)                                                                                 aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGa                                                                   <  7:272310/149‑1 (MQ=255)                                                                                    ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATAAGCAGTGCGCCGAGCa                                                                >  3:237655/1‑149 (MQ=255)                                                                                             cGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGAGTACGTGATACGGAGGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGGGCGCCGAGCATTATGca                                                         >  4:48687/1‑147 (MQ=255)                                                                                              gggCTTGGCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAg                                                      <  4:351612/149‑1 (MQ=255)                                                                                                    gTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATAc                                                >  8:43091/1‑149 (MQ=255)                                                                                                           aTTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACt                                         >  2:224326/1‑149 (MQ=255)                                                                                                            ttCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTg                                        >  1:307974/1‑149 (MQ=255)                                                                                                              ccATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACt                                         <  1:305004/146‑1 (MQ=255)                                                                                                              ccATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACt                                         >  2:305004/1‑146 (MQ=255)                                                                                                                aTGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAtt                                                              >  1:143131/1‑123 (MQ=255)                                                                                                                aTGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCAtt                                                              <  2:143131/123‑1 (MQ=255)                                                                                                                aTGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGccc                                     >  6:47618/1‑148 (MQ=255)                                                                                                                aTGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCa                                    >  3:24498/1‑149 (MQ=255)                                                                                                                      tttCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAg                              >  5:179862/1‑149 (MQ=255)                                                                                                                        tCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGc                             <  4:342032/148‑1 (MQ=255)                                                                                                                           tt‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTg                          >  6:208041/1‑148 (MQ=255)                                                                                                                                  cGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCt                   <  5:208041/149‑1 (MQ=255)                                                                                                                                  cGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCt                   <  6:63393/149‑1 (MQ=255)                                                                                                                                   ggTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtg                  <  3:187548/149‑1 (MQ=255)                                                                                                                                       tacgcgaCGCGCCCTGACCGGGGCGCGTTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                 >  1:221727/7‑147 (MQ=255)                                                                                                                                       tacgcgaCGCGCCCTGACCGGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtga                 <  2:221727/140‑1 (MQ=255)                                                                                                                                          tcAACGCGCCCCGGTAAGGGCGCG‑TCGCGTACGTGATACGGATGTAATAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCg           >  6:343049/2‑149 (MQ=255)                                                                                                                                            gaCGCGCCCTGACCGGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGcc         >  2:376340/2‑149 (MQ=255)                                                                                                                                            aaCGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGcc         <  7:21886/149‑1 (MQ=255)                                                                                                                                              cgcgCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGa       <  8:182778/149‑1 (MQ=255)                                                                                                                                               gcgcCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGaa      <  2:307974/149‑1 (MQ=255)                                                                                                                                                   cccGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGACGCGCTGATGATCGCCGAACTCg  <  3:48687/149‑1 (MQ=255)                                                                                                                                                    ccGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGcc         <  3:190921/141‑1 (MQ=255)                                                                                                                                                    ccGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGcc         >  4:190921/1‑141 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ATGCGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCG  >  NC_002947/156236‑156529

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