Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	577,451	T→C	G18G (GGT→GGC)	essD →	DLP12 prophage; putative phage lysis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	577,451	0	T	C	100.0%	236.6 / NA	65	G18G (GGT→GGC)	essD	DLP12 prophage; putative phage lysis protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/29);  total (36/29)

ACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGTAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAAAGAAGATAAGCGCAAGG  >  NC_000913/577307‑577596                                                                                                                                                 |                                                                                                                                                  aCCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGcc                                                                                                                                             >  1:7064/1‑151 (MQ=255) aCCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGc                                                                                                                                              >  2:1523/1‑150 (MQ=255)    cAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTGGCCTATGGCACCTCAGCAGGCAGTGCCggg                                                                                                                                          <  1:404086/151‑1 (MQ=255)       tACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTAc                                                                                                                                       >  2:55886/1‑151 (MQ=255)          cTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTgg                                                                                                                                    >  1:120752/1‑151 (MQ=255)          cTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTgg                                                                                                                                    >  2:891404/1‑151 (MQ=255)            gAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTgg                                                                                                                                    <  1:443766/149‑1 (MQ=255)            gAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGtt                                                                                                                                  >  1:151179/1‑151 (MQ=255)            gAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGt                                                                                                                                   <  2:244825/150‑1 (MQ=255)             aGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGttt                                                                                                                                 >  1:210177/1‑151 (MQ=255)             aGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGttt                                                                                                                                 >  1:792401/1‑151 (MQ=255)                  ttCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTAc                                                                                                                             >  2:185760/1‑150 (MQ=255)                        gACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAgctgct                                                                                                                      <  2:324190/151‑1 (MQ=255)                        gACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAgctgct                                                                                                                      <  2:519569/151‑1 (MQ=255)                              gAGTCTGGTACGTTTCGTAGTATTCGCTCAATTTACACCCGCTTTGTCTGCGCGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATa                                                                                                                 >  2:669625/1‑150 (MQ=255)                               aGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATaa                                                                                                                >  1:724628/1‑150 (MQ=255)                                    gTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTc                                                                                                            <  2:186024/149‑1 (MQ=255)                                    gTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCAc                                                                                                          <  1:694090/151‑1 (MQ=255)                                      tACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTc                                                                                                            <  1:707481/147‑1 (MQ=255)                                      tACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTc                                                                                                        >  2:27273/1‑151 (MQ=255)                                        cGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTccc                                                                                                      >  1:742853/1‑151 (MQ=255)                                           ttCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTc                                                                                                    >  2:402751/1‑150 (MQ=255)                                           ttCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTc                                                                                                    >  2:438064/1‑150 (MQ=255)                                             cGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTcac                                                                                                  <  2:304583/150‑1 (MQ=255)                                                    ttCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGgc                                                                                           >  1:731805/1‑150 (MQ=255)                                                        cTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGgcagca                                                                                       >  2:95192/1‑150 (MQ=255)                                                           aaTTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCTTGGATCAGTTAACAACGGGAGTCGCCTATGGCCCCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAAt                                                                                     >  1:42218/1‑149 (MQ=255)                                                            aTTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAgg                                                                                  >  2:783786/1‑151 (MQ=255)                                                                acacCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGtgtg                                                                              >  2:186176/1‑151 (MQ=255)                                                                 cacCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGtgtg                                                                              >  1:365186/1‑150 (MQ=255)                                                                    ccGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTgg                                                                          <  1:771234/151‑1 (MQ=255)                                                                       cTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGt                                                                        >  2:805110/1‑150 (MQ=255)                                                                         ttGTCTGCGAGGTTGGGTTATGAAATCCATGGATAAGTTAACAACGGGTTTCGCCTATGGCCCCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTTGGCAGCAATAGGTGTGCTGGGTa                                                                       <  2:463685/149‑1 (MQ=255)                                                                           gTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCt                                                                   <  2:671784/151‑1 (MQ=255)                                                                                  aGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTAttt                                                            <  1:438111/151‑1 (MQ=255)                                                                                  aGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTAttt                                                            >  2:753248/1‑151 (MQ=255)                                                                                      ggggTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGcc                                                        <  2:867192/151‑1 (MQ=255)                                                                                       gggTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCct                                                       >  1:133739/1‑151 (MQ=255)                                                                                       gggTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCct                                                       <  1:647187/151‑1 (MQ=255)                                                                                          ttATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCctgct                                                    >  1:545129/1‑151 (MQ=255)                                                                                                aaTCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGt                                               >  2:802275/1‑150 (MQ=255)                                                                                                aaTCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTa                                              <  2:530602/151‑1 (MQ=255)                                                                                                  tCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTAcc                                            >  2:748903/1‑151 (MQ=255)                                                                                                   ccATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCt                                           >  2:729758/1‑151 (MQ=255)                                                                                                      tGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGa                                         <  2:678093/150‑1 (MQ=255)                                                                                                       ggATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACa                                       <  1:272410/151‑1 (MQ=255)                                                                                                          tAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAAc                                    >  2:802118/1‑151 (MQ=255)                                                                                                           aaGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAAcc                                   >  2:256990/1‑151 (MQ=255)                                                                                                            aGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTGCCTGTCAAACCt                                  <  2:406744/151‑1 (MQ=255)                                                                                                             gTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCtt                                 <  1:905127/151‑1 (MQ=255)                                                                                                             gTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCtt                                 >  1:520072/1‑151 (MQ=255)                                                                                                             gTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCtt                                 <  1:596960/151‑1 (MQ=255)                                                                                                              ttAACAACGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCttt                                >  2:798266/1‑151 (MQ=255)                                                                                                                  caacGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTAttt                            <  1:69798/151‑1 (MQ=255)                                                                                                                   aacGGGTGTCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTc                           >  1:575463/1‑151 (MQ=255)                                                                                                                          gtCGCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGAtt                     >  2:430790/1‑150 (MQ=255)                                                                                                                             gCCTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAaag                 <  2:99239/151‑1 (MQ=255)                                                                                                                              ccTATGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAaaga                <  2:300492/151‑1 (MQ=255)                                                                                                                                 aTGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAaagaa               <  2:526522/149‑1 (MQ=255)                                                                                                                                 aTGGCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAaaga                >  2:644294/1‑148 (MQ=255)                                                                                                                                   ggCACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAAAGAAGATa           <  2:231190/151‑1 (MQ=255)                                                                                                                                     cACCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAAAGAAGATAAg         <  1:685566/151‑1 (MQ=255)                                                                                                                                      aCCTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAAAGAAGATAAgc        <  2:848429/151‑1 (MQ=255)                                                                                                                                        cTCAGCAGGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAAAGAAGATAAgcgc      >  2:383277/1‑151 (MQ=255)                                                                                                                                           agcagGCAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATATGTGTGTTGGGTAGCATGGTATATGGACTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAAAGAAGATAAGCGCAAg   <  2:132578/151‑1 (MQ=25)                                                                                                                                            gcaggcagTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAAAGAAGATAAGCGCAAgg  <  1:810350/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  ACCCAATACCTGAGTTATTCACTGACTCCGAGTCTGTTACGTTTCGTAGTATTCCCTCAATTTACACCCGCTTTGTCTGCGAGGTGGGGTTATGAAATCCATGGATAAGTTAACAACGGGTGTCGCCTATGGCACCTCAGCAGGTAGTGCCGGGTACTGGTTTTTACAGCTGCTAGATAAAGTCACTCCCTCACAGTGGGCAGCAATAGGTGTGCTGGGTAGCCTGGTATTTGGCCTGCTGACGTACCTGACAAACCTTTATTTCAAGATTAAAGAAGATAAGCGCAAGG  >  NC_000913/577307‑577596

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