Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	587,023	A→C	intergenic (‑115/+68)	envY ← / ← ybcH	porin thermoregulatory transcriptional activator/PRK09936 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	587,023	0	A	C	100.0%	378.4 / NA	96	intergenic (‑115/+68)	envY/ybcH	porin thermoregulatory transcriptional activator/PRK09936 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (51/45);  total (51/45)

CCACGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACG  >  NC_000913/586878‑587168                                                                                                                                                  |                                                                                                                                                  ccacGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTg                                                                                                                                              <  2:431604/151‑1 (MQ=255)    cGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATg                                                                                                                                           <  1:394146/151‑1 (MQ=255)    cGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCGTTAGCACGCTTTTTAGCAATTACTCTTGATg                                                                                                                                           <  2:213072/151‑1 (MQ=255)       aaGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCCTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGaa                                                                                                                                        <  1:943100/151‑1 (MQ=255)        aGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAAt                                                                                                                                       <  2:570879/151‑1 (MQ=255)        aGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAAt                                                                                                                                       >  1:707068/1‑151 (MQ=255)         ggTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAAtt                                                                                                                                      >  2:948995/1‑151 (MQ=255)          gTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTc                                                                                                                                     >  2:627203/1‑151 (MQ=255)             cACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtga                                                                                                                                  <  1:303648/151‑1 (MQ=255)             cACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtg                                                                                                                                   >  2:870104/1‑150 (MQ=255)             cACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtg                                                                                                                                   >  1:385840/1‑150 (MQ=255)               ctgctgCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtgat                                                                                                                                 <  1:798135/150‑1 (MQ=255)                tgctgcTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtgatg                                                                                                                                >  1:165941/1‑150 (MQ=255)                 gctgctCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATgag                                                                                                                              <  2:644553/151‑1 (MQ=255)                  ctgctCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGTCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCtgatg                                                                                                                                >  2:1188688/1‑148 (MQ=255)                    gctCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGagag                                                                                                                            <  2:1014034/150‑1 (MQ=255)                      tCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCg                                                                                                                          >  1:140392/1‑150 (MQ=255)                       cAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCg                                                                                                                          <  1:436293/149‑1 (MQ=255)                        aaTTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGaaa                                                                                                                       <  1:806512/151‑1 (MQ=255)                            gCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAagag                                                                                                                    <  1:323955/150‑1 (MQ=255)                             cATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGC‑TTTTAGCAATTACTCTTGATGGAATTTTGATGAGAGCGAAAGAgg                                                                                                                   >  1:216554/1‑149 (MQ=255)                              aTTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGt                                                                                                                  >  2:559167/1‑150 (MQ=255)                                  cGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGcc                                                                                                             <  2:305554/151‑1 (MQ=255)                                     acacctcAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCCGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTTATGAGAGCGAAAGAGGTAAGCCAgg                                                                                                          <  2:173003/148‑1 (MQ=255)                                     aCTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAgg                                                                                                          >  1:151286/1‑151 (MQ=255)                                        cctcAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCg                                                                                                       <  1:953772/151‑1 (MQ=255)                                          tcAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGt                                                                                                      >  2:760686/1‑150 (MQ=255)                                             gATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTAccc                                                                                                  <  1:628244/151‑1 (MQ=255)                                              atatCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCg                                                                                                 >  2:874728/1‑151 (MQ=255)                                               tatCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCg                                                                                                 >  1:888391/1‑150 (MQ=255)                                                atCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGAc                                                                                               <  1:440760/151‑1 (MQ=255)                                                atCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGAc                                                                                               <  1:700117/151‑1 (MQ=255)                                                     aaaCTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTAcc                                                                                          <  2:1155847/151‑1 (MQ=255)                                                        cTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTgg                                                                                       >  1:412536/1‑151 (MQ=255)                                                          ccGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTgg                                                                                       >  2:619212/1‑149 (MQ=255)                                                          ccGCGCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAAAGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGCATTCTGATGAGAGCGAGATAGGTTAGCCAGCTCGTACCCGACTTACCTg                                                                                        >  1:1057866/1‑148 (MQ=255)                                                            gCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTggagg                                                                                    >  2:994686/1‑150 (MQ=255)                                                               cAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTggagg                                                                                    >  2:408523/1‑147 (MQ=255)                                                                aaaGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGAt                                                                                <  2:791680/150‑1 (MQ=255)                                                                       cTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATAc                                                                         >  1:460438/1‑150 (MQ=255)                                                                         aTGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAAt                                                                           <  2:499491/146‑1 (MQ=255)                                                                         aTGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTc                                                                       <  1:1269909/150‑1 (MQ=255)                                                                         aTGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTc                                                                       <  2:730401/150‑1 (MQ=255)                                                                         aTGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCg                                                                      >  2:920064/1‑151 (MQ=255)                                                                             ttCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCgaga                                                                   >  1:347588/1‑150 (MQ=255)                                                                                cTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGc                                                               >  1:1264501/1‑151 (MQ=255)                                                                                  gCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGc                                                               <  1:656896/149‑1 (MQ=255)                                                                                    aTGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTg                                                           <  1:861937/151‑1 (MQ=255)                                                                                     tGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTgc                                                          >  2:1261221/1‑151 (MQ=255)                                                                                      gAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTgcg                                                         <  1:1169976/151‑1 (MQ=255)                                                                                          gATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGAc                                                     >  1:1238778/1‑151 (MQ=255)                                                                                             cGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACgg                                                   <  2:390208/150‑1 (MQ=255)                                                                                                  cGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAAt                                             >  1:512126/1‑151 (MQ=255)                                                                                                   gTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATa                                            <  2:1048851/151‑1 (MQ=255)                                                                                                     tCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAg                                           <  2:1155060/150‑1 (MQ=255)                                                                                                      cGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAgcg                                         >  2:117417/1‑151 (MQ=255)                                                                                                      cGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAgcg                                         >  1:1272023/1‑151 (MQ=255)                                                                                                      cGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAgcg                                         <  1:578999/151‑1 (MQ=255)                                                                                                       gCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAgcgc                                        <  2:284112/151‑1 (MQ=255)                                                                                                        ccGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAgcgc                                        >  2:209394/1‑150 (MQ=255)                                                                                                         cGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAgcgc                                        >  1:884404/1‑149 (MQ=255)                                                                                                          gATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCag                                      <  1:1270357/150‑1 (MQ=255)                                                                                                           aTAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCag                                      >  1:993630/1‑149 (MQ=255)                                                                                                            tAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCag                                      >  1:1261475/1‑148 (MQ=255)                                                                                                            tAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCag                                      >  2:127736/1‑148 (MQ=255)                                                                                                             aaCGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCagaga                                   <  1:805721/150‑1 (MQ=255)                                                                                                             aaCGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCagag                                    >  1:488471/1‑149 (MQ=255)                                                                                                              aCGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCagag                                    >  2:485035/1‑148 (MQ=255)                                                                                                                    ttctttctTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGaaa                               >  1:979949/1‑147 (MQ=255)                                                                                                                    ttctttctTTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAAt                              >  1:307745/1‑148 (MQ=255)                                                                                                                     tctttcttTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGaaa                               <  2:923751/146‑1 (MQ=255)                                                                                                                     tctttcttTAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGaa                                <  1:577051/145‑1 (MQ=255)                                                                                                                         tctttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAgtg                      >  1:72070/1‑151 (MQ=255)                                                                                                                         tctttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAgtg                      >  1:160090/1‑151 (MQ=255)                                                                                                                         tctttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAgt                       <  1:980933/150‑1 (MQ=255)                                                                                                                          ctttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAgcgc                                        <  1:973574/132‑1 (MQ=255)                                                                                                                          ctttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAgcgc                                        >  2:973633/1‑132 (MQ=255)                                                                                                                            ttAGCGCGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTGCCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAgtg                      <  1:184897/148‑1 (MQ=255)                                                                                                                            ttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAgtg                      <  1:30900/148‑1 (MQ=255)                                                                                                                            ttAGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCt                   >  1:625373/1‑151 (MQ=255)                                                                                                                              aGCACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCttt                 <  1:513135/151‑1 (MQ=255)                                                                                                                                cACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTAc               >  1:258037/1‑151 (MQ=255)                                                                                                                                cACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTAc               >  2:664753/1‑151 (MQ=255)                                                                                                                                cACGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTAc               <  1:1028377/151‑1 (MQ=255)                                                                                                                                  cGCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCg             <  2:782517/151‑1 (MQ=255)                                                                                                                                   gCTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGc            >  1:283077/1‑151 (MQ=255)                                                                                                                                    cTTTTTAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGc            >  1:408980/1‑150 (MQ=255)                                                                                                                                     tttttAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCa           >  1:914430/1‑150 (MQ=255)                                                                                                                                      ttttAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCa           >  2:1205318/1‑149 (MQ=255)                                                                                                                                        ttAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAag         >  1:135668/1‑149 (MQ=255)                                                                                                                                         tAGCAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAgg      >  2:779387/1‑151 (MQ=255)                                                                                                                                         tAGCAATTACTCTGGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAagag       <  1:430223/150‑1 (MQ=255)                                                                                                                                            cAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAAc   <  1:251103/151‑1 (MQ=255)                                                                                                                                            cAATTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAAc   >  1:442177/1‑151 (MQ=255)                                                                                                                                             aaTTACTCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACg  >  2:1002913/1‑151 (MQ=255)                                                                                                                                               ttACTCTTGATGGAATTCTGATAAGAGCGAAAGAAGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAAc   <  1:990395/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CCACGCAAGGTTCACTGCTGCTCAATTGCATTTCGCACTCCTCAGATATCAGAAACTCCGCTCAAAGGATCTATGCTTCCTGCATGAGTGATCGGCCCGTTCGCCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACG  >  NC_000913/586878‑587168

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