Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,594,103	T→C	7.9%	intergenic (‑259/+42)	PP_1399 ← / ← kgtP	hypothetical protein/alpha‑ketoglutarate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,594,103	0	T	C	7.9%	192.3 / 3.8	72	intergenic (‑259/+42)	PP_1399/kgtP	hypothetical protein/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base T (25/39);  new base C (3/3);  total (28/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTGGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCAT  >  NC_002947/1593958‑1594252                                                                                                                                                  |                                                                                                                                                       aTTTGGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCAAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑t                                                                                                                                                   >  1:617387/1‑150 (MQ=255)  tttGGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑Tg                                                                                                                                                  >  1:1788835/1‑150 (MQ=255)       gcgcGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATg                                                                                                                                             <  1:1529308/150‑1 (MQ=255)         gcgTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGTGCGCTCCA‑TGCAATggg                                                                                                                                           <  2:342211/150‑1 (MQ=255)            tATGTCGTCCTTCAGGCCGCAAGGCTGGCTTGGGAGCGGTCGCGCGGGCGCCAAAGTGCCGCCACCAATCTCATCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGgcg                                                                                                                                        <  1:430968/150‑1 (MQ=255)             aTGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGgcgc                                                                                                                                       >  2:176568/1‑150 (MQ=255)              tGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGcc                                                                                                                                      <  1:794536/150‑1 (MQ=255)                tCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGcccc                                                                                                                                    <  1:972900/150‑1 (MQ=255)                   tGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCttt                                                                                                                                 <  2:1788835/150‑1 (MQ=255)                   tGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCttt                                                                                                                                 <  1:1667225/150‑1 (MQ=255)                     cTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTa                                                                                                                               <  2:1065747/150‑1 (MQ=255)                      ttCAGGCCGCAAGGCTGGCTTGGGCATAATTACGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCACTTTTAg                                                                                                                              >  1:988758/1‑150 (MQ=255)                      ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAg                                                                                                                              <  1:811203/150‑1 (MQ=255)                      ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAg                                                                                                                              >  1:346078/1‑150 (MQ=255)                        cAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTa                                                                                                                            <  1:1014585/150‑1 (MQ=255)                         aGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTa                                                                                                                            >  2:197804/1‑149 (MQ=255)                         aGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTa                                                                                                                            <  1:197804/149‑1 (MQ=255)                          ggCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACt                                                                                                                          >  2:524837/1‑150 (MQ=255)                           gCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTg                                                                                                                         <  2:1670734/150‑1 (MQ=255)                                     tggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGc                                                                                                               <  1:1478226/150‑1 (MQ=255)                                     tggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGc                                                                                                               <  1:21836/150‑1 (MQ=255)                                       gctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCtt                                                                                                                                  >  2:1713158/1‑129 (MQ=255)                                       gctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCtt                                                                                                                                  <  1:1713158/129‑1 (MQ=255)                                         ttGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGC‑ATGGAGCGCCCCttttctattc                                                                                                                          <  2:1656557/135‑7 (MQ=37)                                         ttGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCA‑TTGCATGGAGCGCCCCttttctattc                                                                                                                          >  1:1656557/1‑129 (MQ=38)                                          tGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTAc                                                                                                          <  1:806753/150‑1 (MQ=255)                                           gggAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACt                                                                                                         <  1:524837/150‑1 (MQ=255)                                                cGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                    <  1:176568/150‑1 (MQ=255)                                                  ttCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGt                                                                                                  >  1:213306/1‑150 (MQ=255)                                                          ggCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCcgcg                                                                                          >  2:1076371/1‑150 (MQ=255)                                                             gcCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGt                                                                                       >  1:1761557/1‑150 (MQ=255)                                                               cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgc                                                                                     >  1:1306291/1‑150 (MQ=255)                                                               cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgc                                                                                     <  1:749755/150‑1 (MQ=255)                                                                aaCGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgct                                                                                    >  1:10451/1‑150 (MQ=255)                                                                     tCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTg                                                                               <  2:1476791/150‑1 (MQ=255)                                                                         ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgc                                                                                     <  2:1280024/140‑1 (MQ=255)                                                                         ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgc                                                                                     >  1:1280024/1‑140 (MQ=255)                                                                         ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgtgt                                                                           >  1:486506/1‑150 (MQ=255)                                                                            gccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGa                                                                                         <  1:547285/133‑1 (MQ=255)                                                                            gccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGa                                                                                         >  2:547285/1‑133 (MQ=255)                                                                              caTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCANGCCAGCNTTACTCAGTGGTNATCCGCGAGTGCTGCTTGGTGTCTNTc                                                                      <  1:1863933/150‑1 (MQ=255)                                                                              caTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTc                                                                      <  2:297946/150‑1 (MQ=255)                                                                                ttGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCAt                                                                    <  2:486506/150‑1 (MQ=255)                                                                                  tcAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATgg                                                                  <  1:789610/149‑1 (MQ=255)                                                                                    aGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTc                                                                >  1:1414221/1‑150 (MQ=255)                                                                                               gTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcag                                                     >  2:1237940/1‑150 (MQ=255)                                                                                                  gAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagcag                                                  >  1:1882197/1‑150 (MQ=255)                                                                                                     ttAATAGGGCGTTAAGAAGGAAAATTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGTTAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGa                                               >  1:595649/1‑150 (MQ=255)                                                                                                        aTAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCa                                            <  2:213306/150‑1 (MQ=255)                                                                                                          aGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAgg                                          >  1:632194/1‑150 (MQ=255)                                                                                                          aGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAgg                                          <  2:1306291/150‑1 (MQ=255)                                                                                                               gTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCttt                                                                       >  2:805340/1‑116 (MQ=255)                                                                                                               gTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCttt                                                                       <  1:805340/116‑1 (MQ=255)                                                                                                                   agaagGACAAGTCTGACTAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAgg                                 <  2:974305/150‑1 (MQ=255)                                                                                                                    gaagGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGc                                >  2:1515628/1‑150 (MQ=255)                                                                                                                    gaagGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGc                                <  1:1431817/150‑1 (MQ=255)                                                                                                                      agGAAAAGTCTGAATAGAAAAGGGGCGCCCCA‑TGCAATGGGGCGCCCCTTATAGTACTGCCAGGCCAGCGTTTCCCCGTCCACCAAAGCCCATGTAACGTTTTGGCGGTCATCGGCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCgg                              <  1:35993/150‑1 (MQ=255)                                                                                                                       ggAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                             <  2:346078/150‑1 (MQ=255)                                                                                                                       ggAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                             <  1:1515628/150‑1 (MQ=255)                                                                                                                         aaaaGTCTGAATAGAAAAGGGGCGCGCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGa                           <  2:329873/150‑1 (MQ=255)                                                                                                                          aaaGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCGGCGAGTGGTACTTGGTGTCTTTAGTTGTGGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGAc                          <  2:66327/150‑1 (MQ=255)                                                                                                                            aGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGt                        <  2:1761557/150‑1 (MQ=255)                                                                                                                              tCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGAGGCGCCCCTTTCAGATCTCAGTTTTTTGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTCGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                      <  2:988758/150‑1 (MQ=255)                                                                                                                              tCTGAATAGAAAAGGGGCGCGCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                      <  2:768775/150‑1 (MQ=255)                                                                                                                                  aaTAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTaaa                                                         >  1:1123678/1‑111 (MQ=255)                                                                                                                                  aaTAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTaaa                                                         <  2:1123678/111‑1 (MQ=255)                                                                                                                                   aTAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAg                   <  1:1834439/148‑1 (MQ=255)                                                                                                                                   aTAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAg                   >  2:1834439/1‑148 (MQ=255)                                                                                                                                      taaaaGGGGCGCCCCATTGC‑ATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                       <  1:1089566/124‑1 (MQ=38)                                                                                                                                      taaaaGGGGCGCCCCA‑TTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                       >  2:1089566/2‑125 (MQ=38)                                                                                                                                            gggCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGct        <  1:524315/150‑1 (MQ=255)                                                                                                                                            gggCGCCCCA‑TTGTATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTCACGTACCAGTAGAAGCCGct        >  2:1830790/1‑150 (MQ=255)                                                                                                                                                  tCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCAt  >  2:371590/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                       ATTTGGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCA‑TGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCAT  >  NC_002947/1593958‑1594252

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