Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	27,434	T→C	6.0%	intergenic (‑281/+165)	PP_0022 ← / ← PP_0023	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	27,434	0	T	C	6.0%	202.3 / 2.3	67	intergenic (‑281/+165)	PP_0022/PP_0023	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base C (2/2);  total (33/34)
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.54e-01

CCCATGTATACCGACGCAGTGGATGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCC  >  NC_002947/27285‑27582                                                                                                                                                      |                                                                                                                                                     cccATGTATACCGACGCAGTGGATGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGatat                                                                                                                                                      <  1:1059602/150‑1 (MQ=255)  ccATGGATACCGACGCAGTGGATGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATc                                                                                                                                                     >  2:1296272/1‑150 (MQ=255)       tataCCGACGCAGTGGATGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATACCAGtct                                                                                                                                                >  1:463159/1‑150 (MQ=255)           ccGACGCAGTGGATGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCtcgtc                                                                                                                                            <  1:1709310/150‑1 (MQ=255)             gACGCAGTGGATGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATACCAGTCTcgtcgt                                                                                                                                          >  1:90239/1‑150 (MQ=255)              aCGCAGTGGATGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCgtcgtc                                                                                                                                         >  2:243549/1‑150 (MQ=255)                       ctGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGaa                                                                                                                                >  2:1448979/2‑150 (MQ=255)                        tGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAg                                                                                                                               <  1:243549/150‑1 (MQ=255)                           gATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTcgtcgt                                                                                                                                          >  2:1281125/1‑136 (MQ=255)                           gATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTcgtcgt                                                                                                                                          <  1:1281125/136‑1 (MQ=255)                           gATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGa                                                                                                                            >  2:1453056/1‑150 (MQ=255)                               aaGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAAc                                                                                                                        <  1:779007/150‑1 (MQ=255)                                 gCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCg                                                                                                                      >  1:1785439/1‑150 (MQ=255)                                 gCAATAACACCAGCCGGAAGGGTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCGGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCg                                                                                                                      >  2:478237/1‑150 (MQ=255)                                    aTACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATACCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATc                                                                                                                   <  2:90239/150‑1 (MQ=255)                                         aCGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAg                                                                                                              >  1:260000/1‑150 (MQ=255)                                         aCGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAg                                                                                                              <  1:1032602/150‑1 (MQ=255)                                         aCGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAg                                                                                                              <  2:1084792/150‑1 (MQ=255)                                         aCGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAg                                                                                                              <  1:1453056/150‑1 (MQ=255)                                              gcgcAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGc                                                                                                         >  2:1283830/1‑150 (MQ=255)                                                 cAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCtt                                                                                                       <  2:1682438/149‑1 (MQ=255)                                                 cAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCtt                                                                                                       >  1:1682438/1‑149 (MQ=255)                                                 cAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTc                                                                                                      <  2:260000/150‑1 (MQ=255)                                                    gCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGa                                                                                                   <  1:607618/150‑1 (MQ=255)                                                     cTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGaa                                                                                                  <  2:144861/150‑1 (MQ=255)                                                         tGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATACCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCa                                                                                              <  2:463159/150‑1 (MQ=255)                                                          ggCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCAt                                                                                             <  2:1743085/150‑1 (MQ=255)                                                               aTCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATcaagc                                                                                        >  1:241709/1‑150 (MQ=255)                                                                  tGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAg                                                                                                                               <  2:1483852/108‑1 (MQ=255)                                                                  tGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAg                                                                                                                               >  1:1483852/1‑108 (MQ=255)                                                                    ccGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTcc                                                                                                                                        >  1:422707/1‑97 (MQ=255)                                                                    ccGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTcc                                                                                                                                        <  2:422707/97‑1 (MQ=255)                                                                      gAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATcaagccaagc                                                                                   >  1:772677/1‑148 (MQ=255)                                                                      gAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATcaagccaagc                                                                                   <  2:772677/148‑1 (MQ=255)                                                                            gTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATcc                                                                                                                  >  1:598882/1‑111 (MQ=255)                                                                            gTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATcc                                                                                                                  <  2:598882/111‑1 (MQ=255)                                                                               gAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAgg                                                                                                             >  2:1208357/1‑113 (MQ=255)                                                                               gAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAgg                                                                                                             <  1:1208357/113‑1 (MQ=255)                                                                                aGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAAcc                                                                                                                       <  1:226508/102‑1 (MQ=255)                                                                                aGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAAcc                                                                                                                       >  2:226508/1‑102 (MQ=255)                                                                                     cTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTcc                                                                                                                                        >  2:1619562/1‑80 (MQ=255)                                                                                     cTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTcc                                                                                                                                        <  1:1619562/80‑1 (MQ=255)                                                                                          tCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCg                                                                                                                 >  1:1598108/1‑98 (MQ=255)                                                                                          tCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCg                                                                                                                 <  2:1598108/98‑1 (MQ=255)                                                                                             aaTCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGAtt                                                          >  1:1041303/1‑150 (MQ=255)                                                                                              aTCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTg                                                         <  2:241709/150‑1 (MQ=255)                                                                                                  cAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTGTCGTCGTCCGCACt                                                                                                                                   >  1:726977/1‑72 (MQ=255)                                                                                                  cAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACt                                                                                                                                   <  2:726977/72‑1 (MQ=255)                                                                                                    gAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGc                                                   <  2:1297191/150‑1 (MQ=255)                                                                                                         acaATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCtt                                                                                                       <  2:1226057/93‑1 (MQ=255)                                                                                                         acaATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCtt                                                                                                       >  1:1226057/1‑93 (MQ=255)                                                                                                                   gCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAgg                                    >  2:1374763/1‑150 (MQ=255)                                                                                                                    cTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATcaagcc                                                                                       <  1:1675086/98‑1 (MQ=255)                                                                                                                    cTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATcaagcc                                                                                       >  2:1675086/1‑98 (MQ=255)                                                                                                                       tCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATcaa                                                                                          <  2:1238239/92‑1 (MQ=255)                                                                                                                       tCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATcaa                                                                                          >  1:1238239/1‑92 (MQ=255)                                                                                                                       tCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGttt                                <  1:478237/150‑1 (MQ=255)                                                                                                                        cAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTg                               <  2:1785439/150‑1 (MQ=255)                                                                                                                          gCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATcaagc                                                                                        <  2:1440217/91‑1 (MQ=255)                                                                                                                          gCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATcaagc                                                                                        >  1:1440217/1‑91 (MQ=255)                                                                                                                                 tCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGcc                      >  2:1514239/1‑150 (MQ=255)                                                                                                                                      tCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATaa                                                                                                                         >  2:1714700/1‑46 (MQ=255)                                                                                                                                      tCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATaa                                                                                                                         <  1:1714700/46‑1 (MQ=255)                                                                                                                                      tCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCaaa                 >  2:492957/1‑150 (MQ=255)                                                                                                                                       cAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCaaaa                >  2:1930469/1‑150 (MQ=255)                                                                                                                                              ctGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAagc         <  1:1283830/150‑1 (MQ=255)                                                                                                                                               tGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGAc                                              >  1:1243588/1‑112 (MQ=255)                                                                                                                                               tGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGAc                                              <  2:1243588/112‑1 (MQ=255)                                                                                                                                                    tatCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTcc   >  1:1767567/1‑150 (MQ=255)                                                                                                                                                     atCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTccc  >  1:1907987/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCCATGTATACCGACGCAGTGGATGCGATGAAGCAATACCACGAGGCGCAGGCTTCTGGCTCATCTGCCGAAGAGGTTGAGCGCCTGCGTCAAATCGCAGAATCACAATTTCGGGCTGTCAGCGAGTATCAGCTCAACGCTCTGGGATATCAGTCTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCC  >  NC_002947/27285‑27582

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