Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,002,446	T→A	5.9%	V796D (GTC→GAC)	PP_1786 →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,002,446	0	T	A	5.9%	204.5 / 2.7	68	V796D (GTC→GAC)	PP_1786	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (35/29);  new base A (2/2);  total (37/31)
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.55e-01

AACGTATTTTTCCCACCGCAGGCGATTGGTGGTGCAACGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGATCCACTTCCATTGCG  >  NC_002947/2002298‑2002592                                                                                                                                                     |                                                                                                                                                   aaCGTATTTTTCCCACCGCAGGCGATTGGTGGTGCAACGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGt                                                                                                                                                    >  4:244571/1‑149 (MQ=255) aaCGTATTTTTCCCACCGCAGGCGATTGGTGGTGCAACGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGt                                                                                                                                                    >  8:211491/1‑149 (MQ=255)            cccACCTCAGGCGATTGGTGGTGCAACGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTAt                                                                                                                                          <  3:208368/148‑1 (MQ=255)                                gtgCAACGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAg                                                                                                                                       <  4:272046/131‑1 (MQ=255)                                gtgCAACGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAg                                                                                                                                       >  3:272046/1‑131 (MQ=255)                                 tgCAACGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGggtctacaggt                                                                                                                    <  2:73943/149‑1 (MQ=255)                                    aaCGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGggtctacaggtcta                                                                                                                 <  4:13606/149‑1 (MQ=255)                                     aCGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTg                                                                                                                                  <  7:161233/131‑1 (MQ=255)                                     aCGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTg                                                                                                                                  >  8:161233/1‑131 (MQ=255)                                       gcgTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGGGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACaa                                                                                                              <  7:147670/149‑1 (MQ=255)                                         gtgtCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAAcc                                                                                                            >  1:229012/1‑149 (MQ=255)                                              gTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGAGGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCtgt                                                                                                        <  4:198194/148‑1 (MQ=255)                                                  cTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGt                                                                                                   <  2:93887/149‑1 (MQ=255)                                                   tGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCtgctgctga                                                                                                  >  5:169499/1‑142 (MQ=255)                                                   tGATAACTTATCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGATGATGTTGAGTGCAGAGAAGCTCATCAGATGACTGTCTACACCTATGAGGGATGCCCGGTCTACAGGTCTACAACACTGTGGCGt                                                                                                   >  5:16972/1‑148 (MQ=255)                                                     aTAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGa                                                                                                 <  7:211491/148‑1 (MQ=255)                                                      tAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCAGGGTCTACAGGTCTACAACCCTGTGGCGTGAg                                                                                                >  1:13311/1‑148 (MQ=255)                                                        aCTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGca                                                                                              <  1:360663/148‑1 (MQ=255)                                                                gagTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCGTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCGATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGAGGCGTGAGCACATGGAGTg                                                                                     >  8:127765/1‑149 (MQ=255)                                                                 agTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTgg                                                                                    >  4:82783/1‑149 (MQ=255)                                                                        aaGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGACTACAcc                                                                                                                                             >  8:186361/1‑85 (MQ=255)                                                                        aaGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGACTACAcc                                                                                                                                             <  7:186361/85‑1 (MQ=255)                                                                        aaGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGACTACAcc                                                                                                                                             >  8:191601/1‑85 (MQ=255)                                                                        aaGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGACTACAcc                                                                                                                                             <  7:191601/85‑1 (MQ=255)                                                                              ttACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGAcc                                                                       >  3:5634/1‑149 (MQ=255)                                                                                aCGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCt                                                                                                          <  2:280101/112‑1 (MQ=255)                                                                                aCGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCt                                                                                                          >  1:280101/1‑112 (MQ=255)                                                                                aCGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAg                                                                                                >  7:74715/1‑122 (MQ=255)                                                                                aCGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAg                                                                                                <  8:74715/122‑1 (MQ=255)                                                                                   ccGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCa                                                                     >  2:45474/1‑146 (MQ=255)                                                                                    ctaGGACCTCGATGTCTGTGTCTTCACTGCTGATTTTGAGTTCAGAGCCGCTCATCAGATTACTGTCTACAACTATGAGGGCTGCCGGGTCTACATGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAgaga                                                                 <  1:298371/147‑1 (MQ=255)                                                                                           cTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATcc                                                                               <  7:61942/128‑1 (MQ=255)                                                                                           cTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATcc                                                                               >  8:61942/1‑128 (MQ=255)                                                                                           cTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGCCCCAGAGATGTCTCg                                                          >  8:178035/1‑149 (MQ=255)                                                                                           cTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCg                                                          >  2:80098/1‑149 (MQ=255)                                                                                           cTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCg                                                          >  7:185493/1‑149 (MQ=255)                                                                                           cTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCg                                                          >  6:108978/1‑149 (MQ=255)                                                                                             cGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGcc                                                        >  5:7969/1‑149 (MQ=255)                                                                                                  tCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGtt                                                    <  1:45474/148‑1 (MQ=255)                                                                                                   cTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACGTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATc                                                                                >  4:67778/1‑119 (MQ=255)                                                                                                   cTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATc                                                                                <  3:67778/119‑1 (MQ=255)                                                                                                    tgtgTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTcc                                                  >  1:210437/1‑148 (MQ=255)                                                                                                    tgtgTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTcc                                                  >  1:351564/1‑148 (MQ=255)                                                                                                     gtgtCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTcc                                                  >  8:254284/1‑147 (MQ=255)                                                                                                      tgtCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCaa                                                <  8:185493/148‑1 (MQ=255)                                                                                                         cTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGt                                            >  8:153079/1‑149 (MQ=255)                                                                                                             aCTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCAGATGGATTGGCATCCTAAAGACCCAGAGATGTATAGCCTGTTCCAAGAGTTCCt                                        <  6:16972/149‑1 (MQ=255)                                                                                                              ctgctgATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTg                                       <  6:7969/149‑1 (MQ=255)                                                                                                              ctgctgATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTg                                       <  2:13311/149‑1 (MQ=255)                                                                                                                 ctgATGTTGAGTGCAGAGCCGCTCATAAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATACTAAAGACCCAGAGATGTATCGCCTGTTCAAAGAGTTCCTGGAg                                    >  5:358464/1‑149 (MQ=255)                                                                                                                  tgATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCAT‑CTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCt                                  >  2:110335/1‑149 (MQ=255)                                                                                                                          aGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGaaaaa                           >  3:69633/1‑149 (MQ=255)                                                                                                                          aGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGaaaaa                           >  3:74415/1‑149 (MQ=255)                                                                                                                            tGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAAACCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAAc                          >  1:309601/1‑148 (MQ=255)                                                                                                                                  gCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGAttt                   <  7:178035/149‑1 (MQ=255)                                                                                                                                    cGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGa                 <  3:244571/149‑1 (MQ=255)                                                                                                                                     gCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTgg                                      >  4:155437/1‑127 (MQ=255)                                                                                                                                     gCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTgg                                      <  3:155437/127‑1 (MQ=255)                                                                                                                                      cTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGATc               <  4:69633/149‑1 (MQ=255)                                                                                                                                      cTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGATc               <  4:74415/149‑1 (MQ=255)                                                                                                                                       tcatcaGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGa                 <  7:316934/146‑1 (MQ=255)                                                                                                                                       tcatcaGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGa                 >  8:316934/1‑146 (MQ=255)                                                                                                                                         atcaGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTTCAACGCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTACTGGAGCTGGAAAAACCCGATTTGATCCAc            >  1:274675/1‑149 (MQ=255)                                                                                                                                         atcaGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAGAACCCGATTTGATCCAc            >  4:217828/1‑149 (MQ=255)                                                                                                                                         atcaGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGATcc              <  5:108978/147‑1 (MQ=255)                                                                                                                                               tGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGATCCACTTCCa       <  2:229012/148‑1 (MQ=255)                                                                                                                                                gACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCaa                                                <  8:60393/106‑1 (MQ=255)                                                                                                                                                gACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCaa                                                >  7:60393/1‑106 (MQ=255)                                                                                                                                                gACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGATCCACTTCCAtt     >  6:86942/1‑149 (MQ=255)                                                                                                                                                  cTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGATCCACTTCCATTGc   >  1:342700/1‑149 (MQ=255)                                                                                                                                                     tCTACACCTATGAGGGCTGCCGTGTCTACAGGTCTACAACGCTGTGTCGTGAGCACATGGATTGGCATCCTAAAGAAACAGAGATGTATCGCCTGTTCCAAGATTTCCTGGAGCTGGAAAAACCCGATTTGATCCACTTCCATTTCg  <  6:358464/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AACGTATTTTTCCCACCGCAGGCGATTGGTGGTGCAACGCGTGTCGTGGCTGATAACTTCTCAGAGTTCAGAAGCGCTTACGCCGAGGACCTCGATGTCTGTGTCTTCACTGCTGATGTTGAGTGCAGAGCCGCTCATCAGATGACTGTCTACACCTATGAGGGCTGCCGGGTCTACAGGTCTACAACCCTGTGGCGTGAGCACATGGATTGGCATCCTAAAGACCCAGAGATGTATCGCCTGTTCCAAGAGTTCCTGGAGCTGGAAAAACCCGATTTGATCCACTTCCATTGCG  >  NC_002947/2002298‑2002592

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