Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	281,599	C→T	100%	G129S (GGC→AGC)	yagA ←	CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	281,599	0	C	T	97.7%	271.5 / ‑2.2	87	G129S (GGC→AGC)	yagA	CP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (38/46);  total (40/47)
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.88e-01

GGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCT  >  NC_000913/281503‑281694                                                                                                 |                                                                                                ggCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCgg                                                                                              >  2:286/1‑100 (MQ=255)  gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCTAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCgg                                                                                              >  2:1324467/1‑99 (MQ=255)  gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGt                                                                                             <  1:118814/100‑1 (MQ=255)  gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGt                                                                                             <  1:631218/100‑1 (MQ=255)      cctCCACCGCCAAAAGGACAGTGGCCCTTAAGATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGCGGCGGGAATGCTCGGTg                                                                                            <  2:1420842/95‑1 (MQ=255)      cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAg                                                                                         >  2:1220518/1‑100 (MQ=255)       gTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgc                                                                                        <  1:705164/100‑1 (MQ=255)       gTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgc                                                                                        >  1:757013/1‑100 (MQ=255)         ccACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgcgc                                                                                      >  1:1927212/1‑100 (MQ=255)         ccACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgcgc                                                                                      >  2:1313531/1‑100 (MQ=255)          cACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGcc                                                                                     <  1:68233/100‑1 (MQ=255)           aCCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGccc                                                                                    <  2:1235096/100‑1 (MQ=255)             cGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCgg                                                                                  >  2:1001406/1‑100 (MQ=255)               ccAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgca                                                                                <  1:270395/100‑1 (MQ=255)                cAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgcag                                                                               >  2:1216383/1‑100 (MQ=255)                 aaaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgcagc                                                                              <  1:130305/100‑1 (MQ=255)                   aaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGcagcag                                                                            <  2:234101/100‑1 (MQ=255)                    nggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAgg                                                                           >  1:1408839/2‑100 (MQ=255)                    aGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAgg                                                                           >  2:1853182/1‑100 (MQ=255)                     ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCGCGGTGAAGCGCCCGGCAGCAGGc                                                                          >  2:32774/1‑100 (MQ=255)                     ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCATGGCGGGAATGCTCTGTGAAGCGCCCGGCAGCAgg                                                                           <  2:1458911/99‑1 (MQ=255)                      gAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGcc                                                                         >  1:2252137/1‑100 (MQ=255)                        aaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCAt                                                                       <  1:1289549/100‑1 (MQ=255)                           tGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGc                                                                    >  2:1360213/1‑100 (MQ=255)                            ggCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCg                                                                   <  2:1938883/100‑1 (MQ=255)                              cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                                                 <  2:1751055/100‑1 (MQ=255)                               ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                                <  1:906638/100‑1 (MQ=255)                               ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                                >  1:1754570/1‑100 (MQ=255)                               ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                                >  1:1754198/1‑100 (MQ=255)                               ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                                <  1:1676374/100‑1 (MQ=255)                               ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                                <  2:15410/100‑1 (MQ=255)                                  tAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCAt                                                             >  1:1406990/1‑100 (MQ=255)                                    aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCa                                                           >  1:1129581/1‑100 (MQ=255)                                    aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCa                                                           <  1:1918635/100‑1 (MQ=255)                                      aTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAgg                                                         <  1:2109128/100‑1 (MQ=255)                                       tCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAgg                                                         >  1:621124/1‑99 (MQ=255)                                       tCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAgg                                                         >  2:1629423/1‑99 (MQ=255)                                          aTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTAt                                                     >  2:212010/1‑100 (MQ=255)                                            cTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTCATgg                                                   >  2:256893/1‑100 (MQ=255)                                             tGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGa                                                  >  1:2128315/1‑100 (MQ=255)                                              gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGan                                                 <  2:1408831/100‑2 (MQ=255)                                              gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                                                 >  1:1208314/1‑100 (MQ=255)                                              gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                                                 >  2:22153/1‑100 (MQ=255)                                               ccAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACg                                                <  1:577583/100‑1 (MQ=255)                                                  gagGCGGTTCGGGGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGt                                              <  2:251375/99‑1 (MQ=255)                                                    ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGc                                            <  2:2043151/99‑1 (MQ=255)                                                        gTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAg                                       <  1:286/100‑1 (MQ=255)                                                           cGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCg                                    <  1:1595020/100‑1 (MQ=255)                                                            ggCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                                   >  2:1714062/1‑100 (MQ=255)                                                            ggCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                                   >  2:2136570/1‑100 (MQ=255)                                                            ggCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                                   >  1:2075812/1‑100 (MQ=255)                                                            ggCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                                   >  1:1514249/1‑100 (MQ=255)                                                                cgtcgtGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCAt                               >  1:1189648/1‑100 (MQ=255)                                                                 gtcgtGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATg                              <  2:1754184/100‑1 (MQ=255)                                                                 gtcgtGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATg                              <  2:1754556/100‑1 (MQ=255)                                                                   cgtGTTCGAACCGGCCCGTGGCGGGACGGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATgg                             <  2:2055826/99‑1 (MQ=255)                                                                   cgtGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgt                            <  1:1574720/100‑1 (MQ=255)                                                                   cgtGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGCTATGGACGGTGCTGAAGGCGGGCATGgt                            >  1:118140/1‑100 (MQ=255)                                                                    gtgtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg                           <  2:1129576/100‑1 (MQ=255)                                                                    gtgtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg                           >  1:2270261/1‑100 (MQ=255)                                                                    gtgtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgt                            <  2:521864/99‑1 (MQ=255)                                                                      gtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg                         >  2:115647/1‑100 (MQ=255)                                                                      gtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg                         >  2:85156/1‑100 (MQ=255)                                                                        tCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcc                       <  2:205120/100‑1 (MQ=255)                                                                         cGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGccc                      >  2:385981/1‑100 (MQ=255)                                                                          gAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcccc                     >  2:3823/1‑100 (MQ=255)                                                                          gAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcccc                     <  2:1874069/100‑1 (MQ=255)                                                                            aCCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTg                   >  1:1446345/1‑100 (MQ=255)                                                                             ccGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTgg                  >  1:1705574/1‑100 (MQ=255)                                                                              cGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGt                 <  1:61132/100‑1 (MQ=255)                                                                                gCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTc                <  2:469301/99‑1 (MQ=255)                                                                                 cccGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCGGGTCCt              >  2:622252/1‑100 (MQ=255)                                                                                   cGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGGAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCg            >  2:1553404/1‑100 (MQ=255)                                                                                    gTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCg            <  2:296737/99‑1 (MQ=255)                                                                                    gTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGa           <  1:2073350/100‑1 (MQ=255)                                                                                     tGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAg          <  1:2289496/100‑1 (MQ=255)                                                                                     tGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAg          <  1:254134/100‑1 (MQ=255)                                                                                     tGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAg          <  1:1248023/100‑1 (MQ=255)                                                                                      ggCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGc         <  2:976359/100‑1 (MQ=255)                                                                                       gCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGcc        <  1:1764503/100‑1 (MQ=255)                                                                                         gggAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAg      <  1:1806516/100‑1 (MQ=255)                                                                                         gggAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAg      <  2:1130970/100‑1 (MQ=255)                                                                                          ggAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAgc     <  2:1891010/100‑1 (MQ=255)                                                                                           gAATGCTCGGTGAAGCGCCCGGCAGCAGGCCGTGGCGGGCCAACAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAgcg    >  1:1609161/1‑100 (MQ=255)                                                                                           gAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAgcg    <  1:1761470/100‑1 (MQ=255)                                                                                            aaTGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAgcgc   <  2:1175658/100‑1 (MQ=255)                                                                                             aTGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCt  <  1:218799/100‑1 (MQ=255)                                                                                                 |                                                                                                GGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCT  >  NC_000913/281503‑281694

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