Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,776,212	T→C	G98G (GGT→GGC)	yfjX →	CP4‑57 prophage; putative antirestriction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,776,212	0	T	C	100.0%	271.9 / NA	76	G98G (GGT→GGC)	yfjX	CP4‑57 prophage; putative antirestriction protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/35);  total (41/35)

GCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCC  >  NC_000913/2776064‑2776356                                                                                                                                                     |                                                                                                                                                 gcATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGgcgc                                                                                                                                                >  1:422088/1‑151 (MQ=255)   aTTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCgg                                                                                                                                              >  2:796132/1‑151 (MQ=255)      ttCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGaa                                                                                                                                            >  1:822077/1‑150 (MQ=255)         gCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGa                                                                                                                                        <  2:120408/151‑1 (MQ=255)               atggaCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCgg                                                                                                                                  <  1:195608/151‑1 (MQ=255)                  gaCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAg                                                                                                                               <  2:332510/151‑1 (MQ=255)                     cGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCg                                                                                                                            >  2:672028/1‑151 (MQ=255)                        ttCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGc                                                                                                                           >  2:411428/1‑149 (MQ=255)                        ttCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGc                                                                                                                           >  1:311823/1‑149 (MQ=255)                         tCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAg                                                                                                                         <  1:822213/150‑1 (MQ=255)                            gtgATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAgg                                                                                                                        <  1:1152984/148‑1 (MQ=255)                                 gaGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTg                                                                                                                <  1:839087/151‑1 (MQ=255)                                   gTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTCATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTc                                                                                                              <  2:954258/151‑1 (MQ=255)                                     aCTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCt                                                                                                             <  2:275938/150‑1 (MQ=255)                                       tGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGcc                                                                                                          >  1:687500/1‑151 (MQ=255)                                       tGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGcc                                                                                                          >  2:801802/1‑151 (MQ=255)                                           ggggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCCGCCTggt                                                                                                      >  1:182754/1‑149 (MQ=255)                                            gggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAtt                                                                                                     >  1:992238/1‑151 (MQ=255)                                             ggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTg                                                                                                    >  1:619616/1‑151 (MQ=255)                                               tATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGc                                                                                                                           >  2:163642/1‑126 (MQ=255)                                               tATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGc                                                                                                                           <  1:163640/126‑1 (MQ=255)                                                    ggTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAatat                                                                                               <  1:673739/149‑1 (MQ=255)                                                    ggTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                                                                                                <  1:1025116/148‑1 (MQ=255)                                                           ttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGGATCGCTGTCTGCCTGATTGAATATAGCCATCa                                                                                       >  2:324029/1‑150 (MQ=255)                                                              cacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcc                                                                                   >  2:1182892/1‑151 (MQ=255)                                                               acGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct                                                                                  <  2:1089169/151‑1 (MQ=255)                                                                cGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct                                                                                  <  1:306253/150‑1 (MQ=255)                                                                  cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgcc                                                                               >  1:573604/1‑151 (MQ=255)                                                                   tCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgcct                                                                              >  2:1123401/1‑151 (MQ=255)                                                                    cAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgcct                                                                              >  1:137450/1‑150 (MQ=255)                                                                    cAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgcct                                                                              >  1:245104/1‑150 (MQ=255)                                                                     aGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTc                                                                             >  1:758979/1‑150 (MQ=255)                                                                      gCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCa                                                                            >  1:728417/1‑150 (MQ=255)                                                                        aaTGGCGGAGCGTTTATGGCTCCTGATGCTGGCGGTGACGGTAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGGTGTCTGCCTGATTGAATATAGTCATCACGCCTGCCTCAcc                                                                          >  2:505365/1‑150 (MQ=255)                                                                        aaTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCg                                                                         <  2:773310/151‑1 (MQ=255)                                                                             cGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGc                                                                    >  2:275521/1‑151 (MQ=255)                                                                             cGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGc                                                                    <  2:758990/151‑1 (MQ=255)                                                                              ggAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCg                                                                   >  2:553118/1‑151 (MQ=255)                                                                                aGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGAt                                                                 >  1:826545/1‑151 (MQ=255)                                                                                 gCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGAt                                                                 <  2:573328/150‑1 (MQ=255)                                                                                  cGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGc                                                               >  1:31113/1‑151 (MQ=255)                                                                                     ttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCAt                                                            >  1:425071/1‑151 (MQ=255)                                                                                       aTGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGa                                                          <  2:569881/151‑1 (MQ=255)                                                                                         ggCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACg                                                        <  1:949158/151‑1 (MQ=255)                                                                                          gCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACgg                                                       >  2:476675/1‑151 (MQ=255)                                                                                             ccTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGa                                                      <  2:683469/149‑1 (MQ=255)                                                                                              cTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGc                                                    >  2:216063/1‑150 (MQ=255)                                                                                                gATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACt                                                 <  1:1163548/151‑1 (MQ=255)                                                                                                 aTGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGc                                                    >  1:764227/1‑147 (MQ=255)                                                                                                 aTGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACt                                                 >  2:629179/1‑150 (MQ=255)                                                                                                 aTGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGGATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACt                                                 >  1:761262/1‑150 (MQ=255)                                                                                                    cTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTAtt                                              >  1:1086556/1‑150 (MQ=255)                                                                                                       aCGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCg                                          <  2:985320/151‑1 (MQ=255)                                                                                                            gACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGc                                     <  1:178920/151‑1 (MQ=255)                                                                                                              cGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCg                                    <  2:1012256/150‑1 (MQ=255)                                                                                                                 tAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGAt                                >  1:1065058/1‑151 (MQ=255)                                                                                                                  aaaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGAtt                               <  2:724108/151‑1 (MQ=255)                                                                                                                  aaaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGAt                                <  2:1049867/150‑1 (MQ=255)                                                                                                                  aaaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGAt                                <  1:646865/150‑1 (MQ=255)                                                                                                                      ggCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGc                           >  2:723776/1‑151 (MQ=255)                                                                                                                         aTTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGCCGGAGCACTATTACCGCCTGCGGGATTACGctct                        <  2:733841/151‑1 (MQ=255)                                                                                                                           ttGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTgc                      >  1:460322/1‑151 (MQ=255)                                                                                                                             gCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTgc                      >  1:685755/1‑149 (MQ=255)                                                                                                                               tCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTgcagc                   >  2:147783/1‑150 (MQ=255)                                                                                                                                cAACGGCACGATTGGCAATGGCGCGGAAATGAGCCCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATGGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTgcagc                   <  1:838967/149‑1 (MQ=255)                                                                                                                                 accGGCATGAAGGGCAATGGCGCGGAAATGAGGGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTCAATATGGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCAc                 <  1:158965/148‑1 (MQ=38)                                                                                                                                 aaCGGCATGAATGGCAATGGCGCGGCAATGAGCGCGGAAGCCGCAGGGATCGCTGGCTGCCTTATTGAATATAACCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCAcc                >  1:137259/1‑151 (MQ=255)                                                                                                                                   cGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCAccc               >  1:1184979/1‑150 (MQ=255)                                                                                                                                   cGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCGCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTgcagc                   <  1:1039308/146‑1 (MQ=255)                                                                                                                                    ggCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCCCGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCt              <  1:1156778/150‑1 (MQ=255)                                                                                                                                    ggCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTg             <  1:561371/151‑1 (MQ=255)                                                                                                                                     gCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGa            <  2:902915/151‑1 (MQ=255)                                                                                                                                       aTGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAAt          >  2:638138/1‑151 (MQ=255)                                                                                                                                       aTGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAAt          >  2:604558/1‑151 (MQ=255)                                                                                                                                           aTGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAg      <  2:1093126/151‑1 (MQ=255)                                                                                                                                            tGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAgc     >  1:254636/1‑151 (MQ=255)                                                                                                                                                aaTGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGcc  <  1:508529/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCC  >  NC_000913/2776064‑2776356

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