Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,105,685	G→T	7.4%	V55V (GTG→GTT)	hisD →	histidinol dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,105,685	0	G	T	7.4%	205.6 / 3.9	65	V55V (GTG→GTT)	hisD	histidinol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base T (3/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTG  >  NC_002947/1105538‑1105823                                                                                                                                                      |                                                                                                                                            aaTTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCtgg‑tg                                                                                                                                             >  8:142676/1‑148 (MQ=255)         gTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCAcc                                                                                                                                    <  4:31664/149‑1 (MQ=255)            tcAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAg                                                                                                                                 >  5:217978/1‑149 (MQ=255)              aaCGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCg                                                                                                                               >  1:166108/1‑149 (MQ=255)              aaCGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑GTGAATTCCACCAGCg                                                                                                                               >  2:115108/1‑149 (MQ=255)                cgccgcTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑GTGAATTCCACCAGCgcc                                                                                                                             <  1:115108/149‑3 (MQ=255)                 gccgcTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGttt                                                                                                                            >  6:109193/1‑149 (MQ=255)                 gccgcTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGttt                                                                                                                            >  3:195428/1‑149 (MQ=255)                     cTGATCCGGAGTTCGCCCGCCCTCTCGCTCCTCTGCTGCGCTGGGAACGTGTGTCCGTTGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGAt                                                                                                                        <  2:166108/149‑1 (MQ=255)                         tCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGgtg                                                                                                                    <  1:172345/149‑1 (MQ=255)                           cGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATg                                                                                                                       >  7:267919/1‑144 (MQ=255)                           cGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATg                                                                                                                       <  8:267919/144‑1 (MQ=255)                           cGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGgtgt                                                                                                                   <  3:78888/148‑1 (MQ=255)                              aTTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGAc                                                                                                               >  7:16818/1‑149 (MQ=255)                                  cGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGcc                                                                                                            <  4:195428/148‑1 (MQ=255)                                     ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGt                                                                                                        <  8:181447/149‑1 (MQ=255)                                     ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGt                                                                                                        >  4:157287/1‑149 (MQ=255)                                     ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGt                                                                                                        >  3:46782/1‑149 (MQ=255)                                     ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑GTGAATTCCACCAGCgccgc                                                                                                                           >  6:156052/1‑126 (MQ=255)                                     ccGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑GTGAATTCCACCAGCgccgc                                                                                                                           <  5:156052/130‑5 (MQ=255)                                      cGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTc                                                                                                       >  7:127919/1‑149 (MQ=255)                                          aTCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGAtt                                                                                                   <  8:127919/149‑1 (MQ=255)                                          aTCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGAtt                                                                                                   <  4:46782/149‑1 (MQ=255)                                          aTCAGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCTGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTTACGCCAAGTCGAtt                                                                                                   <  4:277416/149‑1 (MQ=255)                                               gATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTgac                                                                                                <  7:142676/147‑1 (MQ=255)                                                atcatcTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTgacgac                                                                                             <  4:71304/149‑1 (MQ=255)                                                 tcatcTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTTCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGAcc                                                                                            <  2:130688/149‑1 (MQ=255)                                                  catcTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCt                                                                                           >  2:47094/1‑149 (MQ=255)                                                   atcTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCt                                                                                           <  2:285698/148‑1 (MQ=255)                                                      tgctgAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGCCGCTGG‑TGGAATTCACCCAGCGTTTCGATGgt                                                                                                                     >  7:177804/1‑119 (MQ=255)                                                      tgctgAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGACGCTGG‑TGGAATTCACCCAGCGTTTCGATGgt                                                                                                                     <  8:177804/119‑1 (MQ=255)                                                       gctgAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCCCCCAGCGTTTCGAGGGTGTTGACGCCAAGTCGATTGACGACCTGATcc                                                                                      >  6:95733/1‑149 (MQ=255)                                                        ctgAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCt                                                                                     >  8:145470/1‑149 (MQ=255)                                                                ggAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAACATCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATcgc                                                                              >  1:94652/1‑149 (MQ=255)                                                                        tgtCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCt                                                                     <  8:16818/148‑1 (MQ=255)                                                                          tCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTc                                                                                                       <  7:191351/113‑1 (MQ=255)                                                                          tCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTc                                                                                                       >  8:191351/1‑113 (MQ=255)                                                                          tCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTgg                                                                   >  6:270829/1‑148 (MQ=255)                                                                             gATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAg                                                                                                         <  6:179713/108‑1 (MQ=255)                                                                             gATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAg                                                                                                         >  5:179713/1‑108 (MQ=255)                                                                              aTGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTga                                                                                                 <  8:16192/115‑1 (MQ=255)                                                                              aTGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTga                                                                                                 >  7:16192/1‑115 (MQ=255)                                                                              aTGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTg                                                              <  6:217978/149‑1 (MQ=255)                                                                                  cgcgGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGAGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGATGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGgcgc                                                          >  5:218361/1‑149 (MQ=255)                                                                                      gTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGa                                                       >  4:247184/1‑148 (MQ=255)                                                                                      gTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGa                                                       <  3:247184/148‑1 (MQ=255)                                                                                      gTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGAc                                                      <  7:145470/149‑1 (MQ=255)                                                                                          aTCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGt                                                  <  8:28714/149‑1 (MQ=255)                                                                                                    gCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCgg                                        >  2:204013/1‑149 (MQ=255)                                                                                                     cTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCGACGCGCCTGGAACTGGCGCTGACCCGTATTACCCCg                                         <  1:47094/148‑1 (MQ=255)                                                                                                      tCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTc                                      <  3:157287/149‑1 (MQ=255)                                                                                                        gACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCt                                                                     <  6:178442/116‑1 (MQ=255)                                                                                                        gACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCt                                                                     >  5:178442/1‑116 (MQ=255)                                                                                                               tcaAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGGTTGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAgcgc                                >  1:56290/1‑147 (MQ=255)                                                                                                                caAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGAc                                                                                                               >  8:155985/1‑67 (MQ=255)                                                                                                                caAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGAc                                                                                                               <  7:155985/67‑1 (MQ=255)                                                                                                                caAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAgg                            >  7:248670/1‑149 (MQ=255)                                                                                                                 aaGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGc                           >  7:227329/1‑149 (MQ=255)                                                                                                                     ccGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTg                       <  1:204013/149‑1 (MQ=255)                                                                                                                      cgtgcgtgAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTg                       >  7:120801/1‑148 (MQ=255)                                                                                                                          cgtgAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGaaaa                  <  8:120801/149‑1 (MQ=255)                                                                                                                            tgAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAgg                >  7:204301/1‑149 (MQ=255)                                                                                                                             gAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGgc               >  3:102154/1‑149 (MQ=255)                                                                                                                                    agTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCaa         >  2:49908/2‑149 (MQ=255)                                                                                                                                      tGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAAcc       <  8:227329/149‑1 (MQ=255)                                                                                                                                           aaacgcTGG‑GTGAATTCCACCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGtgtg  >  8:192432/4‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                            AATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCTGAGCTGGGAAAGTGTGTCCGATGACGCGGTCAATCAGCGCGTGCTCGACATCATCAAGGCCGTGCGTGAAC‑GCGGTGATGCGGCGCTGG‑TGGAATTCACCCAGCGTTTCGATGGTGTTGACGCCAAGTCGATTGACGACCTGATCCTCGATCGCG‑CGCGCCTGGAACTGGCGCTGACCCGTATTACCCCGGTCCAGCGCGAGGCCCTGGAAAAGGCCGCCAACCGTGTG  >  NC_002947/1105538‑1105823

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