BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,116,015	A→G	F138L (TTT→CTT)	wcaK ←	putative galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,116,015	0	A	G	99.6%	695.9 / ‑5.3	221	F138L (TTT→CTT)	wcaK	putative galactokinase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (120/100); total (121/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 226 total aligned reads displayed.

GGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGA  >  NC_000913/2115925‑2116104

ggTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCt                                                                                         >  1:3397663/1‑93 (MQ=255)

 gTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTg                                                                                        >  2:2899306/1‑93 (MQ=255)

  ttAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGc                                                                                       >  2:563089/1‑93 (MQ=255)

  atAAATTGCTCAACATAAAACGGGCCAACACTGTAACCAATCATAAACAGCAGCTTATACGCCATAAACGTGCAAAGTGCATGTTAAAGCTGc                                                                                       <  1:7038886/92‑1 (MQ=255)

     aaTTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGc                                                                                    >  2:3245073/1‑93 (MQ=255)

      aTTGCTCATCCTGAAACGGGCCAACACTGTGACTAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCa                                                                                   <  1:7630148/93‑1 (MQ=255)

       ttGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCACAGTGCATGTTCAAGCTGCGGCAc                                                                                  >  2:473468/1‑93 (MQ=255)

       ttGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCAc                                                                                  <  2:3099910/93‑1 (MQ=255)

        tGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACg                                                                                 >  1:703678/1‑93 (MQ=255)

          cTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGcc                                                                               >  1:6735297/1‑93 (MQ=255)

           tCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGAGGAACGCCa                                                                              >  2:822404/1‑93 (MQ=255)

            cATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCAt                                                                             >  2:4986151/1‑93 (MQ=255)

            cATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCAt                                                                             <  2:2269239/93‑1 (MQ=255)

               ccTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATaga                                                                          <  2:895640/93‑1 (MQ=255)

               ccTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATaga                                                                          <  2:5934413/93‑1 (MQ=255)

                cTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATagag                                                                         >  2:1213204/1‑93 (MQ=255)

                cTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATagag                                                                         >  1:2557261/1‑93 (MQ=255)

                 tGGAACGTGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAgaga                                                                        <  1:58201/93‑1 (MQ=255)

                 tGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAgaga                                                                        >  2:5532769/1‑93 (MQ=255)

                 tGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAgaga                                                                        <  1:2248748/93‑1 (MQ=255)

                 tGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAgaga                                                                        <  2:2033661/93‑1 (MQ=255)

                   gAACGGTCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATc                                                                      <  2:559226/93‑1 (MQ=255)

                   gAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATc                                                                      >  2:3774321/1‑93 (MQ=255)

                   gAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATc                                                                      <  2:3675565/93‑1 (MQ=255)

                    aaCGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCg                                                                     <  1:2821533/93‑1 (MQ=255)

                       gggCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACa                                                                  >  1:4637480/1‑93 (MQ=255)

                          gcAACACTGTGACGAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACaaaa                                                               <  1:4809285/92‑1 (MQ=255)

                            aaCACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACaaaaaa                                                             <  1:6488082/93‑1 (MQ=255)

                            aaCACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACaaaaaa                                                             >  1:4759766/1‑93 (MQ=255)

                            aaCACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACaaaaaa                                                             <  1:7404204/93‑1 (MQ=255)

                              cacTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACg                                                           <  2:6360646/93‑1 (MQ=255)

                               acTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGa                                                          <  2:6929244/93‑1 (MQ=255)

                               acTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGa                                                          >  2:60170/1‑93 (MQ=255)

                                cTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGAt                                                         <  1:7104944/93‑1 (MQ=255)

                                cTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGAt                                                         <  1:365144/93‑1 (MQ=255)

                                cTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGAt                                                         >  2:6872679/1‑93 (MQ=255)

                                 tgtgACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATc                                                        >  1:5985632/1‑93 (MQ=255)

                                 tgtgACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATc                                                        >  2:6389830/1‑93 (MQ=255)

                                  gtgACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATcc                                                       <  2:2952372/93‑1 (MQ=255)

                                  gtgACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATcc                                                       <  1:2971547/93‑1 (MQ=255)

                                  gtgACCAATCATAAACAGCGGCTTTTGCGCCATAAACGTGCAAACTTCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATcc                                                       <  1:7442188/93‑1 (MQ=255)

                                   tgACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATccg                                                      >  1:2597890/1‑93 (MQ=255)

                                    gACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATccgc                                                     >  2:3842226/1‑93 (MQ=255)

                                     aCCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATccgcc                                                    <  2:4266837/93‑1 (MQ=255)

                                       cAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATccgccgc                                                  >  1:3241428/1‑92 (MQ=255)

                                        aaTCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGAc                                                 >  1:7022298/1‑93 (MQ=255)

                                        aaTCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGAc                                                 >  1:5325199/1‑93 (MQ=255)

                                         aTCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGAcc                                                >  2:7078130/1‑93 (MQ=255)

                                           cATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTg                                              <  2:7521407/93‑1 (MQ=255)

                                           cATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTg                                              <  2:733766/93‑1 (MQ=255)

                                           cATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTg                                              <  2:1445985/93‑1 (MQ=255)

                                             tAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGa                                            >  2:2519808/1‑93 (MQ=255)

                                              aaaCAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGAt                                           >  2:5403345/1‑93 (MQ=255)

                                                aCAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATaa                                         <  1:3194395/93‑1 (MQ=255)

                                                aCAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATaa                                         >  2:5276990/1‑93 (MQ=255)

                                                aCAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATaa                                         <  1:5436469/93‑1 (MQ=255)

                                                 cAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAAt                                        <  1:2358828/93‑1 (MQ=255)

                                                 cAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAAt                                        >  1:3761495/1‑93 (MQ=255)

                                                 cAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAAt                                        >  1:7515472/1‑93 (MQ=255)

                                                  aGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATg                                       <  1:4037571/93‑1 (MQ=255)

                                                   gCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATgg                                      <  2:481362/93‑1 (MQ=255)

                                                   gCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATgg                                      >  1:1163833/1‑93 (MQ=255)

                                                    cGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGc                                     <  1:1459712/93‑1 (MQ=255)

                                                      gCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGcgt                                   <  1:6796603/93‑1 (MQ=255)

                                                      gCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGcgt                                   >  2:6156890/1‑93 (MQ=255)

                                                      gCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGcgt                                   >  2:4424157/1‑93 (MQ=255)

                                                       cTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGcgtc                                  >  1:6359904/1‑93 (MQ=255)

                                                         ttttCGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGcgtcgt                                >  1:1655336/1‑93 (MQ=255)

                                                             cGCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAAcc                            <  2:1610638/93‑1 (MQ=255)

                                                              gCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCt                           <  2:3627001/93‑1 (MQ=255)

                                                              gCCATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCt                           >  1:7209721/1‑93 (MQ=255)

                                                                cATAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGa                         <  1:5789109/93‑1 (MQ=255)

                                                                 aTAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAc                        <  2:1958520/93‑1 (MQ=255)

                                                                 aTAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAc                        >  2:6225191/1‑93 (MQ=255)

                                                                  tAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAca                       >  1:2335235/1‑93 (MQ=255)

                                                                  tAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAca                       >  2:4078620/1‑93 (MQ=255)

                                                                  tAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAca                       >  2:6350239/1‑93 (MQ=255)

                                                                  tAAACGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAca                       <  2:1697265/93‑1 (MQ=255)

                                                                     aCGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAcagca                    <  2:7736834/93‑1 (MQ=255)

                                                                     aCGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAcagca                    <  2:6632783/93‑1 (MQ=255)

                                                                     aCGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAcagca                    <  1:3685005/93‑1 (MQ=255)

                                                                     aCGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAcagca                    <  1:4917404/93‑1 (MQ=255)

                                                                     aCGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAcagca                    <  2:5191433/93‑1 (MQ=255)

                                                                      cGTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGAcagcag                   <  2:5804645/93‑1 (MQ=255)

                                                                       gTGCAAAGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAgg                  >  1:6318209/1‑93 (MQ=255)

                                                                           aaaGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCa              <  1:2052612/93‑1 (MQ=255)

                                                                           aaaGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCa              >  1:5366528/1‑93 (MQ=255)

                                                                           aaaGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCa              >  1:728454/1‑93 (MQ=255)

                                                                             aGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACg            >  1:7806791/1‑93 (MQ=255)

                                                                             aGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACg            <  1:5106639/93‑1 (MQ=255)

                                                                             aGTGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACg            >  2:876725/1‑93 (MQ=255)

                                                                               tGCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGaa          <  2:7484979/93‑1 (MQ=255)

                                                                               ggCATGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGaa          <  1:4923885/92‑1 (MQ=255)

                                                                                  aTGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTc       <  2:5323803/93‑1 (MQ=255)

                                                                                   tGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg      >  2:4005397/1‑93 (MQ=255)

                                                                                   tGTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCg      >  1:4488704/1‑93 (MQ=255)

                                                                                    gTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCgg     >  1:3347376/1‑93 (MQ=255)

                                                                                    gTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCgg     >  2:6279656/1‑93 (MQ=255)

                                                                                    gTTCAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCgg     >  1:4155951/1‑93 (MQ=255)

                                                                                       cAAGCTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGa  >  1:1405189/1‑93 (MQ=255)

GGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGA  >  NC_000913/2115925‑2116104

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