Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,156,605	T→G	7.1%	F31L (TTT→TTG)	yihS →	aldose‑ketose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,156,605	0	T	G	7.1%	213.7 / 4.9	70	F31L (TTT→TTG)	yihS	aldose‑ketose isomerase
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base G (2/3);  total (32/38)
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.47e-01

CGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTAGC  >  NC_002947/1156457‑1156749                                                                                                                                                      |                                                                                                                                                   cGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCttt                                                                                                                                                    <  3:246942/149‑1 (MQ=255)   ccATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTg                                                                                                                                                   >  3:5925/1‑148 (MQ=255)   ccATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTg                                                                                                                                                   <  4:5925/148‑1 (MQ=255)          cccGCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCAGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGcc                                                                                                                                           >  2:22124/1‑149 (MQ=255)             gCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagc                                                                                                                                        <  3:143226/149‑4 (MQ=255)                   tAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagcag                                                                                                                                      >  1:244008/1‑140 (MQ=255)                   tAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagcag                                                                                                                                      <  2:244008/145‑6 (MQ=255)                   tAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagca                                                                                                                                       >  1:184351/1‑140 (MQ=255)                   tAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagca                                                                                                                                       <  2:184351/144‑5 (MQ=255)                     gAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAc‑c                                                                                                                               <  7:194726/149‑1 (MQ=255)                       aTTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGCAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAc‑cc                                                                                                                              >  4:40114/1‑148 (MQ=255)                           gCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACgg                                                                                                                         >  3:83111/1‑149 (MQ=255)                              aaCGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGttt                                                                                                                      <  1:343665/149‑1 (MQ=255)                                cGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCg                                                                                                                    >  6:80116/1‑149 (MQ=255)                                 gACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTTGCCTAAGCACCAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCg                                                                                                                    <  1:86479/149‑1 (MQ=255)                                        aaGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCttt                                                                                                                                                    >  8:76821/1‑110 (MQ=255)                                        aaGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCttt                                                                                                                                                    <  7:76821/110‑1 (MQ=255)                                        aaGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTAGGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTg                                                                                                            >  4:195591/1‑149 (MQ=255)                                           gATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGAt                                                                                                         <  3:170367/149‑1 (MQ=255)                                                   tGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATgg                                                                                                 >  3:156465/1‑149 (MQ=255)                                                       ccATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCaa                                                                                             >  7:29796/1‑149 (MQ=255)                                                       ccATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCaa                                                                                             >  3:307099/1‑149 (MQ=255)                                                           gACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCt                                                                                                             <  7:298105/129‑1 (MQ=255)                                                           gACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCt                                                                                                             >  8:298105/1‑129 (MQ=255)                                                             ccccccTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CACCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGcc                                                                                       >  4:118028/1‑149 (MQ=255)                                                              cccccTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGccc                                                                                      <  1:269188/149‑1 (MQ=255)                                                                     aaCCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTACTCCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTg                                                                                <  1:86826/149‑1 (MQ=255)                                                                        cTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCg                                                                            >  1:363745/1‑149 (MQ=255)                                                                         tGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGt                                                                           >  3:141004/1‑149 (MQ=255)                                                                          gTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTg                                                                          >  4:232752/1‑149 (MQ=255)                                                                          gTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAAAGGGCAACTCGCCCCCGGTGCCCGTg                                                                          >  4:237599/1‑149 (MQ=255)                                                                           tATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGt                                                                           <  5:80116/147‑1 (MQ=255)                                                                           tATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑c                                                                        <  5:36321/149‑1 (MQ=255)                                                                             tgtgAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGGGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑cc                                                                       <  5:321175/148‑1 (MQ=255)                                                                               tgAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGa                                                                     >  6:120735/1‑148 (MQ=255)                                                                                gAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAAAGGGGAACTCCGCCCCGGACCCCCTC‑CCGaat                                                                   >  2:147358/1‑148 (MQ=255)                                                                                 aaCAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCg                                                                                   <  3:101118/133‑1 (MQ=255)                                                                                 aaCAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCg                                                                                   >  4:101118/1‑133 (MQ=255)                                                                                   cAGTTGGCTCAATGCCCCGGCCCACACCTCCTGGCGCGTGGCTGAAGCA‑CAGCGCCTGCT‑GCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCAt                                                               <  6:211033/149‑1 (MQ=255)                                                                                     gTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATc                                                              <  4:156465/149‑1 (MQ=255)                                                                                     gTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATc                                                              <  3:195591/149‑1 (MQ=255)                                                                                           tCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGATACCATCAACAcc                                                        >  7:49553/1‑149 (MQ=255)                                                                                           tCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCCTCaccacc                                                        >  7:43890/1‑149 (MQ=255)                                                                                             aaagCCCCGGCCCACACCGCCAGGCGCCTGCCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCCCGAAACTAC‑CCGACGGTTACGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CAGAAACCATCAACAccgc                                                      >  8:193812/4‑149 (MQ=255)                                                                                                 ccccGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACAccgcccgc                                                  <  5:238052/149‑1 (MQ=255)                                                                                                       cccACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTGCCCGAAACCATCAACACCGCCCGCATGAc                                             <  1:22124/149‑1 (MQ=255)                                                                                                           cacCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCAGGTGCCCGTG‑CCGCAACCATCAACACCGCCCGCATGACCCAtt                                        >  1:113328/1‑149 (MQ=255)                                                                                                            acCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCG‑CCCCGGTGCCCGTG‑CCGAAACCATCAACAccgcccgc                                                  <  2:311888/137‑1 (MQ=255)                                                                                                              cgcctggcgcctggCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCt                                     >  3:360430/1‑149 (MQ=255)                                                                                                              cgcctggcgcctggCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCt                                     >  5:15963/1‑149 (MQ=255)                                                                                                                cctggcgcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCttt                                   <  1:147358/149‑1 (MQ=255)                                                                                                                cctggcgcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCttt                                   <  5:82703/149‑1 (MQ=255)                                                                                                                   ggcgcctggcCGAAGCA‑CCGCGCCGGCTGGCCTGTGCCCAGGCCGCGAAACTAC‑CCGACGGGTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGcc                                <  8:43890/149‑1 (MQ=255)                                                                                                                   gccgcgtggcTGAACCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTCC‑CCGACGGCTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGcc                                <  8:49553/147‑1 (MQ=255)                                                                                                                    gcgcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGccc                               >  3:41643/1‑149 (MQ=255)                                                                                                                    gcgcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGc                                 <  4:760/147‑1 (MQ=255)                                                                                                                      gcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTg                                  <  1:366006/144‑1 (MQ=255)                                                                                                                      gcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTg                                  >  2:366006/1‑144 (MQ=255)                                                                                                                      gcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTg                             <  6:15963/149‑1 (MQ=255)                                                                                                                       cctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTgg                            <  3:118028/149‑1 (MQ=255)                                                                                                                             tGAAGCA‑CAGCGCCTGCTGGCCTTTGCCTAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATAGCTTTGCCCTGGCCCAcc                      >  5:356990/1‑149 (MQ=255)                                                                                                                                 gCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTg                    >  3:113002/1‑147 (MQ=255)                                                                                                                                 gCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCa                  <  5:120735/149‑1 (MQ=255)                                                                                                                                  ca‑caGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAg                 <  4:141004/149‑1 (MQ=255)                                                                                                                                   a‑caGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAgg                <  4:360430/149‑1 (MQ=255)                                                                                                                                   a‑caGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAgg                <  3:237599/149‑1 (MQ=255)                                                                                                                                   a‑caGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAgg                <  3:232752/149‑1 (MQ=255)                                                                                                                                         gcCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATc           <  2:363745/149‑1 (MQ=255)                                                                                                                                            tgctgGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAg                 <  6:327560/140‑1 (MQ=255)                                                                                                                                            tgctgGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAg                 >  5:327560/1‑140 (MQ=255)                                                                                                                                                 gCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTa    <  4:113002/148‑1 (MQ=255)                                                                                                                                                 gCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTAg   <  4:41643/149‑1 (MQ=255)                                                                                                                                                  ccTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTAGc  >  5:341980/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CGCCATGACCCCGCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTAGC  >  NC_002947/1156457‑1156749

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