Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,477,022	T→G	9.4%	D203A (GAC→GCC)	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,022	0	T	G	9.4%	162.2 / 9.5	63	D203A (GAC→GCC)	dsdC	DNA‑binding transcriptional dual regulator DsdC
Reads supporting (aligned to +/- strand):  ref base T (19/38);  new base G (4/2);  total (23/40)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGAT  >  NC_000913/2476875‑2477170                                                                                                                                                    |                                                                                                                                                       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                                                                                                                                         <  2:964585/150‑1 (MQ=255)             tGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTa                                                                                                                                         <  1:640753/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                                                                                                                                    <  1:1446786/150‑1 (MQ=255)                  cATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTg                                                                                                                                    <  2:624189/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                 >  1:1353226/1‑150 (MQ=255)                                                                                                                                    gCTCCATGACAGACAGGCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                 >  2:673262/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              >  1:1014146/1‑150 (MQ=255)                                                                                                                                       ccATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc              >  2:225453/1‑150 (MQ=255)                                                                                                                                         aTGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactc            <  1:995573/150‑5 (MQ=255)                                                                                                                                          tGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc           >  2:403244/1‑150 (MQ=255)                                                                                                                                          tGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc           <  1:219472/150‑1 (MQ=255)                                                                                                                                          tGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc           >  1:1413854/1‑150 (MQ=255)                                                                                                                                           gCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca          <  1:225453/150‑1 (MQ=255)                                                                                                                                            ccgttctgtcTT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCTTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa         <  1:580892/146‑1 (MQ=255)                                                                                                                                                 ctgtcAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAg    <  1:637688/150‑1 (MQ=255)                                                                                                                                                 ctgtcAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAg    >  1:883531/1‑150 (MQ=255)                                                                                                                                                   gtcAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGAGCTGACGGCGca                >  1:1348418/1‑136 (MQ=255)                                                                                                                                                   gtcAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGAGCTGACGGCGca                <  2:1348418/136‑1 (MQ=255)                                                                                                                                                   gtcAT‑GGAGCAACGTACAGTGACACA‑GGATAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTACGGGCTACATAATGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGAt  >  2:471949/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGAT  >  NC_000913/2476875‑2477170

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: ‑