Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,156,364	G→A	6.0%	D157D (GAC→GAT)	PP_3659 ←	transcriptional regulator, AraC family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,156,364	0	G	A	6.0%	191.4 / 2.6	66	D157D (GAC→GAT)	PP_3659	transcriptional regulator, AraC family
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base A (2/2);  total (33/33)
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

CAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGTCATCAAGGGCAGTCCT  >  NC_002947/4156215‑4156513                                                                                                                                                      |                                                                                                                                                       cAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTg                                                                                                                                                        >  1:1640908/1‑150 (MQ=255)  aaaCCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGCATGCGCCTGt                                                                                                                                                       <  2:968272/150‑1 (MQ=255)  aaaCCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGt                                                                                                                                                       >  1:1310160/1‑150 (MQ=255)     ccTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGTGCCTGTCGc                                                                                                                                                    >  1:583526/1‑150 (MQ=255)              ccAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCAcagccagc                                                                                                                                           <  1:1080422/150‑1 (MQ=255)              ccAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCAcagccagc                                                                                                                                           <  2:806635/150‑1 (MQ=255)                    aGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCagcagagca                                                                                                                                     >  1:812529/1‑150 (MQ=255)                       tcttctCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCCCAGCCAGCAGAGCAGGc                                                                                                                                  >  2:76273/1‑150 (MQ=255)                       tcttctCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGc                                                                                                                                  <  2:477695/150‑1 (MQ=255)                        cttctCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGc‑c                                                                                                                                <  2:1864478/150‑1 (MQ=255)                              gAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCAc                                                                                                                          >  1:1358181/1‑150 (MQ=255)                                acacAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGa                                                                                                                        <  1:1873706/150‑1 (MQ=255)                                  acaGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTGACGATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGc                                                                                                                      >  1:417024/1‑150 (MQ=255)                                    aGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCAc                                                                                                                    <  2:1640908/150‑1 (MQ=255)                                       ggAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCAccag                                                                                                                 >  1:255936/1‑150 (MQ=255)                                           gTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCAccagccag                                                                                                             >  1:1250731/1‑150 (MQ=255)                                                       cTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGTGCCTGTCGCACAGCCAGCAGAGCGGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACaa                                                                                                 <  2:583526/150‑1 (MQ=255)                                                          cACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAac                                                                                              >  2:1601266/1‑150 (MQ=255)                                                                 tttCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTc                                                                                       <  1:505433/150‑1 (MQ=255)                                                                  ttCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCg                                                                                      <  2:529969/150‑1 (MQ=255)                                                                     cTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCaa                                                                                   >  1:812490/1‑150 (MQ=255)                                                                      tCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCcagca                                                                                                                                          <  2:596233/95‑1 (MQ=255)                                                                      tCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCcagca                                                                                                                                          >  1:596233/1‑95 (MQ=255)                                                                           aaaGCGCATTGGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGCGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCAcc                                                                             >  2:788234/1‑150 (MQ=255)                                                                            aaGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACaaaa                                                                                               >  2:442249/1‑131 (MQ=255)                                                                            aaGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACaaaa                                                                                               <  1:442249/131‑1 (MQ=255)                                                                              gcgcATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATg                                                                          >  2:1339955/1‑150 (MQ=255)                                                                                       aTAGGCGATGTCGTCAACCGGCTGAGGAGGACGCAGCGATAACGTCACAATAGGAATGCGCCTGTCGCACAGTCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCATGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCattaat                                                                 <  2:191176/150‑1 (MQ=255)                                                                                         aGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCatgaatga                                                               <  2:812529/150‑1 (MQ=255)                                                                                           tCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATa                                                             >  1:826543/1‑150 (MQ=255)                                                                                              cTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAAtt                                                          <  1:1540892/150‑1 (MQ=255)                                                                                                   gtcCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATCAATGATAATTgcagc                                                     <  2:255936/150‑1 (MQ=255)                                                                                                    tcCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTgcagca                                                    >  1:1345447/1‑150 (MQ=255)                                                                                                      cACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATg                                                  >  2:1757388/1‑150 (MQ=255)                                                                                                              gaggagGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCt                                          <  2:1250731/150‑1 (MQ=255)                                                                                                               aggaggACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTc                                         <  1:629834/150‑1 (MQ=255)                                                                                                                 gaggaCGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAAc                                                                                  <  1:444966/107‑1 (MQ=255)                                                                                                                 gaggaCGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAAc                                                                                  >  2:444966/1‑107 (MQ=255)                                                                                                                 gaggaCGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGctc                                               >  2:321099/1‑142 (MQ=255)                                                                                                                 gaggaCGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGctc                                               <  1:321099/142‑1 (MQ=255)                                                                                                                       gCATCGATAACTTCACAATAGGAATGCGCCTGTCGCCCCGCC‑GTTCCTCGTGCGCGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCa                                 <  1:1868470/150‑1 (MQ=255)                                                                                                                       gCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGGCTTATCTCCCTCggctgcga                                 >  1:574771/1‑143 (MQ=255)                                                                                                                        cATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCAc                                <  1:452923/150‑1 (MQ=255)                                                                                                                         aTCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCAcca                                                                                                                  <  2:1306869/67‑1 (MQ=255)                                                                                                                         aTCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCAcca                                                                                                                  >  1:1306869/1‑67 (MQ=255)                                                                                                                         aTCGATAACTTCACAATAGGAATGCGCCTATCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCAcca                                                                                                                  >  2:517824/1‑67 (MQ=255)                                                                                                                         aTCGATAACTTCACAATAGGAATGCGCCTATCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCAcca                                                                                                                  >  2:512273/1‑67 (MQ=255)                                                                                                                         aTCGATAACTTCACAATAGGAATGCGCCTATCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCAcca                                                                                                                  <  1:517824/67‑1 (MQ=255)                                                                                                                         aTCGATAACTTCACAATAGGAATGCGCCTATCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCAcca                                                                                                                  <  1:512273/67‑1 (MQ=255)                                                                                                                            gATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGAc                            >  2:232373/1‑150 (MQ=255)                                                                                                                                aCTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGt                        <  1:1601266/150‑1 (MQ=255)                                                                                                                                 cTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTa                       >  2:190115/1‑150 (MQ=255)                                                                                                                                 cTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTa                       >  1:714923/1‑150 (MQ=255)                                                                                                                                 cTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTa                       <  1:1339955/150‑1 (MQ=255)                                                                                                                                  ttCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTTCACCGACACGTAg                      >  1:853752/1‑150 (MQ=255)                                                                                                                                       aaTAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGtc                 >  2:1614068/1‑150 (MQ=255)                                                                                                                                        aTAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGtca                <  2:714923/150‑1 (MQ=255)                                                                                                                                        aTAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGtca                <  2:1345447/150‑1 (MQ=255)                                                                                                                                         tAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGtcat               <  1:1115879/150‑1 (MQ=255)                                                                                                                                            gAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGTCATCaa            >  2:1030657/1‑150 (MQ=255)                                                                                                                                                  gcCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCatgaatga                                                               <  2:745755/93‑1 (MQ=255)                                                                                                                                                  gcCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCatgaatga                                                               >  1:745755/1‑93 (MQ=255)                                                                                                                                                  gcCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTa                       <  1:434774/133‑1 (MQ=255)                                                                                                                                                  gcCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTa                       >  2:434774/1‑133 (MQ=255)                                                                                                                                                  gcCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGTCATCAAGGGCAg      <  2:1310160/150‑1 (MQ=255)                                                                                                                                                  gcCTGTCGCACAGCCAGCAGACCCGAC‑GGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGTCATCAAGGGCAg      >  1:1698559/1‑150 (MQ=255)                                                                                                                                                    cTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGTCATCAAGGGCAGTc    >  1:371170/1‑150 (MQ=255)                                                                                                                                                      gTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGTCATCAAGGGCAGTCCt  <  1:1025615/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGCATCGATAACTTCACAATAGGAATGCGCCTGTCGCACAGCCAGCAGAGCAGGC‑CGTGCACGAGCACCAGCCAGGCACCGTAACAAAACACCCTTCGCAACTCACCATGATCATGAATGATAATTGCAGCATGCTCTCCCTCGCATTCCACCGACCCGTAGATGTCATCAAGGGCAGTCCT  >  NC_002947/4156215‑4156513

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