Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,610,696	G→C	5.2%	S20S (TCG→TCC)	emrE →	small multidrug resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,610,696	0	G	C	5.2%	250.8 / 3.5	77	S20S (TCG→TCC)	emrE	small multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base C (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

ATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGC  >  NC_002947/5610549‑5610844                                                                                                                                                     |                                                                                                                                                     aTCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGAt                                                                                                                                                    <  2:20604/150‑1 (MQ=255)          tctcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAaggcggt                                                                                                                                           >  1:807004‑M2/1‑143 (MQ=255)                atatCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAgg                                                                                                                                     >  2:210879/1‑150 (MQ=255)                 tatCGATAGGTCATGACTTGACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTCA‑GAAGGCGCAGGCGCATGAATGCCTATACCTACCTCGCCCTTGACATATGCGTCGACGGCGTCGCCACCGCCTCGATGAAGGCGGTGAAAGGc                                                                                                                                    <  1:875989/150‑1 (MQ=255)                          gTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACg                                                                                                                           >  1:28987/1‑150 (MQ=255)                              tGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGt                                                                                                                       <  1:125716/150‑1 (MQ=255)                               gACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATTCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGaaa                                                                                                                                       >  1:1009075/1‑133 (MQ=255)                               gACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATTCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGaaa                                                                                                                                       <  2:1009075/133‑1 (MQ=255)                                       cccATCCAGGCAACGTCTAGAATGCCCTCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGCTGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActgc                                                                                                              >  1:1118963/1‑150 (MQ=255)                                             cAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCAGCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActgctgcggg                                                                                                        >  1:260579/1‑147 (MQ=255)                                               ggCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtg                                                                                                      <  1:363095/150‑1 (MQ=255)                                                gCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtgg                                                                                                     <  1:1077729/150‑1 (MQ=255)                                                     gTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGtt                                                                                                >  2:242001/1‑150 (MQ=255)                                                      tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTg                                                                                               >  1:492457/1‑150 (MQ=255)                                                       cTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTgg                                                                                              >  1:803298/1‑150 (MQ=255)                                                       cTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTgg                                                                                              <  2:128543/150‑1 (MQ=255)                                                            aTGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGCAAAACCTAAGCACACCGTTGCCACTGCTGCTGATGGTGGTTGGCTATg                                                                                         <  2:1137217/150‑1 (MQ=255)                                                            aTGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATCAAGGCGGTGAAAGGCCt                                                                                                                                  <  2:617627/109‑1 (MQ=255)                                                            aTGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATAAAGGCGGTGAAAGGCCt                                                                                                                                  >  1:617627/1‑109 (MQ=255)                                                                ccGCTCAGTAACGGAGGATCCTGC‑CATCAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAGGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGTCACTGCTGCTGATGGTGGTTGGCTATGGcct                                                                                     >  1:794102/1‑148 (MQ=255)                                                                 cGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAg                                                                                                                                      <  2:885537/100‑1 (MQ=255)                                                                 cGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAg                                                                                                                                      >  1:885537/1‑100 (MQ=255)                                                                        tAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGCCTATGGCATTGCGTTCt                                                                             <  2:1142073/150‑1 (MQ=255)                                                                        tAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCACCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCt                                                                             <  2:1118963/150‑1 (MQ=255)                                                                          aCGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGa                                                                                                                                <  2:605245/97‑1 (MQ=255)                                                                          aCGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGa                                                                                                                                >  1:605245/1‑97 (MQ=255)                                                                          aCGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTg                                                                                                                     <  1:291613/108‑1 (MQ=255)                                                                          aCGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTg                                                                                                                     >  2:291613/1‑108 (MQ=255)                                                                                gaTCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTg                                                                     >  1:1057070/1‑150 (MQ=255)                                                                                    cTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGGGAAAGGCCTGAGCACGCCGTTGCCACTGCAGCTGATGGTGGTTTGGTATGGCTTTGCGTTCTGGATGCTGACGc                                                                 >  1:984451/1‑150 (MQ=255)                                                                                     tGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGCTGCTGACGCt                                                                >  2:456306/1‑150 (MQ=255)                                                                                       c‑cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtgg                                                              <  1:1105883/150‑1 (MQ=255)                                                                                         cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGTCGGtgtt                                                             >  2:991121/1‑148 (MQ=255)                                                                                         cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggt                                                             <  2:1057070/150‑1 (MQ=255)                                                                                         cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGAATCAGCTCCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCttgt                                                             <  1:991121/150‑1 (MQ=255)                                                                                          aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGt                                                                                                                       >  1:159543/1‑91 (MQ=255)                                                                                          aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGt                                                                                                                       <  2:159543/91‑1 (MQ=255)                                                                                          aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                                            <  2:803298/150‑1 (MQ=255)                                                                                           tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACg                                                                  <  2:552849/143‑1 (MQ=255)                                                                                           tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACg                                                                  >  1:552849/1‑143 (MQ=255)                                                                                              aTGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGc                                                                 >  1:972008/1‑141 (MQ=255)                                                                                              aTGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGc                                                                 <  2:972008/141‑1 (MQ=255)                                                                                              aTGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTcc                                                        >  2:235321/1‑150 (MQ=255)                                                                                                 ccTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgca                                                     <  1:235321/150‑1 (MQ=255)                                                                                                 ccTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgca                                                     <  1:242001/150‑1 (MQ=255)                                                                                                  cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGaggcggtggc                                                                                                                                        <  1:443885‑M2/66‑11 (MQ=255)                                                                                                  cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGaggcggtggc                                                                                                                                        >  2:443885‑M2/1‑56 (MQ=255)                                                                                                  cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATgg                                                                                        >  1:516543/1‑114 (MQ=255)                                                                                                  cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATgg                                                                                        <  2:516543/114‑1 (MQ=255)                                                                                                  cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                                                    >  2:707452/1‑150 (MQ=255)                                                                                                      acctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCAtt                                                >  2:20524/1‑150 (MQ=255)                                                                                                        ctaccgCGCCGTTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGGGAAACGCCTGAGGACGCCGTTGCTACTGCGGCGGATGGTGGTTGGCTCTGGCATTGCGTTCTCGATGCTGACGTTGGGGGTCTTGAGCATTcc                                              >  2:504128/1‑150 (MQ=255)                                                                                                           cctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGaggcggtggcga                                                                                                                                      <  1:1037805‑M2/59‑13 (MQ=255)                                                                                                           cctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTTCATCGaggcggtggcga                                                                                                                                      >  2:1037805‑M2/1‑47 (MQ=255)                                                                                                           cctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGt                                           >  1:431573/1‑150 (MQ=255)                                                                                                            ctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTg                                          >  1:108395/1‑150 (MQ=255)                                                                                                            ctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTg                                          >  2:691775/1‑150 (MQ=255)                                                                                                            ctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTg                                          <  1:66000/150‑1 (MQ=255)                                                                                                             tCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGa                                                                    >  1:976954/1‑123 (MQ=255)                                                                                                             tCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGa                                                                    <  2:976954/123‑1 (MQ=255)                                                                                                               gCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgggg                                       >  1:779895/1‑150 (MQ=255)                                                                                                                    tGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGc                                  <  2:492457/150‑1 (MQ=255)                                                                                                                    ggCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGc                                  <  2:151078/149‑1 (MQ=255)                                                                                                                    ggCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGc                                  <  2:141815/149‑1 (MQ=255)                                                                                                                       cATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTa                               <  1:210879/150‑1 (MQ=255)                                                                                                                              gcCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATc                        >  2:1089717/1‑150 (MQ=255)                                                                                                                               ccGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCt                       >  2:900350/1‑150 (MQ=255)                                                                                                                                cGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTg                      >  1:235076/1‑150 (MQ=255)                                                                                                                                        tCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGc              >  2:78910/1‑150 (MQ=255)                                                                                                                                           ccACCGCCTCGATGAaggcggtggcgatgacttcggcgcagatggcaatggcgaggta                                                                                                       >  2:868978‑M2/1‑15 (MQ=255)                                                                                                                                           ccACCGCCTCGATGAaggcggtggcgatgacttcggcgcagatggcaatggcgaggta                                                                                                       <  1:868978‑M2/58‑44 (MQ=255)                                                                                                                                           ccACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTgg           <  2:431573/150‑1 (MQ=255)                                                                                                                                            cACCGCCTCGATGAaggcggtggcgatgacttcggcgcagatggcaatggcgaggtaggtataggcattcatggcaggatcctccgttactgagcgggcattctagacgttgcctggatggggtaaagtcattacctatcgatatgggag          <  1:750219‑M2/150‑137 (MQ=255)                                                                                                                                              ccGCCTCGATGAaggcggtggcgatgacttcggcgcagatggcaatggcgaggtaggtataggcattcatggcag                                                                                   >  1:1185426‑M2/1‑12 (MQ=255)                                                                                                                                              ccGCCTCGATGAaggcggtggcgatgacttcggcgcagatggcaatggcgaggtaggtataggcattcatggcag                                                                                   <  2:1185426‑M2/75‑64 (MQ=255)                                                                                                                                                 ccTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGATGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATtg     <  2:175457/150‑1 (MQ=255)                                                                                                                                                    cGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGc  <  2:361415/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGC‑CATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGC  >  NC_002947/5610549‑5610844

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