Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,114,502	C→T	*465* (TAG→TAA)	wcaM ←	colanic acid biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,114,502	0	C	T	100.0%	250.7 / NA	77	*465* (TAG→TAA)	wcaM	colanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/42);  total (35/42)

TGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGT  >  NC_000913/2114367‑2114631                                                                                                                                        |                                                                                                                                  tGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTAc                                                                                                                                >  1:446705/1‑139 (MQ=255)   ttCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTc                                                                                                                              >  1:276121/1‑139 (MQ=255)     cGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                                            <  2:485674/139‑1 (MQ=255)      gATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                                           <  1:824739/139‑1 (MQ=255)         ttttttCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAg                                                                                                                        <  1:327444/139‑1 (MQ=255)               cGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                                  <  2:918154/139‑1 (MQ=255)                    gCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATc                                                                                                             <  2:627941/139‑1 (MQ=255)                       tCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTaa                                                                                                          <  1:923155/139‑1 (MQ=255)                          aTTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGt                                                                                                       <  2:939679/139‑1 (MQ=255)                            tACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                                     <  2:528568/139‑1 (MQ=255)                              cccTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAt                                                                                                   <  2:141514/139‑1 (MQ=255)                               ccTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtc                                                                                                  >  2:578557/1‑139 (MQ=255)                                    aTTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             <  2:232794/139‑1 (MQ=255)                                     ttGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                                  >  2:1007311/1‑117 (MQ=255)                                     ttGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                                  <  1:1007311/117‑1 (MQ=255)                                     ttGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTcc                                                                                            <  2:100814/139‑1 (MQ=255)                                     ttGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTcc                                                                                            <  2:276121/139‑1 (MQ=255)                                       gATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCcac                                                                                          <  2:491923/139‑1 (MQ=255)                                        aTTACTGCATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTaa                                                                                                    <  2:818953/128‑1 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATc                                                                                                             <  2:896464/119‑1 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATc                                                                                                             >  1:896464/1‑119 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTaa                                                                                                    >  1:818953/1‑128 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCcaca                                                                                         >  2:1105378/1‑139 (MQ=255)                                         ttACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             <  2:699898/134‑1 (MQ=255)                                         ttACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             >  1:699898/1‑134 (MQ=255)                                         ttACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        <  1:472845/139‑1 (MQ=255)                                          tACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAtt                                                                                       <  1:1105378/139‑1 (MQ=255)                                               aaTTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGt                                                                                                       >  2:767130/1‑118 (MQ=255)                                               aaTTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGt                                                                                                       <  1:767130/118‑1 (MQ=255)                                                   aCTTGTGATGTTACGCCGCCTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                              <  2:599291/139‑1 (MQ=255)                                                       gtgATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             >  2:45054/1‑120 (MQ=255)                                                       gtgATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             <  1:45054/120‑1 (MQ=255)                                                       gtgATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGa                                                                          >  2:119215/1‑139 (MQ=255)                                                         gATGTTATGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGAtt                                                                        >  2:500624/1‑139 (MQ=255)                                                            gTTATGCTGCTTTCTTGTTGATTGCAGTAGCATTGTTCCTAAGTACGCCTCCCTCTTTCCAGGAATTGTCGCAAATTTTCTCCCTCAGTTCCGGCATTCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAAt                                                                     <  1:500624/139‑1 (MQ=255)                                                                  ccgcTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTc                                                               >  1:964613/1‑139 (MQ=255)                                                                     cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGAt                                                                         >  1:343071/1‑126 (MQ=255)                                                                     cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGAt                                                                         <  2:343071/126‑1 (MQ=255)                                                                     cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGat                                                            >  1:200179/1‑139 (MQ=255)                                                                                  tGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                     <  2:110624/101‑1 (MQ=255)                                                                                  tGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                                     >  1:110624/1‑101 (MQ=255)                                                                                       tAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTcc                                          >  2:812773/1‑139 (MQ=255)                                                                                       tAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTcc                                          >  1:709549/1‑139 (MQ=255)                                                                                        aGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCt                                         <  2:709549/139‑1 (MQ=255)                                                                                          cATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCttt                                       >  1:77571/1‑139 (MQ=255)                                                                                                   tAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCAc                              <  1:812773/139‑1 (MQ=255)                                                                                                   tAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCAc                              <  2:200179/139‑1 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttc                                                                                               >  1:1180080/1‑68 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttc                                                                                               <  2:1180080/68‑1 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             >  1:1030444/1‑70 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                                             <  2:1030444/70‑1 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGat                                                            <  2:654946/103‑1 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGat                                                            >  1:654946/1‑103 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACg                                                     <  1:1096271/110‑1 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACg                                                     >  2:1096271/1‑110 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATg                        >  1:560579/1‑139 (MQ=255)                                                                                                         tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATg                        >  2:634265/1‑139 (MQ=255)                                                                                                           aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATc                                                                                                             >  2:416801/1‑52 (MQ=255)                                                                                                           aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATc                                                                                                             <  1:416801/52‑1 (MQ=255)                                                                                                           aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTaa                                                                                                          >  1:599639/1‑55 (MQ=255)                                                                                                           aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTaa                                                                                                          <  2:599639/55‑1 (MQ=255)                                                                                                           aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAata                                                                                <  2:1102796/81‑1 (MQ=255)                                                                                                           aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAata                                                                                >  1:1102796/1‑81 (MQ=255)                                                                                                            cTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                                            <  2:334420/36‑1 (MQ=255)                                                                                                            cTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                                            >  1:334420/1‑36 (MQ=255)                                                                                                            cTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTa                                                                       >  1:300358/1‑89 (MQ=255)                                                                                                            cTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTa                                                                       <  2:300358/89‑1 (MQ=255)                                                                                                              ccATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACAt                   <  1:578557/139‑1 (MQ=255)                                                                                                               cATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTcc                                                                                            >  1:996972/1‑65 (MQ=255)                                                                                                               cATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTcc                                                                                            <  2:996972/65‑1 (MQ=255)                                                                                                                           aaTGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTc                                           <  1:584342/102‑1 (MQ=255)                                                                                                                           aaTGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTc                                           >  2:584342/1‑102 (MQ=255)                                                                                                                             tGGTCGCAAATTTACTCTCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGAcagca    <  2:446705/139‑1 (MQ=255)                                                                                                                             tGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGAcagca    >  1:953286/1‑139 (MQ=255)                                                                                                                               gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        >  1:493514/1‑53 (MQ=255)                                                                                                                               gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                                        <  2:493514/53‑1 (MQ=255)                                                                                                                               gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAgt  >  2:1074875/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                  TGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGT  >  NC_000913/2114367‑2114631

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