Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,430,519	T→C	S491P (TCA→CCA)	stfR →	Rac prophage; putative tail fiber protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,430,519	0	T	C	100.0%	228.1 / NA	64	S491P (TCA→CCA)	stfR	Rac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/35);  total (29/35)

CGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTAAC  >  NC_000913/1430371‑1430665                                                                                                                                                     |                                                                                                                                                   cgcCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCa                                                                                                                                                  <  2:688731/151‑1 (MQ=255)  gcCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAg                                                                                                                                                 <  1:63496/151‑1 (MQ=255)   ccGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGc                                                                                                                                                >  2:505783/1‑151 (MQ=255)    cGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGcc                                                                                                                                               >  2:609736/1‑151 (MQ=255)        aCAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATg                                                                                                                                           <  1:360028/151‑1 (MQ=255)         cagcagCTAGACGGGCAGAGGATATTGCATCCGCCGTGGCTCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTa                                                                                                                                             <  2:563117/148‑1 (MQ=255)               cTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATg                                                                                                                                     <  1:44223/150‑1 (MQ=255)                 aaaCGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCa                                                                                                                                   <  1:187026/150‑1 (MQ=255)                  aaCGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGc                                                                                                                                    >  2:639747/1‑148 (MQ=255)                    cGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACG‑AAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAgaga                                                                                                                               >  2:136932/1‑150 (MQ=255)                      ggCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGaaa                                                                                                                             <  1:271127/151‑1 (MQ=255)                        cAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACg                                                                                                                           >  2:99308/1‑151 (MQ=255)                         agagGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGt                                                                                                                          <  2:312394/151‑1 (MQ=255)                           agGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCt                                                                                                                        <  2:375684/151‑1 (MQ=255)                            ggATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCt                                                                                                                        <  1:357638/150‑1 (MQ=255)                             gATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGc                                                                                                                      <  1:357135/151‑1 (MQ=255)                                  tGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGaa                                                                                                                  <  1:60970/150‑1 (MQ=255)                                   gCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGaaa                                                                                                                 <  2:385385/150‑1 (MQ=255)                                   gCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGaa                                                                                                                  <  2:289351/149‑1 (MQ=255)                                   gCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGaa                                                                                                                  >  2:279819/1‑149 (MQ=255)                                   gCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACG‑AAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGaaa                                                                                                                 <  1:136899/149‑1 (MQ=255)                                    cATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGaa                                                                                                                  >  1:437737/1‑148 (MQ=255)                                     aTCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGAc                                                                                                              <  1:115442/151‑1 (MQ=255)                                     aTCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGAc                                                                                                              >  2:268106/1‑151 (MQ=255)                                       ccgccgTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGAcc                                                                                                             <  1:562715/150‑1 (MQ=255)                                        cgccgTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAg                                                                                                           <  2:563119/151‑1 (MQ=255)                                           cgTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAAc                                                                                                        <  2:180064/151‑1 (MQ=255)                                            gTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAAc                                                                                                        <  1:208851/150‑1 (MQ=255)                                                     gAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGat                                                                                               >  1:294001/1‑150 (MQ=255)                                                       ggATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGAtata                                                                                            <  1:636625/151‑1 (MQ=255)                                                        gATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATAc                                                                                           <  2:559386/151‑1 (MQ=255)                                                         aTGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATAcc                                                                                          <  2:432039/151‑1 (MQ=255)                                                           gCAAGTACGACGAAAAAGGGGATAGTACAGCTCATCAGGGCGACCAACAGTACGCCTGAAAAGCTAGCGGCAACGCCAAAGGCAGTAAACCCAACCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGCCGCGGATATACCCg                                                                                        >  1:112210/1‑151 (MQ=255)                                                             aaGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGa                                                                                       <  1:666208/150‑1 (MQ=255)                                                               gTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATa                                                                                     >  1:603906/1‑150 (MQ=255)                                                                 acgacgAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAgg                                                                                  <  1:179056/151‑1 (MQ=255)                                                                 acgacgAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAg                                                                                   <  2:278935/150‑1 (MQ=255)                                                                          aaGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTcc                                                                         >  2:279818/1‑151 (MQ=255)                                                                          aaGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTcc                                                                         >  2:304057/1‑151 (MQ=255)                                                                                       aGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTa                                                             <  2:466520/150‑1 (MQ=255)                                                                                       aGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTa                                                             >  2:496930/1‑150 (MQ=255)                                                                                        gCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAc                                                           <  2:310433/151‑1 (MQ=255)                                                                                          tCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAcgc                                                         >  2:113978/1‑151 (MQ=255)                                                                                           cagcagTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAcg                                                          <  1:371754/149‑1 (MQ=255)                                                                                           cagcagTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAAcg                                                          <  2:601346/149‑1 (MQ=255)                                                                                               agTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTc                                                     >  1:695384/1‑150 (MQ=255)                                                                                                    gACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTaa                                                >  2:119629/1‑150 (MQ=255)                                                                                                     aCCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTaaa                                               >  1:359858/1‑150 (MQ=255)                                                                                                      ccAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAAc                                             >  2:545389/1‑151 (MQ=255)                                                                                                      ccAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAAc                                             >  2:635339/1‑151 (MQ=255)                                                                                                        aaCAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAAc                                             >  1:175709/1‑149 (MQ=255)                                                                                                         aCAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAAc                                             >  2:670777/1‑148 (MQ=255)                                                                                                           aGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACt                                        <  1:312666/151‑1 (MQ=255)                                                                                                                gTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCt                                   <  1:416048/151‑1 (MQ=255)                                                                                                                        cGCTGGCGGCAACGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCg                           <  1:474141/151‑1 (MQ=255)                                                                                                                                  aaCGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGt                 >  2:567711/1‑151 (MQ=255)                                                                                                                                  aaCGCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGt                 >  2:573210/1‑151 (MQ=255)                                                                                                                                     gCCAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTAt              >  1:3966/1‑151 (MQ=255)                                                                                                                                       cAAAGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTAtgt            >  1:325301/1‑151 (MQ=255)                                                                                                                                         aaGGCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGc          >  1:449963/1‑151 (MQ=255)                                                                                                                                            gCAGTAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCggg       >  1:552493/1‑151 (MQ=255)                                                                                                                                                tAAAACCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTa    >  2:729141/1‑150 (MQ=255)                                                                                                                                                 aaaaCCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTAAc  <  2:385564/151‑1 (MQ=255)                                                                                                                                                 aaaaCCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTAAc  <  2:64698/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CGCCGAGACAGCAGCTAAACGGGCAGAGGATATTGCATCCGCCGTGGCGCTTGAGGATGCAAGTACGACGAAAAAGGGGATAGTACAGCTCAGCAGTGCGACCAACAGTACGTCTGAAACGCTGGCGGCAACGCCAAAGGCAGTAAAATCAGCCTATGACAATGCAGAGAAACGTCTGCAGAAAGACCAGAACGGCGCTGATATACCCGATAAGGGATGCTTCCTGAACAACATTAACGCGGTCAGTAAAACAGACTTTGCTGATAAGCGTGGTATGCGTTATGTGCGGGTTAAC  >  NC_000913/1430371‑1430665

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