Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,444,691	Δ1 bp	5.6%	coding (141/642 nt)	nicF ←	maleamate amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,444,691	0	G	.	5.6%	235.2 / 6.9	72	coding (141/642 nt)	nicF	maleamate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base . (2/2);  total (36/36)
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.85e-01

TCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGT  >  NC_002947/4444545‑4444837                                                                                                                                                    |                                                                                                                                                   tCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGc                                                                                                                                                   <  7:206069/148‑1 (MQ=255)   cTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGCGTGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGccc                                                                                                                                                 >  1:304411/1‑149 (MQ=255)   cTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGcccc                                                                                                                                                >  4:52200/1‑148 (MQ=255)   cTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGcccc                                                                                                                                                >  3:225386/1‑148 (MQ=255)   cTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGccc                                                                                                                                                 <  8:27157/148‑1 (MQ=255)    ttCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCgg                                                                                                                                               >  5:337783/1‑149 (MQ=255)    ttCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCgg                                                                                                                                               >  7:127667/1‑149 (MQ=255)       atcaCCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCgggg                                                                                                                                             <  6:237031/148‑1 (MQ=255)          aCCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCCGGGCgca                                                                                                                                         <  4:225386/149‑3 (MQ=255)                gcCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTTCCGCGCCCACGCCCGGGGCGTAGAGAgg                                                                                                                                   >  8:317416/1‑149 (MQ=255)                 ccTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGt                                                                                                                                  >  8:184354/1‑149 (MQ=255)                   ttGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCAC‑CCCggg                                                                                                                                              >  5:360419/1‑134 (MQ=255)                   ttGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCAC‑CCCggg                                                                                                                                              <  8:118788/134‑1 (MQ=255)                   ttGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCAC‑CCCggg                                                                                                                                              <  6:360419/134‑1 (MQ=255)                   ttGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCAC‑CCCggg                                                                                                                                              >  7:118788/1‑134 (MQ=255)                        cACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑GGGTAGCGGATGGTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCg                                                                                                                           <  4:304818/149‑1 (MQ=255)                         aCCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCgg                                                                                                                          >  5:85076/1‑149 (MQ=255)                         aCCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCgg                                                                                                                          <  2:304411/149‑1 (MQ=255)                           ccACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCgggg                                                                                                                        >  2:47683/1‑149 (MQ=255)                               acGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCCGGGCGCACCGCTGTACGCCCcgggcgtg                                                                                                                    <  3:52200/149‑5 (MQ=255)                                  ccgccgTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGgtggtg                                                                                                                   >  5:174205/1‑147 (MQ=255)                                          aGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGc                                                                                                          <  6:85076/148‑1 (MQ=255)                                          aGCAAAATGCGGAGGC‑TGGTAGCGCATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCa                                                                                                         >  7:337214/1‑149 (MQ=255)                                              aaaTGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTCCATg                                                                                                       >  8:158222/1‑147 (MQ=255)                                                 tGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTc                                                                                                  <  1:47683/149‑1 (MQ=255)                                                    ggAGGC‑TGGTAGCGGATGTTGGTTTGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGGTCTACTGCCGCGACCACGCACGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGat                                                                                               <  7:184354/149‑1 (MQ=255)                                                     gAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatc                                                                                              <  6:203255/149‑1 (MQ=255)                                                        gC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGatcatc                                                                                           <  3:49061/149‑1 (MQ=255)                                                         c‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGTGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGAtcatca                                                                                          <  5:116484/149‑1 (MQ=255)                                                           tGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCAt                                                                                                        >  2:191166/1‑134 (MQ=255)                                                           tGGTAGCGGATCTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCAt                                                                                                        <  1:191166/134‑1 (MQ=255)                                                             gTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATcagca                                                                                       <  1:296603/149‑1 (MQ=255)                                                                gCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCC‑GGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGc                                                                                   >  4:30836/1‑149 (MQ=255)                                                                      gTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCt                                                                              >  5:310370/1‑149 (MQ=255)                                                                          gTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGc                                                                           <  2:12674/148‑1 (MQ=255)                                                                          gTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGc                                                                           <  3:191845/148‑1 (MQ=255)                                                                                accaccAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCa                                                                      >  7:57109/1‑147 (MQ=255)                                                                                 ccaccaAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGat                                                                                               >  7:134788/1‑121 (MQ=255)                                                                                 ccaccaAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGat                                                                                               <  8:134788/121‑1 (MQ=255)                                                                                  caccaAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCAGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGAGCGCCCGGGGTGGTGTAGCCCTGCAGGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGc                                                                  >  1:368638/1‑149 (MQ=255)                                                                                     caAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGa                                                               >  7:89036/1‑149 (MQ=255)                                                                                            ccGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCt                                                         <  2:318911/148‑1 (MQ=255)                                                                                             cGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTgg                                                       >  3:250829/1‑149 (MQ=255)                                                                                              gCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCaaa                                                                    >  5:163085/1‑135 (MQ=255)                                                                                              gCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCaaa                                                                    <  6:163085/135‑1 (MQ=255)                                                                                                aGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGccc                                                    <  7:69016/149‑1 (MQ=255)                                                                                                    gcgGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCcaca                                                >  6:64220/1‑149 (MQ=255)                                                                                                    gcgGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCcaca                                                >  4:370938/1‑149 (MQ=255)                                                                                                      gggCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGcccgcgact                                              >  5:101841/1‑143 (MQ=255)                                                                                                       ggCCAGCGCCAGCAGGCGGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCacacac                                             <  8:337214/149‑1 (MQ=255)                                                                                                        gccagcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACAcc                                            <  8:270231/149‑1 (MQ=255)                                                                                                            gcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTgg                                        <  6:310370/149‑1 (MQ=255)                                                                                                            gcgccagcAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTgg                                        <  4:250829/149‑1 (MQ=255)                                                                                                             cgccagcAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGt                                       <  6:337783/149‑1 (MQ=255)                                                                                                                cagcagGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGt                                    >  6:333279/1‑149 (MQ=255)                                                                                                                 agcagGCCGCCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGtt                                   >  8:184851/1‑149 (MQ=255)                                                                                                                               cTGCTCGACTGCCGCGACCAGGCCCGGGGCGTACAGCAGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCttt                     <  6:166832/149‑1 (MQ=255)                                                                                                                                 gCTCGACTGCCGCGACCACGCCCGGGGCGTACCGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGGCGATCCTCAGCAGGGCCGGAGGGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGCTAGATCTCCCGGGAGCatt                     >  6:312343/1‑147 (MQ=255)                                                                                                                                  cTCGTCTGCCGCGACCACGCCCGGGGCGTACAGCGGGGCGCGCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTtg                    <  8:238382/147‑1 (MQ=255)                                                                                                                                  cTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTtgtg                  <  7:158222/149‑1 (MQ=255)                                                                                                                                  cTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTtgt                   <  3:370938/148‑1 (MQ=255)                                                                                                                                  cTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTtgt                   <  5:64220/148‑1 (MQ=255)                                                                                                                                       cTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCC‑GGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGc            <  3:30836/149‑1 (MQ=255)                                                                                                                                          ccGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCt           >  8:316991/1‑148 (MQ=255)                                                                                                                                          ccGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTc          <  7:184851/149‑1 (MQ=255)                                                                                                                                           cgcgACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCt                                          <  5:75303/116‑1 (MQ=255)                                                                                                                                           cgcgACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCt                                          >  6:75303/1‑116 (MQ=255)                                                                                                                                           cgcgACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCa         >  6:159306/1‑149 (MQ=255)                                                                                                                                            gcgACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCAt        >  6:133297/1‑149 (MQ=255)                                                                                                                                                 cACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCg   >  3:37729/1‑149 (MQ=255)                                                                                                                                                 cACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAACAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCg   >  6:374141/1‑149 (MQ=255)                                                                                                                                                  aCGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGt  >  5:22922/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TCCTTCATCACCGGCGCCTTGCGCACCCACACGCCGCCGTCAGCAAAATGCGGAGGC‑TGGTAGCGGATGTTGGTATGCACCACCAAGGTACCGCAGTCGCGGGCCAGCGCCAGCAGGCCGGCGGCCTGCTCGACTGCCGCGACCACGCCCGGGGCGTACAGCGGTGCGCCCGGGGTGGTGTAGCCCTGCATGAAGTCGATCATCAGCAGGGCCGGCTTGCGGCCAAAGCCGATGCGCTGGCCCCACACACCCTGGTAGTTGTCACGGGCGCTTTGTGCGTCGCTCATGTCCGT  >  NC_002947/4444545‑4444837

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑