Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,477,020	T→A	9.5%	R204* (AGA→TGA) ‡	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,020	0	T	A	9.5%	175.6 / 8.5	63	R204* (AGA→TGA) ‡	dsdC	DNA‑binding transcriptional dual regulator DsdC
Reads supporting (aligned to +/- strand):  ref base T (19/38);  new base A (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 = 5.69e-01

CATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAG  >  NC_000913/2476871‑2477168                                                                                                                                                      |                                                                                                                                                       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                                                                                                                                       <  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                                                                                                                                  <  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               >  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            >  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:219472/150‑1 (MQ=255)                                                                                                                                              tGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc         >  1:1413854/1‑150 (MQ=255)                                                                                                                                              tGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc         >  2:403244/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)                                                                                                                                                      |                                                                                                                                                       CATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAT‑GGAGCAACGTACAGTGACACA‑GGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAG  >  NC_000913/2476871‑2477168

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