Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	156,388	A→G	7.2%	intergenic (‑13/+30)	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,388	0	A	G	7.2%	199.1 / 4.4	69	intergenic (‑13/+30)	ygdQ/citN	transport protein/citrate transporter
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base G (2/3);  total (34/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.95e-01

CAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCG  >  NC_002947/156247‑156529                                                                                                                                               |                                                                                                                                               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                                                                                                                                   >  7:365943/1‑149 (MQ=255)      tCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcg‑tcgcg                                                                                                                                   >  2:249780/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                                                                                                                                <  8:364233/114‑1 (MQ=255)                                            tcAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTAc                                                                                                                                >  7:364233/1‑114 (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)                                                                                                                                               |                                                                                                                                               CAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTT‑CGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCG‑TCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTCG  >  NC_002947/156247‑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: ‑