Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,273,906	T→C	intergenic (‑35/‑219)	uvrA ← / → ssb	ATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC/single‑stranded DNA‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,273,906	0	T	C	100.0%	207.7 / NA	59	intergenic (‑35/‑219)	uvrA/ssb	ATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC/single‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/33);  total (26/33)

GGATTTGCCAGAACCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGTGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAA  >  NC_000913/4273758‑4274053                                                                                                                                                     |                                                                                                                                                    ggATTTGCCAGAACCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCcggcgg                                                                                                                                                   <  1:198966/151‑1 (MQ=255) ggATTTGCCAGAACCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCcggcgg                                                                                                                                                   <  2:436181/151‑1 (MQ=255)  gATTTGCCAGAACCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCggcggc                                                                                                                                                  >  2:87620/1‑151 (MQ=255)  gATTTGCCAGAACCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCggcggc                                                                                                                                                  >  2:61878/1‑151 (MQ=255)         cAGAACCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGc                                                                                                                                           >  1:504495/1‑151 (MQ=255)          aGAACCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCt                                                                                                                                          <  1:366986/151‑1 (MQ=255)          aGAACCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCt                                                                                                                                          <  1:394579/151‑1 (MQ=255)             aCCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTAcc                                                                                                                                       >  1:732119/1‑151 (MQ=255)                cGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGc                                                                                                                                    >  1:258443/1‑151 (MQ=255)                    aGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGAt                                                                                                                                 <  2:698998/150‑1 (MQ=255)                    aGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATc                                                                                                                                <  2:633902/151‑1 (MQ=255)                     gCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATca                                                                                                                               <  2:566557/151‑1 (MQ=255)                      cccGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATcac                                                                                                                              >  1:531556/1‑151 (MQ=255)                        cGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACaa                                                                                                                            <  2:336172/151‑1 (MQ=255)                          gTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAAc                                                                                                                          <  1:259769/151‑1 (MQ=255)                                aCAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTa                                                                                                                     >  2:467139/1‑150 (MQ=255)                                 cAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAAt                                                                                                                   >  2:543697/1‑151 (MQ=255)                                  aaTGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAAtt                                                                                                                  >  1:236840/1‑151 (MQ=255)                                   aTGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAAtt                                                                                                                  >  2:549330/1‑150 (MQ=255)                                   aTGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTa                                                                                                                 <  2:738603/151‑1 (MQ=255)                                     gAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTAt                                                                                                                <  2:376902/150‑1 (MQ=255)                                      aGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATg                                                                                                               >  2:382604/1‑150 (MQ=255)                                          tGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGacac                                                                                                           >  2:820805/1‑150 (MQ=255)                                           gTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGAcaca                                                                                                          <  2:427837/150‑1 (MQ=255)                                           gTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACaa                                                                                                         <  2:494746/151‑1 (MQ=255)                                              gcgGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACaaa                                                                                                        <  1:410126/149‑1 (MQ=255)                                                gggggATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTg                                                                                                     >  2:415576/1‑150 (MQ=255)                                                   ggATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCt                                                                                                 <  1:68285/151‑1 (MQ=255)                                                       aaCGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCtgaa                                                                                              <  1:475150/150‑1 (MQ=255)                                                             gTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAatat                                                                                       >  2:567763/1‑151 (MQ=255)                                                                 aTGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAg                                                                                   <  2:53102/151‑1 (MQ=255)                                                                    tttttGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTAt                                                                                >  2:269552/1‑151 (MQ=255)                                                                     ttttGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTAt                                                                                <  2:133409/150‑1 (MQ=255)                                                                       ttGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTg                                                                              <  1:582066/150‑1 (MQ=255)                                                                           gaTTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATg                                                                         >  2:134520/1‑151 (MQ=255)                                                                            aTTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGc                                                                        >  1:322993/1‑151 (MQ=255)                                                                            aTTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGc                                                                        <  2:236864/151‑1 (MQ=255)                                                                               aTGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCAt                                                                      <  1:742121/150‑1 (MQ=255)                                                                                 gggTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTa                                                                    <  2:761034/150‑1 (MQ=255)                                                                                   gTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTAccc                                                                 <  2:66367/151‑1 (MQ=255)                                                                                   gTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTAccc                                                                 <  1:61855/151‑1 (MQ=255)                                                                                     gCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCGCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCgg                                                               <  2:817733/151‑1 (MQ=255)                                                                                      cGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGa                                                              <  1:254776/151‑1 (MQ=255)                                                                                           gcCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGtgtt                                                         >  2:315504/1‑151 (MQ=255)                                                                                             cccGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGtgttgt                                                       <  2:487822/151‑1 (MQ=255)                                                                                                  aCTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAAc                                                  >  2:41844/1‑151 (MQ=255)                                                                                                    ttCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACa                                                 >  2:69511/1‑150 (MQ=255)                                                                                                                 aTTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGCTGCTACCGGCGATCACAAACGGATAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAg                                     >  1:426726/1‑149 (MQ=255)                                                                                                                   tCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGgttt                                 >  1:440848/1‑151 (MQ=255)                                                                                                                         ttCCCGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTc                            <  1:344321/150‑1 (MQ=255)                                                                                                                           cccGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCg                         >  1:131509/1‑151 (MQ=255)                                                                                                                           cccGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCg                         <  1:605825/151‑1 (MQ=255)                                                                                                                            ccGGATTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCgg                        >  1:431069/1‑151 (MQ=255)                                                                                                                                aTTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAAcc                    <  1:48593/151‑1 (MQ=255)                                                                                                                                aTTAAACGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAAcc                    <  1:315480/151‑1 (MQ=255)                                                                                                                                    aaCGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAgg                >  1:781528/1‑151 (MQ=255)                                                                                                                                     aCGCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGt               <  1:563799/151‑1 (MQ=255)                                                                                                                                       gCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTca             <  1:855593/151‑1 (MQ=255)                                                                                                                                       gCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTca             <  1:737831/151‑1 (MQ=255)                                                                                                                                       gCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTca             <  1:526801/151‑1 (MQ=255)                                                                                                                                       gCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTc              >  1:786823/1‑150 (MQ=255)                                                                                                                                       gCTTTTTTGCCCGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTc              >  1:319433/1‑150 (MQ=255)                                                                                                                                                 ccGGCGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGt    >  2:333564/1‑150 (MQ=255)                                                                                                                                                  cggcggCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATAAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTaa  <  2:23458/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGATTTGCCAGAACCCGAAAGCCCGGTCACGACAATGAGCTTGTCGCGGGGGATAACGAGGTTGATGTTTTTGAGATTATGGGTGCGGGCGCCCCGAACTTCGATCTTATCCATTCACCTTTCCCGGATTAAACGCTTTTTTGCCCGGTGGCATGGTGCTACCGGCGATCACAAACGGTTAATTATGACACAAATTGACCTGAATGAATATACAGTATTGGAATGCATTACCCGGAGTGTTGTGTAACAATGTCTGGCCAGGTTTGTTTCCCGGAACCGAGGTCACAACATAGTAA  >  NC_000913/4273758‑4274053

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