Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,577,655	T→G	6.9%	C24G (TGC→GGC)	PP_3158 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,577,655	0	T	G	6.9%	180.4 / 5.0	58	C24G (TGC→GGC)	PP_3158	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base G (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACCTGCTGAT  >  NC_002947/3577508‑3577794                                                                                                                                                       |                                                                                                                                            cGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCtgtg                                                                                                                                            >  1:109937/1‑149 (MQ=255)    aCTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTgcgc                                                                                                                                         >  7:179686/1‑149 (MQ=255)    aCTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTgcgc                                                                                                                                         >  6:20488/1‑149 (MQ=255)        tttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTg                                                                                                                                      >  7:105352/1‑148 (MQ=255)         ttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGt                                                                                                                                     <  8:179686/148‑1 (MQ=255)         ttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTa                                                                                                                                    >  6:160556/1‑149 (MQ=255)          tCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTAc                                                                                                                                   >  3:244080/1‑149 (MQ=255)            ccATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCTACCCTGTTCCGATGGTTGCAAGGTACTCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACg                                                                                                                                  <  8:295010/148‑1 (MQ=255)            ccATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACg                                                                                                                                  >  8:270317/1‑148 (MQ=255)             cATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGCGCACAGGCATGCCAGGGCACCTGtcagta                                                                                                                                <  6:219781/149‑6 (MQ=38)              aTTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACgcag                                                                                                                               <  1:266494/149‑1 (MQ=255)                    gggCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGcagcagcag                                                                                                                         <  2:109937/149‑1 (MQ=255)                       ctctATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGc                                                                                                                      >  4:114196/1‑149 (MQ=255)                          tATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGcc                                                                                                                     <  6:256792/147‑1 (MQ=255)                          tATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAAGTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGt                                                                                                                   >  1:87417/1‑149 (MQ=255)                            tGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCg                                                                                                                    <  2:271396/146‑1 (MQ=255)                            tGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGt                                                                                                                   <  7:39925/147‑1 (MQ=255)                             gCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTGCTGCGTACAGGCGCACAGGCATGCCAGGGCACCTGTCAgtagcaggcgg                                                                                                                        >  3:137270/1‑140 (MQ=18)                             gCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTGCTGCGTACAGGCGCACAGGCATGCCAGGGCACCTGTCAgtagcaggcgg                                                                                                                        <  4:137270/144‑5 (MQ=18)                                  gCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATgg                                                                                                            <  4:244080/148‑1 (MQ=255)                                  gCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                           <  6:61465/149‑1 (MQ=255)                                       gCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                           >  3:185055/1‑144 (MQ=255)                                       gCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                           <  4:185055/144‑1 (MQ=255)                                             cTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcg                                                                                                <  7:270317/149‑1 (MQ=255)                                                  tCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCGTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaac                                                                                           >  7:218030/1‑149 (MQ=255)                                                  tCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaac                                                                                           >  4:183299/1‑149 (MQ=255)                                                    ccGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacgg                                                                                         >  6:242441/1‑149 (MQ=255)                                                    ccGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacg                                                                                          <  5:160556/148‑1 (MQ=255)                                                       gcgcAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGTCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacggca                                                                                       >  1:379792/1‑148 (MQ=255)                                                           aGCCGCCCAATATCGACCCGTAGCCGAAGGTAGCAAGCTACGCGAATTTCCTACCCCGCCTGCTACTG‑‑‑ACAGTTGCCCTGGCATGCCTGTGCGCCTTTACGCATCAGCAGGTCCGTGACATGGTCAACCAGCTCGGCAACGTCAAAgcc                                                                                  <  7:166119/149‑1 (MQ=255)                                                                   aaTATCGACCCGGTGCCGATGGTTGCAAGTTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCt                                                                          >  6:215466/1‑149 (MQ=255)                                                                      atCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGtt                                                                       >  6:374055/1‑149 (MQ=255)                                                                         gACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCa                                                                    <  8:105352/149‑1 (MQ=255)                                                                           cccGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGTCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGa                                                                  <  2:379792/149‑1 (MQ=255)                                                                                 gCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCg                                                            <  7:72046/149‑1 (MQ=255)                                                                                          tGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacgg                                                                                         <  2:281141/111‑1 (MQ=255)                                                                                            cAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTg                                                 <  3:114196/149‑1 (MQ=255)                                                                                             aaGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTgg                                                <  5:215466/149‑1 (MQ=255)                                                                                             aaGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTgg                                                >  2:102131/1‑149 (MQ=255)                                                                                                   cgcgAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGc                                          >  1:1461/1‑149 (MQ=255)                                                                                                       aaTTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTc                                      <  7:120340/149‑1 (MQ=255)                                                                                                              tccCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATgcgc                               >  6:15830/1‑149 (MQ=255)                                                                                                                ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCTCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGTCCGCCTGGCCACGCTGTCGATgcgc                               <  8:138593/147‑1 (MQ=255)                                                                                                                ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATgcgc                               <  8:218030/147‑1 (MQ=255)                                                                                                                ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGcgcg                              >  4:365268/1‑148 (MQ=255)                                                                                                                      cTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCt                       >  1:93426/1‑149 (MQ=255)                                                                                                                         ccaCTG‑‑‑ACCCGGGCCCTTGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGGGACATGGGCAACCACCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa                    <  5:144422/147‑1 (MQ=255)                                                                                                                          tACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa                    >  2:128631/1‑148 (MQ=255)                                                                                                                          tACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGACTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACAGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa                    >  4:93652/1‑148 (MQ=255)                                                                                                                               cgtaCAGGCGCACAGGCATGCCAGGGCACCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAg                 >  7:43058/4‑149 (MQ=255)                                                                                                                                  aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGc                                                                               >  3:69522/1‑84 (MQ=255)                                                                                                                                  aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGc                                                                               <  4:69522/84‑1 (MQ=255)                                                                                                                                  aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGc                <  2:1461/147‑1 (MQ=255)                                                                                                                                   cAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGGCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTg             <  8:376579/149‑1 (MQ=255)                                                                                                                                         cccTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTAGctc       >  7:140944/1‑148 (MQ=255)                                                                                                                                         cccTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACctg       <  5:15830/149‑1 (MQ=255)                                                                                                                                              gCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACCTGCtgat  >  8:314229/1‑149 (MQ=255)                                                                                                                                                    cTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGCCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCa                                                               >  6:333367/1‑82 (MQ=255)                                                                                                                                                    cTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGCCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCa                                                               <  5:333367/82‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                            CGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACCTGCTGAT  >  NC_002947/3577508‑3577794

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