Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,150,087	A→G	V60A (GTA→GCA)	eptC ←	LPS heptose I phosphoethanolamine transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,150,087	0	A	G	100.0%	249.5 / NA	66	V60A (GTA→GCA)	eptC	LPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/34);  total (32/34)

TGTTGGTTTCGAACATCACAAACAGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAA  >  NC_000913/4149941‑4150233                                                                                                                                                   |                                                                                                                                                   tgttggtTTCGAACATCACAAACAGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCgg                                                                                                                                                >  1:871307/1‑151 (MQ=255)           gAACATCACAAACAGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACcac                                                                                                                                      <  2:94639/151‑1 (MQ=255)           gAACATCACAAACAGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACcac                                                                                                                                      >  2:45128/1‑151 (MQ=255)            aaCATCACAAACAGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACcac                                                                                                                                      >  1:21377/1‑150 (MQ=255)            aaCATCACAAACAGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACcac                                                                                                                                      >  1:309576/1‑150 (MQ=255)                    aaaCAGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTg                                                                                                                             <  2:813598/151‑1 (MQ=255)                         gAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATa                                                                                                                         >  2:254116/1‑150 (MQ=255)                          aaCGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAAc                                                                                                                       >  1:22795/1‑151 (MQ=255)                          aaCGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAAc                                                                                                                       >  2:1041786/1‑151 (MQ=255)                             gCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGc                                                                                                                    <  2:904783/151‑1 (MQ=255)                                  gCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAg                                                                                                                 <  1:16777/149‑1 (MQ=255)                                  gCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTc                                                                                                               <  1:986169/151‑1 (MQ=255)                                   cGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTcg                                                                                                              <  2:859061/151‑1 (MQ=255)                                   cGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTcg                                                                                                              >  2:595800/1‑151 (MQ=255)                                       aaCTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGaa                                                                                                          <  1:101646/151‑1 (MQ=255)                                          tCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACt                                                                                                                              >  1:757579/1‑128 (MQ=255)                                          tCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACt                                                                                                                              <  2:757554/128‑1 (MQ=255)                                              gtccGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCAATTACTGCGGCAATAATTTTAATCCGCTTCGAAAAGAGGAATGCCGGGATCACCCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCg                                                                                                     <  2:657546/147‑1 (MQ=255)                                              gACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGtt                                                                                                   <  2:76980/151‑1 (MQ=255)                                                ccGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGtt                                                                                                   >  1:870899/1‑149 (MQ=255)                                                ccGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGtt                                                                                                   >  2:154879/1‑149 (MQ=255)                                                ccGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAg                                                                                                 <  2:948482/151‑1 (MQ=255)                                                  gTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTg                                                                                               <  2:635044/151‑1 (MQ=255)                                                    aGATGTCGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGcc                                                                                             <  2:593835/151‑1 (MQ=255)                                                        gACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACt                                                                                          <  1:24955/150‑1 (MQ=255)                                                        gACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACt                                                                                          <  1:1000282/150‑1 (MQ=255)                                                        gACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACt                                                                                          <  1:333617/150‑1 (MQ=255)                                                         aCGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACtat                                                                                        <  1:148498/151‑1 (MQ=255)                                                               tagCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCAc                                                                                  >  2:896289/1‑151 (MQ=255)                                                                     aGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGt                                                                            <  2:554319/151‑1 (MQ=255)                                                                      gcgcCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTa                                                                           <  1:486138/151‑1 (MQ=255)                                                                      gcgcCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTa                                                                           <  1:175505/151‑1 (MQ=255)                                                                         ccgccAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAaat                                                                        <  2:620175/151‑1 (MQ=255)                                                                          cgccAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAaat                                                                        >  2:1020468/1‑150 (MQ=255)                                                                             cAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATg                                                                    <  2:464115/151‑1 (MQ=255)                                                                             cAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATg                                                                    >  2:780948/1‑151 (MQ=255)                                                                                 gagGCCGCCCATAGCACCACGCCGATTACTTCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCt                                                                >  1:884533/1‑151 (MQ=255)                                                                                    gccgccCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTa                                                             <  2:766157/151‑1 (MQ=255)                                                                                               caccacGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTagag                                                  <  2:142681/151‑1 (MQ=255)                                                                                               caccacGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTag                                                    >  2:973500/1‑149 (MQ=255)                                                                                                accacGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTag                                                    >  1:840189/1‑148 (MQ=255)                                                                                                  cacGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGaaa                                               <  2:354407/151‑1 (MQ=255)                                                                                                  cacGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGaaa                                               <  2:607987/151‑1 (MQ=255)                                                                                                    cGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGaaaaa                                             <  1:161770/151‑1 (MQ=255)                                                                                                         aTTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCaaa                                        <  2:533185/151‑1 (MQ=255)                                                                                                          ttACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCaaa                                        >  2:689770/1‑150 (MQ=255)                                                                                                              tGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGt                                     >  2:812182/1‑149 (MQ=255)                                                                                                               gCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTa                                    <  1:134344/149‑1 (MQ=255)                                                                                                                cGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTaga                                                   <  2:318871/133‑1 (MQ=255)                                                                                                                cGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTaga                                                   >  1:318881/1‑133 (MQ=255)                                                                                                                cGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTagag                                 <  1:432684/151‑1 (MQ=255)                                                                                                                 ggCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGc                                <  1:119931/151‑1 (MQ=255)                                                                                                                 ggCAATAATTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGc                                >  1:275372/1‑151 (MQ=255)                                                                                                                    aataatTTTAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGt                             >  1:408259/1‑151 (MQ=255)                                                                                                                            tAATCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGACTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGAAGAGCAGTGCCCAAcc                     >  1:447519/1‑151 (MQ=255)                                                                                                                             aaTCCGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAAccc                    >  2:862767/1‑151 (MQ=255)                                                                                                                                ccGCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAgg                 >  1:820709/1‑151 (MQ=255)                                                                                                                                  gCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGc               >  2:620919/1‑151 (MQ=255)                                                                                                                                  gCTTCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGc               >  1:302142/1‑151 (MQ=255)                                                                                                                                    ttCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCtt             >  2:876454/1‑151 (MQ=255)                                                                                                                                    ttCGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAAT‑ATTGCCTGTAGCAGGGGAGAAAAAAACCAAAAGGAGAGCAGTGCCCAACCCCGGGCttt            >  2:86205/1‑151 (MQ=255)                                                                                                                                      cGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTc           >  2:118275/1‑151 (MQ=255)                                                                                                                                      cGGAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTc           <  2:323087/151‑1 (MQ=255)                                                                                                                                       ggAAAGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTcc          >  1:386338/1‑151 (MQ=255)                                                                                                                                         aaaGAGGAATGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAg        <  1:812249/151‑1 (MQ=255)                                                                                                                                                aaTGCCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTaaaa  >  2:327813/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TGTTGGTTTCGAACATCACAAACAGAACGCTCTGCGAGAACTCCTGACCGTAGATGACGTAGTAGCACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAA  >  NC_000913/4149941‑4150233

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