Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,741,967	T→C	8.0%	intergenic (‑35/+378)	PP_3305 ← / ← PP_3306	TerC family membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,741,967	0	T	C	8.0%	180.7 / 5.1	65	intergenic (‑35/+378)	PP_3305/PP_3306	TerC family membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base C (2/3);  total (30/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 = 9.09e-01

GCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTAACTTCGCAC  >  NC_002947/3741821‑3742114                                                                                                                                                    |                                                                                                                                                    gCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCATAGTTGTTCAACTACGa                                                                                                                                                   <  3:179842/149‑1 (MQ=255)     acCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTcc                                                                                                                                               <  4:127001/149‑1 (MQ=255)          aGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGAc                                                                                                                                            <  8:287572/147‑1 (MQ=255)              aaGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTTGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCg                                                                                                                                      >  2:219919/1‑149 (MQ=255)              aaGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGc                                                                                                                                       <  6:394681/148‑1 (MQ=255)              aaGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCg                                                                                                                                      >  7:151408/1‑149 (MQ=255)               aGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGCTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGATCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTAAGACATGGCgg                                                                                                                                     >  5:407099/1‑149 (MQ=255)                 cgcCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGAGGta                                                                                                                                   >  8:293577/1‑149 (MQ=255)                                cacaACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAgggg                                                                                                                    >  7:304253/1‑149 (MQ=255)                                    acaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCg                                                                                                                <  2:105674/149‑1 (MQ=255)                                    acaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCg                                                                                                                >  2:165429/1‑149 (MQ=255)                                         gAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGt                                                                                                            <  8:151408/148‑1 (MQ=255)                                              aGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGaa                                                                                                      <  2:406440/149‑1 (MQ=255)                                                  aagccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATaaa                                                                                                  <  3:305646/149‑1 (MQ=255)                                                   agccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTc                                                                                                           >  4:406854/1‑139 (MQ=255)                                                   agccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTc                                                                                                           <  3:406854/139‑1 (MQ=255)                                                     ccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGaa                                                                                               <  8:304253/149‑1 (MQ=255)                                                      caaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACAGTCAGAGGGTCCCGGCGTCACGAATAAAGAAt                                                                                              >  6:257411/1‑149 (MQ=255)                                                       aaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTc                                                                                                           <  2:235421/135‑1 (MQ=255)                                                       aaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTc                                                                                                           >  1:235421/1‑135 (MQ=255)                                                       aaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCATAGTTGTTCAACTACGACtatac                                                                                                                                             <  6:195662/101‑5 (MQ=38)                                                       aaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCATAGTTGTTCAACTACGACtatac                                                                                                                                             >  5:195662/1‑97 (MQ=38)                                                            caGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCTACATGGCGGTATAACCGTCAGAGGGGCCCGGAGTCACGAATAAACAATTTATc                                                                                         >  5:349756/1‑148 (MQ=255)                                                                    aGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAg                                                                                >  7:164144/1‑149 (MQ=255)                                                                          aaGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCt                                                                          >  8:112310/1‑149 (MQ=255)                                                                            gAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTg                                                                         >  6:176768/1‑148 (MQ=255)                                                                            gAACGTCTTCGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGc                                                                        <  1:219919/149‑1 (MQ=255)                                                                                 tCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                                   <  1:165429/149‑1 (MQ=255)                                                                                 tCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                                   >  8:274642/1‑149 (MQ=255)                                                                                   ttGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGa                                                                 <  7:112310/149‑1 (MQ=255)                                                                                   ttGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTACGACATGGCGGTATAACCGTCAGAGGGGCACGGCGTAACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAg                                                                  >  7:337445/1‑148 (MQ=255)                                                                                   ttGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCATAGTTGTTCAACTACGACt                                                                                                                                                 <  1:237261/69‑1 (MQ=37)                                                                                   ttGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCATAGTTGTTCAACTACGACt                                                                                                                                                 >  2:237261/1‑69 (MQ=37)                                                                                       tAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAAc                                                             <  5:176768/149‑1 (MQ=255)                                                                                          agaTACGTACTAAGGGCAGACATGAGCCCCTCCTAGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCagag                                                                                                                       >  1:263084/1‑89 (MQ=255)                                                                                          agaTACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCagag                                                                                                                       <  2:263084/88‑1 (MQ=255)                                                                                                gTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACaaa                                                    >  7:37547/1‑149 (MQ=255)                                                                                                  aCTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATc                                                  >  8:30133/1‑149 (MQ=255)                                                                                                   cTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGAGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATc                                                  <  7:293577/148‑1 (MQ=255)                                                                                                      aGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACTAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGt                                              <  5:257411/149‑1 (MQ=255)                                                                                                          cAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTAAGACTCCGACATGGCTGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                                          <  3:186/149‑1 (MQ=255)                                                                                                           cgagaTGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                                          <  5:361861/144‑1 (MQ=255)                                                                                                             aCATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATc                                                  <  4:82884/138‑1 (MQ=255)                                                                                                             aCATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATc                                                  >  3:82884/1‑138 (MQ=255)                                                                                                                 gAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGa                                                                 >  4:336981/1‑119 (MQ=255)                                                                                                                 gAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGa                                                                 <  3:336981/119‑1 (MQ=255)                                                                                                                    ccccTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTa                                 >  7:219269/1‑148 (MQ=255)                                                                                                                          ct‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCt                          >  3:146506/1‑149 (MQ=255)                                                                                                                              tataGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTg                        <  2:406078/148‑1 (MQ=255)                                                                                                                              tataGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTg                        <  8:164144/148‑1 (MQ=255)                                                                                                                                taGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCAAGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAg                     <  3:49283/149‑1 (MQ=255)                                                                                                                                   tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTg                   >  5:11229/1‑148 (MQ=255)                                                                                                                                        gTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGaaaaa              <  7:274642/148‑1 (MQ=255)                                                                                                                                           gAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAaacta           >  7:257661/1‑148 (MQ=255)                                                                                                                                            aacaacTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAAt                                                   >  4:184247/1‑107 (MQ=255)                                                                                                                                            aacaacTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAAt                                                   <  3:184247/107‑1 (MQ=255)                                                                                                                                            aacaacTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAaactaac         <  3:304606/149‑1 (MQ=255)                                                                                                                                            aaacacTATGACTCCGACATGGCGGTATATCCGTCAGAGGGGCCCGTCGTCACTAATAAAGAATTTATCGAGCGCTAGTAGGCTGCAAGTAGAAAAGATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAaactaa          <  6:407099/144‑1 (MQ=255)                                                                                                                                                acTATTACTCCTTCATGGCGTTATAACCGTCAGAGGGGCCCGGCGTCACGATTAAATAATTTATCGAGCGCGAGTAGGGTGCAATTATAAAACATTTACAAATCTTTTCCATGACAAGCTAAAACGCTTGAAGTGAAAAACTAACttgg     <  6:349756/149‑3 (MQ=255)                                                                                                                                                 cTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACaa                                     >  1:285432/1‑116 (MQ=255)                                                                                                                                                 cTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACaa                                     <  2:285432/116‑1 (MQ=255)                                                                                                                                                 cTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCt                                  <  7:140797/119‑1 (MQ=255)                                                                                                                                                 cTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCt                                  >  8:140797/1‑119 (MQ=255)                                                                                                                                                  tATGACTCCGTCATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTAACTTCGca   <  7:381619/149‑1 (MQ=255)                                                                                                                                                   aTGACTCCGCCATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTAACTTCGcac  >  4:261545/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCT‑GTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTAACTTCGCAC  >  NC_002947/3741821‑3742114

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