Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,955,965	T→A	6.1%	A201A (GCT→GCA)	PP_1753 →	DMT superfamily permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,955,965	0	T	A	6.1%	234.9 / 2.8	77	A201A (GCT→GCA)	PP_1753	DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (34/38);  new base A (2/3);  total (36/41)
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.45e-01

TGCTGACTTCATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGTCGCGGGTA  >  NC_002947/1955817‑1956112                                                                                                                                                       |                                                                                                                                                    tgctgACTTCATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGCCCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAgctgcgct                                                                                                                                                     >  6:348822/1‑149 (MQ=255) tgctgACTTCATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAgctgcgct                                                                                                                                                     >  6:261532/1‑149 (MQ=255)  gctgACTTCATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAgctgcgctg                                                                                                                                                    <  3:54931/149‑1 (MQ=255)        ttCATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAgctgcgctgcgctgc                                                                                                                                              >  6:190740/1‑149 (MQ=255)          cATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCgctgct                                                                                                                                             <  5:82764/148‑1 (MQ=255)          cATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCgctgct                                                                                                                                             <  3:182272/148‑1 (MQ=255)          cATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCgctgct                                                                                                                                             <  2:268447/148‑1 (MQ=255)          cATTGGGGACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCgctgct                                                                                                                                             <  2:87301/148‑1 (MQ=255)           aTTGGGTACCCACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCAGCGCAGCGCAGCTGa                                                                                                                                           <  2:317624/149‑1 (MQ=255)              gggTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGTTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCAGCGCAGCGCAGCTGagg                                                                                                                                         <  2:117508/148‑3 (MQ=255)              gggTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCAGCGCAGCGCAGCTGagg                                                                                                                                         <  2:56044/148‑3 (MQ=255)                 tACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCtt                                                                                                                                     <  6:103856/149‑1 (MQ=255)                 tACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCtt                                                                                                                                     <  1:281964/149‑1 (MQ=255)                   ccTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTgg                                                                                                                                   <  6:312197/149‑1 (MQ=255)                     tACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGt                                                                                                                                 >  6:131140/1‑149 (MQ=255)                      acacCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTg                                                                                                                                >  6:362799/1‑149 (MQ=255)                      acacCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGACCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTggg                                                                                                                                  >  1:357860/1‑147 (MQ=255)                      acacCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGACCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTggg                                                                                                                                  <  2:357860/147‑1 (MQ=255)                       caccacCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGc                                                                                                                               <  6:361912/149‑1 (MQ=255)                             cGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGc                                                                                                                         <  5:190740/149‑1 (MQ=255)                               gCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTg                                                                                                                       >  8:75178/1‑149 (MQ=255)                                 gtgtTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGa                                                                                                                     <  5:261532/149‑1 (MQ=255)                                    ttGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGc                                                                                                                  >  6:226117/1‑149 (MQ=255)                                     tGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGcc                                                                                                                 >  7:281539/1‑149 (MQ=255)                                      gACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACTCCCTGAGCCAGAAGCAGTTGCTGACCGTGTTGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAgctgcgctgcgctgc                                                                                                                                              <  1:258961/119‑1 (MQ=255)                                      gACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAgctgcgctgcgctgc                                                                                                                                              >  2:258961/1‑119 (MQ=255)                                      gACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCt                                                                                                                >  7:12574/1‑149 (MQ=255)                                       aCCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTTATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCtt                                                                                                               <  8:58272/149‑1 (MQ=255)                                         cATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGAACGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTg                                                                                                              >  3:68981/1‑148 (MQ=255)                                             ttGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACt                                                                                                         >  1:354518/1‑149 (MQ=255)                                               gCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGa                                                                                                       >  5:370375/1‑149 (MQ=255)                                                     gACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGCTCTACCTCAGCTGCGCTGCGCTGCTGCCACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACACCgg                                                                                                 >  8:164425/1‑149 (MQ=255)                                                        cAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGc                                                                                              <  8:180886/149‑1 (MQ=255)                                                            aTCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCAGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAgg                                                                                          >  5:181563/1‑149 (MQ=255)                                                             tCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAggc                                                                                         >  6:84899/1‑149 (MQ=255)                                                               tGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGCGACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAggct                                                                                        >  1:86093/1‑149 (MQ=255)                                                                 ggTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGt                                                                                                >  8:225103/1‑138 (MQ=255)                                                                 ggTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGt                                                                                                <  7:225103/138‑1 (MQ=255)                                                                       cTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGgctgctg                                                                                <  2:354518/148‑1 (MQ=255)                                                                        tACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAGCCAGGTAATGATGGTGAACTACATCAGATGCGATGCGATGCTGACACCTTGGTTGCACCCGCTGGAGGCCTTGCAAATGACCCCGGTGCAAGGCTGgctgctgc                                                                               >  1:278250/1‑149 (MQ=255)                                                                        tACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGAc                                                                                                      <  6:126896/125‑1 (MQ=255)                                                                        tACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGAc                                                                                                      >  5:126896/1‑125 (MQ=255)                                                                        tACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGgctgctgct                                                                              >  1:68723/1‑149 (MQ=255)                                                                             ccTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCt                                                                         <  5:131140/149‑1 (MQ=255)                                                                             ccTGAGCCAGAAGCAGTTGCTGAACGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCt                                                                         <  4:68981/149‑1 (MQ=255)                                                                                    cAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCt                                                                  >  1:88572/1‑149 (MQ=255)                                                                                     aGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTCCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTg                                                                 <  1:276697/149‑1 (MQ=255)                                                                                      gAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTg                                                                 >  3:99655/1‑148 (MQ=255)                                                                                      gAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTg                                                                 >  7:336320/1‑148 (MQ=255)                                                                                      gAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTg                                                                 >  7:79665/1‑148 (MQ=255)                                                                                      gAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGa                                                                <  6:370375/149‑1 (MQ=255)                                                                                       aaGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGaa                                                               <  7:75178/149‑1 (MQ=255)                                                                                       aaGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGaa                                                               >  1:192022/1‑149 (MQ=255)                                                                                           aGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCt                                                                                                                <  1:101630/96‑1 (MQ=255)                                                                                           aGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCt                                                                                                                >  2:101630/1‑96 (MQ=255)                                                                                           aGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACAcc                                                           <  5:84899/149‑1 (MQ=255)                                                                                            gTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACAccc                                                          >  2:143771/1‑149 (MQ=255)                                                                                            gTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACAccc                                                          <  2:192298/149‑1 (MQ=255)                                                                                             ttGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCt                                                         >  6:62004/1‑149 (MQ=255)                                                                                                cTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCtggt                                                      <  5:226117/149‑1 (MQ=255)                                                                                                  gACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCtggtgg                                                    >  8:53600/1‑149 (MQ=255)                                                                                                     cGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCtt                                                 <  2:86093/149‑1 (MQ=255)                                                                                                     cGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATATACCTCAGCTGCCCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTCCTGCTGGCCTGTTGTCTGAACACCCAGGTGGCt                                                  >  1:175540/1‑148 (MQ=255)                                                                                                           gC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGgcg                                           >  3:356112/1‑149 (MQ=255)                                                                                                            c‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGgcgc                                          <  8:281539/149‑1 (MQ=255)                                                                                                                  gcag‑caGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTc                                      <  4:99655/148‑1 (MQ=255)                                                                                                                   cag‑cagGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTTCAAGTCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGc                                    <  2:43308/149‑1 (MQ=255)                                                                                                                    ag‑cagGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGcc                                   >  5:355690/1‑149 (MQ=255)                                                                                                                       cagGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCTACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCg                                  >  5:38168/1‑148 (MQ=255)                                                                                                                        agGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGaa                                <  2:192022/149‑1 (MQ=255)                                                                                                                        agGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGaa                                <  2:88572/149‑1 (MQ=255)                                                                                                                          gTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGc                                                                           >  1:379971/1‑104 (MQ=255)                                                                                                                          gTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGc                                                                           <  2:379971/104‑1 (MQ=255)                                                                                                                             atgatgGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCt                           <  6:355690/149‑1 (MQ=255)                                                                                                                              tgatgGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTg                          <  2:314275/149‑1 (MQ=255)                                                                                                                              tgatgGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGC‑TTCGGCGaacgcgccg                         <  7:65592/149‑8 (MQ=255)                                                                                                                               gatgGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTgg                         >  5:168462/1‑149 (MQ=255)                                                                                                                                atgGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGgc                        >  8:335382/1‑149 (MQ=255)                                                                                                                                        cTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTgg                 >  4:334402/1‑148 (MQ=255)                                                                                                                                         ggtgtCAGCAGCGCAGCGCAGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTgg                 >  7:308054/5‑147 (MQ=255)                                                                                                                                           ccTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTg                  <  6:168462/144‑1 (MQ=255)                                                                                                                                           ccTCAGCTGCGCTGCGCTGATGACAGATTTGGTGCACCCGCTGGAGGACTTTCAACTGACCCCTGAGCAAGGCTGGCTGCTGATGGCCTGTTGTCTGAACACCATTGTGTCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGa               <  4:281938/147‑1 (MQ=255)                                                                                                                                               aGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGt          <  4:356112/148‑1 (MQ=255)                                                                                                                                                  tgcgctgcgctgcTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGTcgcg      >  6:245358/1‑149 (MQ=255)                                                                                                                                                      ctgcgctgcTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGTCGCGGGTa  <  5:362799/149‑1 (MQ=255)                                                                                                                                                       tgcgctgcTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGTCGCGGGt   >  2:17042/1‑147 (MQ=255)                                                                                                                                                       |                                                                                                                                                    TGCTGACTTCATTGGGTACCTACACCACCGGCGTGTTGACCATCTTGCTGGCGACCAGCATCTGGGTGTTCTACGCCCTGAGCCAGAAGCAGTTGCTGACCGTGTGGC‑ACTCGCAG‑CAGGTAATGATGGTGATCTACCTCAGCTGCGCTGCGCTGCTGACACCTTGGGTGCACCCGCTGGAGGCCTTGCAACTGACCCCGGTGCAAGGCTGGCTGCTGCTGGCCTGTTGTCTGAACACCCTGGTGGCTTACGGCGCGTTCGCCGAAGCGCTGGCGCACTGGGAGGCGTCGCGGGTA  >  NC_002947/1955817‑1956112

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