Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,701,319	C→T	6.0%	A26T (GCG→ACG)	kdpC ←	potassium‑transporting ATPase subunit C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,701,319	0	C	T	6.0%	197.3 / 2.2	66	A26T (GCG→ACG)	kdpC	potassium‑transporting ATPase subunit C
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base T (2/2);  total (31/35)
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

CTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCA  >  NC_002947/4701175‑4701466                                                                                                                                                  |                                                                                                                                                    cTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTTGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGGAAGCGACTTGGGCAATACCGGTCACCGCCAGc                                                                                                                                                 >  8:36524/1‑149 (MQ=255) cTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGc                                                                                                                                                 >  2:163016/1‑149 (MQ=255)        gggCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACa                                                                                                                                          <  8:312403/149‑1 (MQ=255)          gCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAgg                                                                                                                                        <  3:213512/149‑1 (MQ=255)                     aaCCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCAc                                                                                                                             >  7:75973/1‑149 (MQ=255)                       ccAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCAc                                                                                                                             >  2:84900/1‑147 (MQ=255)                       ccAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCAc                                                                                                                             >  7:220855/1‑147 (MQ=255)                       ccAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCAGCGCGT‑ACCAGGCTCCCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATCCCGGTCACCGCCAGCGTCTACAGTGCGACCGTCAc                                                                                                                             >  8:303917/1‑147 (MQ=255)                         aGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAgg                                                                                                                         <  1:73468/149‑1 (MQ=255)                         aGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAgg                                                                                                                         <  1:163016/149‑1 (MQ=255)                          ttgaTCGCCCAGGACTTCCAGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAggg                                                                                                                        >  3:366772/4‑149 (MQ=255)                          ttgaTCGCCCAGGACCTCCAGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAggg                                                                                                                        >  5:333799/4‑149 (MQ=255)                            cATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtc                                                                                                                      <  7:219660/149‑1 (MQ=255)                              tCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtcat                                                                                                                    <  8:220855/149‑1 (MQ=255)                                gCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtcatc                                                                                                                   >  8:259065/1‑148 (MQ=255)                                gCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGtcatc                                                                                                                   <  6:335823/148‑1 (MQ=255)                                  ccTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGTCAGCGGATACAGGGCGCCCGTCACCAGGGTCATcaa                                                                                                                 >  5:9649/1‑148 (MQ=255)                                        aGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGTGACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGc                                                                                                           >  1:276881/1‑149 (MQ=255)                                        aGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGcc                                                                                                          <  1:84900/149‑1 (MQ=255)                                          tCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAg                                                                                                        <  8:180229/149‑1 (MQ=255)                                                gCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCaa                                                                                                  >  7:317675/1‑149 (MQ=255)                                                gCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCaa                                                                                                  >  6:369723/1‑149 (MQ=255)                                                  gATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATg                                                                                                >  1:314879/1‑149 (MQ=255)                                                  gATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATAACGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATg                                                                                                >  2:332130/1‑149 (MQ=255)                                                      aaTGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCgg                                                                                            <  6:14188/149‑1 (MQ=255)                                                        tGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGc                                                                                          <  4:366772/149‑1 (MQ=255)                                                        tGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGc                                                                                          >  7:319150/1‑149 (MQ=255)                                                                cccGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATaa                                                                                  <  1:204054/149‑1 (MQ=255)                                                                cccGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTCGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATaa                                                                                  <  6:106966/149‑1 (MQ=255)                                                                  cGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGTCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAg                                                                                 >  3:69097/1‑148 (MQ=255)                                                                   gCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAg                                                                                 <  6:260484/147‑1 (MQ=255)                                                                   gCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAg                                                                                 >  5:260484/1‑147 (MQ=255)                                                                     aCCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCa                                                                               <  3:43617/147‑1 (MQ=255)                                                                       cTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGt                                                                             <  6:92611/147‑1 (MQ=255)                                                                       cTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGt                                                                             >  5:92611/1‑147 (MQ=255)                                                                           ccTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCAcc                                                                                                                            <  1:270168/96‑1 (MQ=255)                                                                           ccTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCAcc                                                                                                                            >  2:270168/1‑96 (MQ=255)                                                                               gcgcATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATggcggc                                                                    <  8:319150/148‑1 (MQ=255)                                                                                 gcATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGcctc                                                                 <  8:75973/149‑1 (MQ=255)                                                                                   aTCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGcctcct                                                               <  8:319845/149‑1 (MQ=255)                                                                                   aTCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGcctcct                                                               <  2:276881/149‑1 (MQ=255)                                                                                   aTCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGcctcct                                                               <  7:259065/149‑1 (MQ=255)                                                                                      gcgT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGTCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCa                                                             <  4:69097/148‑1 (MQ=255)                                                                                      cccgaaCCAGGCCAACGGTAGCCTGGTACGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGcctcct                                                               >  3:97366/6‑147 (MQ=255)                                                                                      cccgaaCCAGGCCAACGGTAGCCTGGTACGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGcctcct                                                               >  3:92178/6‑147 (MQ=255)                                                                                             cAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAgg                                                      >  6:52814/1‑149 (MQ=255)                                                                                             cAGGCTAAAGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACATGGCGACCGTCACCAGGGTCCTCAACAAAGCCAGGATCAATGCCGGGCGGAGATAAGCAGTCATGGCGGCCTCCTAAGACAAgg                                                      >  2:133124/1‑149 (MQ=255)                                                                                                    ccGTTGGCCTGGTTCGGGAAAGCGAATTGGGCAATACCGGTCACCGTCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAgg                                                <  6:9649/148‑1 (MQ=255)                                                                                                           ccTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGag                                        >  2:171702/1‑149 (MQ=255)                                                                                                           ccTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGACCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGag                                        <  8:317675/149‑1 (MQ=255)                                                                                                               gTTCGGGAAAGCGACTTTGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCATGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGGTCAGGGCGTTGagcagc                                    <  7:303917/149‑1 (MQ=255)                                                                                                               gTTCGGGAAAGCGACTTGGGCAATACCGGGCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGagcagc                                    <  2:238156/149‑1 (MQ=255)                                                                                                                ttCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGagcag                                     >  4:291526/1‑147 (MQ=255)                                                                                                                   gggAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGTGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCGGGt                                <  2:314879/149‑1 (MQ=255)                                                                                                                         gCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCa                          <  5:52814/149‑1 (MQ=255)                                                                                                                         gCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCa                          <  1:332130/149‑1 (MQ=255)                                                                                                                           gACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGc                        <  6:333799/149‑1 (MQ=255)                                                                                                                                gggCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGAt                   >  7:322592/1‑149 (MQ=255)                                                                                                                                 ggCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGAt                   <  7:36524/148‑1 (MQ=255)                                                                                                                                    aaTACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCg               <  8:322592/149‑1 (MQ=255)                                                                                                                                        ccGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCa            >  6:236213/1‑148 (MQ=255)                                                                                                                                         cGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCaaa          >  7:275285/1‑149 (MQ=255)                                                                                                                                           gTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATcc                 >  3:273792/1‑140 (MQ=255)                                                                                                                                           gTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATcc                 <  4:273792/140‑1 (MQ=255)                                                                                                                                               ccGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGcac    <  5:236213/149‑1 (MQ=255)                                                                                                                                                cGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAG‑ATGATCCCGGCAAACGGcacca  >  5:327483/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CTGCCGAGGGCCGCGAGTGGAACCAGCCATCGCCCTGGAAGTCCTGGGCGATCAATGCCGAACCCCGCACCTGGCCTTGCGCATCGCGT‑ACCAGGCTACCGTTGGCCTGGTTCGGGAAAGCGACTTGGGCAATACCGGTCACCGCCAGCGGATACAGGGCGCCCGTCACCAGGGTCATCAACAAAGCCAGGCTCAATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCA  >  NC_002947/4701175‑4701466

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