Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,831,780	T→C	5.9%	L220P (CTG→CCG)	yfiP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,831,780	0	T	C	5.9%	194.4 / 2.1	67	L220P (CTG→CCG)	yfiP	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base C (2/2);  total (30/37)
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.41e-01

GCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCG  >  NC_002947/1831634‑1831928                                                                                                                                                    |                                                                                                                                                     gcggcgCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGc                                                                                                                                                    >  1:436509/1‑149 (MQ=255)   ggcgCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCg                                                                                                                                                  >  2:142100/1‑149 (MQ=255)      gcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCtt                                                                                                                                               <  5:409908/149‑1 (MQ=255)          cAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                                            <  2:59236/148‑1 (MQ=255)          cAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                                            <  2:52972/148‑1 (MQ=255)          cAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTgc                                                                                                                                             <  7:244390/147‑2 (MQ=255)           aaGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGag                                                                                                                                            >  4:400367/1‑146 (MQ=255)            aGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAg                                                                                                                                         >  8:303679/1‑149 (MQ=255)             gcgcGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTAGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGc                                                                                                                                        <  2:304603/149‑1 (MQ=255)                  atgatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgac                                                                                                                                   >  4:312653/1‑149 (MQ=255)                   tgatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgacg                                                                                                                                  >  7:336974/1‑149 (MQ=255)                    gatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgacga                                                                                                                                 >  5:29631/1‑149 (MQ=255)                     atCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgacgac                                                                                                                                <  3:195144/149‑1 (MQ=255)                        aTTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGacac                                                                                                                              <  2:436509/148‑1 (MQ=255)                            cTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTCAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcacccac                                                                                                                          <  3:10333/148‑1 (MQ=255)                                   gCAGAAGTTGCGGCTATGTTCCTGGACCTTTCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                                  <  8:336974/149‑1 (MQ=255)                                   gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                                  <  7:258867/149‑1 (MQ=255)                                   gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                                  >  2:253072/1‑149 (MQ=255)                                    cAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                                 >  5:78921/1‑149 (MQ=255)                                     aGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGc                                                                                                                <  4:76191/149‑1 (MQ=255)                                      gAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTGCCg                                                                                                                                           <  1:223492/121‑1 (MQ=255)                                      gAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTGCCg                                                                                                                                           >  2:223492/1‑121 (MQ=255)                                                     tgCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGt                                                                                                 <  8:298361/148‑1 (MQ=255)                                                        cTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTcc                                                                                               >  4:224794/1‑147 (MQ=255)                                                           gACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCTAAggg                                                                                          <  2:352058/149‑1 (MQ=255)                                                           gACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAggg                                                                                          <  2:357492/149‑1 (MQ=255)                                                           gACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAg                                                                                            <  2:142434/147‑1 (MQ=255)                                                              cTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGGCATGCATTACATGAGCGGCAAGCGCCGGATGCCGCTTGACGACCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgggg                                                                                         >  1:119559/1‑148 (MQ=255)                                                               ttGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGca                                                                                       >  5:60860/1‑148 (MQ=255)                                                                     gCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAAccc                                                                                 <  7:303679/148‑1 (MQ=255)                                                                       gATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTg                                                                              <  2:89303/149‑1 (MQ=255)                                                                       gATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTg                                                                              <  1:253072/149‑1 (MQ=255)                                                                        aTGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACCACCCCTGc                                                                             >  5:56872/1‑149 (MQ=255)                                                                         tGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCt                                                                            <  6:78921/149‑1 (MQ=255)                                                                           gCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCtt                                                                           <  6:253111/148‑1 (MQ=255)                                                                             cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTgg                                                                        <  1:323004/149‑1 (MQ=255)                                                                                  ttCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCa                                                                    >  2:150746/1‑148 (MQ=255)                                                                                     gCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATaa                                                                 <  3:224794/148‑1 (MQ=255)                                                                                       aGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGa                                                              <  2:119559/149‑1 (MQ=255)                                                                                         gcgcTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACg                                                            >  4:361421/1‑149 (MQ=255)                                                                                                aTGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTgg                                                     <  6:148606/149‑1 (MQ=255)                                                                                                    cTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGc                                                 <  3:400367/149‑1 (MQ=255)                                                                                                     tATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCta                                                                                                     <  5:146483/96‑1 (MQ=255)                                                                                                     tATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCta                                                                                                     >  6:146483/1‑96 (MQ=255)                                                                                                          ggATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGt                                           >  6:403356/1‑149 (MQ=255)                                                                                                           gATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGt                                           <  1:150746/148‑1 (MQ=255)                                                                                                           gATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTg                                          >  5:232501/1‑149 (MQ=255)                                                                                                                gtTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGc                                       >  6:200699/1‑147 (MQ=255)                                                                                                                 ttCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGc                                       >  1:171682/1‑146 (MQ=255)                                                                                                                       ccTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCCGGGTGGCCggccggga                               >  8:69292/1‑142 (MQ=255)                                                                                                                           cATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCTAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAgg                            <  7:215364/148‑1 (MQ=255)                                                                                                                           cATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAggg                           <  5:412932/149‑1 (MQ=255)                                                                                                                                 cTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGc                                                                                        <  5:217856/82‑1 (MQ=255)                                                                                                                                 cTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGc                                                                                        >  6:217856/1‑82 (MQ=255)                                                                                                                                  tGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTgg                    >  1:7430/1‑149 (MQ=255)                                                                                                                                  caaGCGGCAGCCGCCGCTTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTgg                    >  4:236238/3‑149 (MQ=255)                                                                                                                                    aGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGc                   <  5:403356/148‑1 (MQ=255)                                                                                                                                        tcAAGCGGCAGCCGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAt              >  6:120308/2‑149 (MQ=255)                                                                                                                                        gCAAGCGGCGGCTGCCGGTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAt              >  5:422417/1‑149 (MQ=255)                                                                                                                                        gCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATaa                                                                 >  5:217897/1‑98 (MQ=255)                                                                                                                                        gCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATaa                                                                 <  6:217897/98‑1 (MQ=255)                                                                                                                                        gCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAt              <  6:60860/149‑1 (MQ=255)                                                                                                                                          aaGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAtct            >  7:56703/1‑149 (MQ=255)                                                                                                                                            gcggcgGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCTAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTcc          <  5:374773/149‑1 (MQ=255)                                                                                                                                             cggcggCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTccc         <  2:7430/149‑1 (MQ=255)                                                                                                                                                   cTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGAcc   <  6:29631/149‑1 (MQ=255)                                                                                                                                                   cTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCATTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGCCGATTGGCTTGATCTCCCTTGAcc   >  4:301449/1‑149 (MQ=255)                                                                                                                                                    tGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCg  >  6:288634/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCA‑GCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCG  >  NC_002947/1831634‑1831928

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