Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,477,016	T→A	12.0%	Q205L (CAG→CTG)	dsdC ←	DNA‑binding transcriptional dual regulator DsdC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,477,016	0	T	A	12.0%	152.1 / 9.3	63	Q205L (CAG→CTG)	dsdC	DNA‑binding transcriptional dual regulator DsdC
Reads supporting (aligned to +/- strand):  ref base T (13/38);  new base A (6/2);  total (23/40)
Fisher's exact test for biased strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCA  >  NC_000913/2476868‑2477162                                                                                                                                                     |                                                                                                                                                     tcCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctg                                                                                                                                                     <  2:599843/150‑1 (MQ=255) tcCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctg                                                                                                                                                     <  2:1022339/150‑1 (MQ=255) tcCCATCGCCACCCAAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctg                                                                                                                                                     <  1:1258979/150‑1 (MQ=255)  cccATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgt                                                                                                                                                    <  1:1424032/150‑1 (MQ=255)   ccATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtc                                                                                                                                                   <  2:555003/150‑1 (MQ=255)   ccATCGCCACCCAAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGCAGTCGTTGAACCATGAcagtc                                                                                                                                                   >  2:995573/1‑150 (MQ=255)    cATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGAATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtct                                                                                                                                                  >  2:1381760/1‑150 (MQ=255)        gCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCa                                                                                                                                              >  1:966027/1‑150 (MQ=255)         ccACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAt                                                                                                                                             <  1:116612/150‑1 (MQ=255)           aCCACAATATGATTCATCGCGGGAATAACAGGTAAATCAGAAGGATCATAGCCAATTCCAGAAGATGTCTGCAAATTAACCTCATATTGTTGCGCCCAACTATGGCATTCATGCGTTCCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑gt                                                                                                                                          <  2:977610/150‑2 (MQ=255)            ccccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGa                                                                                                                                         <  2:788330/150‑1 (MQ=255)                aaTATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGGCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCAT‑GGAGCaa                                                                                                                                     <  1:661667/150‑1 (MQ=255)                  tatGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCAT‑GGAGCAACg                                                                                                                                   <  2:386974/150‑1 (MQ=255)                    tGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTa                                                                                                                                 <  1:640753/150‑1 (MQ=255)                    tGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTa                                                                                                                                 <  2:964585/150‑1 (MQ=255)                      aTTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑TGAGCAACGTACa                                                                                                                               >  2:1112694/1‑150 (MQ=255)                        tCATCGCGGCAATAACAGCTAAATCAGAACGATCAAATCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGt                                                                                                                             <  2:1284113/150‑1 (MQ=255)                        tCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGGTACAATCCGGTATGTGAA‑GGAGCAACCTaaact                                                                                                                             >  2:1487709/1‑146 (MQ=255)                        tCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGt                                                                                                                             <  1:286642/150‑1 (MQ=255)                         cATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTg                                                                                                                            <  2:624189/150‑1 (MQ=255)                         cATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTg                                                                                                                            <  1:1446786/150‑1 (MQ=255)                           tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGac                                                                                                                          <  1:1000198/150‑1 (MQ=255)                            cgcgGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTCCATGGCGGTGTGTGAT‑GGGGCAACGTACAGTGgca                                                                                                                         >  1:981407/1‑150 (MQ=255)                             gcgGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGacac                                                                                                                        >  2:36548/1‑150 (MQ=255)                             gcgGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCGCCAAGCCTGTCTGTCAT‑GGAGCAACGTACAGTGccac                                                                                                                        >  2:1047620/1‑150 (MQ=255)                                     aaCAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAAt                                                                                                               <  1:1037839/150‑1 (MQ=255)                                      aCAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAAtt                                                                                                              <  2:836455/150‑1 (MQ=255)                                       cAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTa                                                                                                             <  1:1665445/150‑1 (MQ=255)                                            aaaTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCgtg                                                                                                        >  2:683892/1‑150 (MQ=255)                                              aTCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTGTACAATCCTGTGTGTGATGGGGGCAAC‑TAAAGTGACACAGGGTTAATTaactggt                                                                                                       >  1:1271853/1‑144 (MQ=255)                                                    aCGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTaa                                                                                                >  1:1195276/1‑150 (MQ=255)                                                        tCAAAGCCAATTCCCGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGGCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCa                                                                                            <  1:1439917/150‑1 (MQ=255)                                                            aGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTc                                                                                        <  1:186888/150‑1 (MQ=255)                                                              ccAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCtt                                                                                      <  1:683892/150‑1 (MQ=255)                                                                aaTTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTACATGCCGGTCTGTAAT‑GGGGCAACCTACAGTGACACA‑CGTTAATTACCGGGGTGGTTAAAGCCTGTCGTTg                                                                                    >  1:283177/1‑150 (MQ=255)                                                                     cAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCAAGCCAGTGTGTGAA‑GGAGCAACGTTCAGTGACACA‑GGTTAATTAACGGGTTGGGTAAAGCATGTGTTTGAGCgg                                                                               >  2:470380/1‑149 (MQ=255)                                                                         agaTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTCCAAGCCGGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTTATTAACGTGTTGGTTAAAGCATGTGTTTGAGCGTAtta                                                                           >  2:1336622/1‑149 (MQ=255)                                                                          gaTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTcc                                                                          <  2:1276434/150‑1 (MQ=255)                                                                              tCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGc                                                                      <  2:948645/150‑1 (MQ=255)                                                                                 gCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACGGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTAc                                                                   <  2:471049/150‑1 (MQ=255)                                                                                         aCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGc                                                           <  1:1422443/150‑1 (MQ=255)                                                                                             cATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GCTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGa                                                       <  2:786867/150‑1 (MQ=255)                                                                                             aaTAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGa                                                       <  1:1381760/149‑1 (MQ=255)                                                                                                   gTTGCGCCGAACTATGCCATTCATCCGTAGCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttctt                                                 <  1:1187853/150‑1 (MQ=255)                                                                                                         cccAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatcca                                           <  2:961805/150‑1 (MQ=255)                                                                                                                  gCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGt                                  <  1:1487709/150‑1 (MQ=255)                                                                                                                   ccATTCATCCGTACCGGAGTCGTTTCTACAATCCAGTATGTTAT‑GGAGCAACGTACAATTAAAAA‑CGTTAATTACCGTGTTGGTTTAAGAATGTATTTGAGCGTATTACGGGCGACAAAAAGGCAGGATTTATTAATACAAAAGaaaatt                                 >  2:654235/1‑147 (MQ=255)                                                                                                                     aTTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgat                               <  1:1047058/150‑1 (MQ=255)                                                                                                                       tCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtgatga                             <  1:36548/150‑1 (MQ=255)                                                                                                                            cGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtgag                        <  1:1336622/150‑1 (MQ=255)                                                                                                                                cccgactCGTTGCTCCATGCCTGTGTGTCAT‑GGAGCAAAGTAGAGTTACACA‑GGTTAATTACCGTGTTTGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATAATGGCAGGAGTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgg                    <  2:1245168/144‑2 (MQ=255)                                                                                                                                ccGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgc                    <  1:508988/150‑1 (MQ=255)                                                                                                                                  ggAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgcgc                  <  1:61537/150‑1 (MQ=255)                                                                                                                                  ggAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgcgc                  <  2:1195276/150‑1 (MQ=255)                                                                                                                                     gTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGa               >  1:821031/1‑150 (MQ=255)                                                                                                                                         ttGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACggn           <  1:654235/150‑2 (MQ=255)                                                                                                                                         ttGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgc           <  1:1217332/150‑1 (MQ=255)                                                                                                                                          tGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg          <  1:1288828/150‑1 (MQ=255)                                                                                                                                           gCTCCATGACAGACAGGCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc         >  2:1552216/1‑150 (MQ=255)                                                                                                                                           gCTCCATGACAGACAGGCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc         >  2:673262/1‑150 (MQ=255)                                                                                                                                           gCTCCATGACAGACAGGCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc         >  1:1353226/1‑150 (MQ=255)                                                                                                                                           gCTCCATGACAGACAGGCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc         >  1:860627/1‑150 (MQ=255)                                                                                                                                            cTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca        >  2:484622/1‑150 (MQ=255)                                                                                                                                              ccATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc      >  2:225453/1‑150 (MQ=255)                                                                                                                                              ccATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc      >  1:1014146/1‑150 (MQ=255)                                                                                                                                                aTGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactc    <  1:995573/150‑5 (MQ=255)                                                                                                                                                 tGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc   >  1:1413854/1‑150 (MQ=255)                                                                                                                                                 tGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc   <  1:219472/150‑1 (MQ=255)                                                                                                                                                 tGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc   >  2:403244/1‑150 (MQ=255)                                                                                                                                                  gCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca  <  1:225453/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCA  >  NC_000913/2476868‑2477162

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