Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,477,014	C→A	6.8%	A206S (GCA→TCA)	dsdC ←	dsd operon activator; autorepressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,477,014	0	C	A	6.8%	206.3 / 2.8	77	A206S (GCA→TCA)	dsdC	dsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (13/58);  new base A (4/2);  total (17/60)
Fisher's exact test for biased strand distribution p-value = 1.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAA  >  NC_000913/2476866‑2477163                                                                                                                                                     |                                                                                                                                                       cttcCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGcc                                                                                                                                                       <  2:1620575/150‑1 (MQ=255) cttcCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGcc                                                                                                                                                       <  1:1774735/150‑1 (MQ=255)   tcCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctg                                                                                                                                                     <  1:614799/150‑1 (MQ=255)    cccATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCcggt                                                                                                                                                    >  2:1162270/1‑150 (MQ=255)     ccATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCCTGCctgtg                                                                                                                                                   >  2:93670/1‑149 (MQ=255)       aTCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCACAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtctg                                                                                                                                                 <  2:1744650/150‑1 (MQ=255)       aTCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtctg                                                                                                                                                 <  1:1141607/150‑1 (MQ=255)              ccccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGa                                                                                                                                          <  1:1809051/150‑1 (MQ=255)               cccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAg                                                                                                                                         <  2:426866/150‑1 (MQ=255)                   atatCATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCAAc                                                                                                                                     <  2:197362/150‑1 (MQ=255)                      tGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCTATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACTGTATGTCAAGGAGCaaacaa                                                                                                                                  >  1:72695/1‑146 (MQ=255)                      tGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTa                                                                                                                                  <  2:1095677/150‑1 (MQ=255)                           cATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATTCCTGTCTGTCATGGAGCAACGTACAGTg                                                                                                                             >  2:1082178/1‑150 (MQ=255)                             tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGac                                                                                                                           <  2:1718648/150‑1 (MQ=255)                                 ggCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAg                                                                                                                       <  2:1511090/150‑1 (MQ=255)                                  gCAATAACAGCTAAATCAGAAAGATCAAAGCCAATTCCAGACGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCTTACCGGAGTCGTTGCTCCATGCCGGTCTGTCATGGAGCAACGTACAGTGACTCAgg                                                                                                                      <  1:946401/150‑1 (MQ=255)                                         cAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTa                                                                                                              <  1:557873/150‑1 (MQ=255)                                           gCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTAcc                                                                                                            <  1:369024/150‑1 (MQ=255)                                            cTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCg                                                                                                           <  1:260388/150‑1 (MQ=255)                                            cTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCg                                                                                                           <  2:1309700/150‑1 (MQ=255)                                             tAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCgt                                                                                                          >  2:278887/1‑150 (MQ=255)                                             tAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCAACATTCCAGTATGTGCTGGAGGCACGTAAAATGCCCCCAG‑TTAATTAACgt                                                                                                          >  1:1119806/1‑150 (MQ=255)                                              aaaTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCgtg                                                                                                         <  1:697619/150‑1 (MQ=255)                                                 tCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTg                                                                                                      <  2:116555/150‑1 (MQ=255)                                                  cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTgg                                                                                                     <  1:909133/150‑1 (MQ=255)                                                   aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCAAACTATGCCATTCATCCGTACCGGAGTCGTTTATACATGCCTATAAGTCAAGGTGCAACGTAACCTGACACAGGTTTAATTACC‑TGTTGGt                                                                                                    >  1:1102725/1‑150 (MQ=255)                                                      aCGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTaa                                                                                                 <  2:755995/150‑1 (MQ=255)                                                       cGATCAAAGCCAATTCCAGAAGCTGTCGGCAAATTAACCGCATACTGTTGCGCCCAACTATGCCATTCAACCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTaaa                                                                                                <  2:1955138/150‑1 (MQ=255)                                                        gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAg                                                                                               <  2:491286/150‑1 (MQ=255)                                                        gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAg                                                                                               <  1:607116/150‑1 (MQ=255)                                                        gATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAg                                                                                               <  1:1299552/150‑1 (MQ=255)                                                            aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATg                                                                                           <  2:605878/150‑1 (MQ=255)                                                                    ttCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAg                                                                                   <  1:1467309/150‑1 (MQ=255)                                                                    ttCCAGAAGATGTCGGCAAAGTAACCGCATAATGTTGCGCCCAACTCTGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGGCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAg                                                                                   <  1:1282285/150‑1 (MQ=255)                                                                      ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                 <  1:162552/150‑1 (MQ=255)                                                                              tGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACATACAGTGACACAGG‑GTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCgg                                                                         <  1:1449482/150‑1 (MQ=255)                                                                               gtcgtcAAATTAACCGCATAATGTTGCGCCCAACTATGACATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCggg                                                                        <  2:474408/150‑1 (MQ=255)                                                                               gTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCggg                                                                        <  1:1003466/150‑1 (MQ=255)                                                                                   gCAAATTAACCGCATAATGTTGCGCCCATCTATGCCATTCATCCGTTCCGTAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTAc                                                                    <  1:1684242/150‑1 (MQ=255)                                                                                    cAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACa                                                                   <  1:203708/150‑1 (MQ=255)                                                                                           aCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCCTGCCGGTGGGTCATGGAGCAACGGACAGTGACACAGG‑GTAATTACCCTGTTGGGTAAAGCATGTCTTTGAGCGGATTCCGGGCGACATACTggg                                                            >  2:1164702/1‑149 (MQ=255)                                                                                             cGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAg                                                          <  2:243132/150‑1 (MQ=255)                                                                                                    tGTTGCTCCCAACTATGACATTCATCCGTCCCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttct                                                   <  2:1989981/150‑1 (MQ=255)                                                                                                     tttGCGCCCAACTATGCCTTTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttctt                                                  <  1:720037/149‑1 (MQ=255)                                                                                                        gcgcCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcat                                               <  1:721983/150‑1 (MQ=255)                                                                                                         cgcCCAACTATGCCATTCATCCGTACCGGAGTCGGTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatc                                              <  2:1221377/150‑1 (MQ=255)                                                                                                           cccAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatcca                                            <  1:1962763/150‑1 (MQ=255)                                                                                                             cAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatccatc                                          <  1:1703842/150‑1 (MQ=255)                                                                                                                    gCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGt                                   <  1:390809/150‑1 (MQ=255)                                                                                                                      cATTCATCCGTACCGGAGTCGTTTCTCCATGCCGGTCTGTCATGGGGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATTTCTTTTAGCGTATTCCGGGCTACATACAGGCAGGATTTGTTCATCCAACAGAAAGtga                                 >  2:1027054/1‑150 (MQ=255)                                                                                                                        ttCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                               <  2:1427019/150‑1 (MQ=255)                                                                                                                        ttCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGGCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatg                               <  2:993207/150‑1 (MQ=255)                                                                                                                        ttCATCCGTAACGGAGTCGTTGCTCCATGCCTCTGTGTCATGGAGCAACGGACATTGACACAGG‑ATAATTAGCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAACGtgatg                               <  2:324209/150‑1 (MQ=255)                                                                                                                          cATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgag                             <  2:245172/150‑1 (MQ=255)                                                                                                                            tCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                           <  1:1947153/150‑1 (MQ=255)                                                                                                                            tCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                           <  1:511324/150‑1 (MQ=255)                                                                                                                                     ccggcgctGCTCCATGACTCTCTTCCCTCGAGCAACGTACAGTGACACCTC‑ATAATAACCTTGTTGGTTACAGCATGTCTCTGAGGGTATTCCGGGCTACATACTCGCAGGAATTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCt                  <  2:345582/143‑1 (MQ=255)                                                                                                                                       gTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGa                >  1:45561/1‑150 (MQ=255)                                                                                                                                         cGTTGCTCCGGGCCGGGCGGGCATGGTGCGGCGGAGGGTGACACAGG‑GTAATTACCGTGTTGGTTGAAGCATGTGTTTGAGCGGGTTGCGGGGTACGTACTGGCAGGATTTCTTCATCCATCAGCAAGTGATGTGTGCGTTGCGCGGACg              >  2:2009076/1‑150 (MQ=255)                                                                                                                                         cGTTGCTCCATTCCTGTCTGTCATGGAGCAACGTACAGTGACACTGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCTGGCTACATACTGTCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg              <  2:199642/150‑1 (MQ=255)                                                                                                                                         cGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg              <  2:242661/150‑1 (MQ=255)                                                                                                                                         cGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑GTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg              <  2:1889282/150‑1 (MQ=255)                                                                                                                                          gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg             >  1:915921/1‑150 (MQ=255)                                                                                                                                           ttGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgc            >  2:1374224/1‑150 (MQ=255)                                                                                                                                            tGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg           >  1:1318750/1‑150 (MQ=255)                                                                                                                                             gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGGCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc          <  1:278887/150‑1 (MQ=255)                                                                                                                                             gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGCGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc          >  2:879958/1‑150 (MQ=255)                                                                                                                                             gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc          >  2:1895662/1‑150 (MQ=255)                                                                                                                                             gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc          >  1:369957/1‑150 (MQ=255)                                                                                                                                              cTCCATGCCTGTCTGTCATGGAGCAACGTACAGAGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca         <  2:743433/150‑1 (MQ=255)                                                                                                                                               tCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcat        <  2:1815982/150‑1 (MQ=255)                                                                                                                                               gccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcat        <  1:566571/149‑1 (MQ=255)                                                                                                                                                ccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc       <  1:504362/150‑1 (MQ=255)                                                                                                                                                 cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca      <  2:1858493/150‑1 (MQ=255)                                                                                                                                                 cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca      <  1:1082178/150‑1 (MQ=255)                                                                                                                                                 cATGCCTGTCTGTCATGGAGCAACGTAAAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca      <  2:72695/150‑1 (MQ=255)                                                                                                                                                   tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc    <  1:820097/150‑1 (MQ=255)                                                                                                                                                   tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc    <  1:1944663/150‑1 (MQ=255)                                                                                                                                                   tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactca    <  2:1718035/150‑6 (MQ=255)                                                                                                                                                     ccTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa  <  1:663593/150‑1 (MQ=255)                                                                                                                                                     ccTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa  <  2:1907144/150‑1 (MQ=255)                                                                                                                                                     ccTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa  >  1:1212826/1‑150 (MQ=255)                                                                                                                                                     ccTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGACATTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTACATCCAACACAAAGTGATGAGTGAGTTGCGCTGACGGAGCATAATCaa  >  1:1219157/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                       CTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGG‑TTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAA  >  NC_000913/2476866‑2477163

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