BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,896,607	T→C	T139A (ACG→GCG)	b3718 ←	putative isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,896,607	0	T	C	99.6%	1624.5 / ‑5.6	500	T139A (ACG→GCG)	b3718	putative isomerase
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (267/231); total (269/231)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 503 total aligned reads displayed.

AAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCAT  >  NC_000913/3896517‑3896696

aaaaGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCAt                                                                                         >  1:4807502/1‑93 (MQ=255)

aaaaGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCAt                                                                                         >  2:9067197/1‑93 (MQ=255)

  aaGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTc                                                                                       >  2:3878085/1‑93 (MQ=255)

   aGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCg                                                                                      <  2:4965839/93‑1 (MQ=255)

   aGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCg                                                                                      >  1:3776734/1‑93 (MQ=255)

     tCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGc                                                                                    <  1:6628345/93‑1 (MQ=255)

     tCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGc                                                                                    >  2:7852053/1‑93 (MQ=255)

     tCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGc                                                                                    <  2:1973505/93‑1 (MQ=255)

       gccgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGcca                                                                                  >  1:8259108/1‑91 (MQ=255)

       gccgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCaa                                                                                  <  1:3460459/93‑1 (MQ=255)

        ccgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCaaa                                                                                 <  2:6486965/93‑1 (MQ=255)

         cgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAAt                                                                                <  1:2693635/93‑1 (MQ=255)

         cgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAAt                                                                                <  1:9677339/93‑1 (MQ=255)

         cgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAAt                                                                                <  2:8150751/93‑1 (MQ=255)

          gccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAAtt                                                                               <  2:4442323/93‑1 (MQ=255)

          gccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAAtt                                                                               >  1:1500591/1‑93 (MQ=255)

           cccAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTg                                                                              >  2:1049914/1‑93 (MQ=255)

             cAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGcc                                                                            <  2:6048465/93‑1 (MQ=255)

              aaCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCCCGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCca                                                                           >  1:7082541/1‑93 (MQ=255)

              aaCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCca                                                                           <  2:4239679/93‑1 (MQ=255)

               aCTCGCCATGCGCCACAATATCATCCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCcac                                                                          <  2:5530853/93‑1 (MQ=255)

               aCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCcac                                                                          <  1:2225633/93‑1 (MQ=255)

               aCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCcac                                                                          <  1:4303798/93‑1 (MQ=255)

                cTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCcaca                                                                         <  1:3200933/93‑1 (MQ=255)

                  cGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGAAAATGGGTCGCATTCGGCAAATTGCCACaaa                                                                       <  1:5575110/93‑1 (MQ=255)

                   gCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACaaaa                                                                      >  1:1576008/1‑93 (MQ=255)

                     cATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAAt                                                                    >  1:9884506/1‑93 (MQ=255)

                     cATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAAt                                                                    <  1:8654145/93‑1 (MQ=255)

                      aTGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCCAATTGCCACAAAAATg                                                                   >  2:3350304/1‑93 (MQ=255)

                      aTGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATg                                                                   >  2:5878359/1‑93 (MQ=255)

                       tGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGa                                                                  >  1:8861725/1‑93 (MQ=255)

                        gcgcCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGAc                                                                 <  2:611680/93‑1 (MQ=255)

                        gcgcCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGAc                                                                 <  1:453872/93‑1 (MQ=255)

                         cgcCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGAcc                                                                <  2:2759807/93‑1 (MQ=255)

                         cgcCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGAcc                                                                <  2:8891469/93‑1 (MQ=255)

                          gcCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCa                                                               >  1:5046092/1‑93 (MQ=255)

                           ccACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCAt                                                              <  2:7196950/93‑1 (MQ=255)

                           ccACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCAt                                                              >  1:7483308/1‑93 (MQ=255)

                            cacaATATCAACCATCTCTCCGTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATc                                                             <  2:8641148/93‑1 (MQ=255)

                             acaATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCt                                                            >  2:2917799/1‑93 (MQ=255)

                             acaATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCt                                                            <  2:1563173/93‑1 (MQ=255)

                              caATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTg                                                           >  2:3565479/1‑93 (MQ=255)

                                 tatCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCAc                                                        >  2:9749946/1‑93 (MQ=255)

                                  atCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCAcc                                                       >  1:4098714/1‑93 (MQ=255)

                                    cAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTa                                                     <  2:8286059/93‑1 (MQ=255)

                                      aCCATCTCTCCCTGGATCGGGACCTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCAccctat                                                   <  2:5632048/93‑5 (MQ=255)

                                      aCCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAAt                                                   <  1:4430630/93‑1 (MQ=255)

                                          tctctcCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATTACCATCTGCACCTAATCCCa                                               <  1:9466951/93‑1 (MQ=255)

                                          tctctcCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCa                                               >  2:4145559/1‑93 (MQ=255)

                                           ctctcCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAg                                              <  2:6710371/93‑1 (MQ=255)

                                           ctctcCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAg                                              <  2:4820193/93‑1 (MQ=255)

                                           ctctcCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAg                                              <  1:1921388/93‑1 (MQ=255)

                                              tcCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcac                                           >  2:4639991/1‑93 (MQ=255)

                                              tcCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcac                                           <  1:912672/93‑1 (MQ=255)

                                               cccTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcacc                                          >  2:4387685/1‑93 (MQ=255)

                                                 cTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcaccac                                        >  1:2592480/1‑93 (MQ=255)

                                                 cTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcaccac                                        >  2:251898/1‑93 (MQ=255)

                                                  tGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCaccacc                                       >  1:613384/1‑93 (MQ=255)

                                                    gATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCaa                                     >  1:2757880/1‑93 (MQ=255)

                                                     aTCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCaaa                                    <  2:8070536/93‑1 (MQ=255)

                                                     aTCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCaaa                                    >  2:135013/1‑93 (MQ=255)

                                                       cGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATc                                  <  1:7998043/93‑1 (MQ=255)

                                                        gggAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATcc                                 <  1:1679929/93‑1 (MQ=255)

                                                        gggAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATcc                                 >  1:9733902/1‑93 (MQ=255)

                                                         ggAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCa                                >  1:1230281/1‑93 (MQ=255)

                                                         ggAACTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCa                                >  2:2217718/1‑93 (MQ=255)

                                                           aaCTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGt                              <  1:4228658/93‑1 (MQ=255)

                                                            aCTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTc                             <  2:7456288/93‑1 (MQ=255)

                                                             cTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTcc                            >  2:9529087/1‑93 (MQ=255)

                                                             cTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTcc                            <  2:1088413/93‑1 (MQ=255)

                                                             cTCCACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTcc                            >  2:5764005/1‑93 (MQ=255)

                                                               ccACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGc                          >  1:6419881/1‑93 (MQ=255)

                                                                cACGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGcc                         <  1:2393434/93‑1 (MQ=255)

                                                                 aCGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCt                        >  2:2340587/1‑93 (MQ=255)

                                                                  cGGTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCtt                       >  2:7665261/1‑93 (MQ=255)

                                                                    gTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCa                     <  2:3782506/93‑1 (MQ=255)

                                                                    gTTTGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCa                     >  1:5956284/1‑93 (MQ=255)

                                                                     tttGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAc                    <  1:6938001/93‑1 (MQ=255)

                                                                      ttGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcg                   >  1:4696810/1‑93 (MQ=255)

                                                                       tGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcgc                  >  2:9326303/1‑93 (MQ=255)

                                                                       tGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcgc                  >  1:9264103/1‑93 (MQ=255)

                                                                       tGCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcgc                  >  2:1703183/1‑93 (MQ=255)

                                                                        gCTCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcgcg                 >  1:7293466/1‑93 (MQ=255)

                                                                          tCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGcc               >  1:7683456/1‑93 (MQ=255)

                                                                          tCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGcc               <  1:5671415/93‑1 (MQ=255)

                                                                          tCATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGcc               >  1:5420654/1‑93 (MQ=255)

                                                                           cATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCACATCCAGTCCGCCTTCACGCGCCa              >  2:6859514/1‑93 (MQ=255)

                                                                           cATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCa              <  1:8435235/93‑1 (MQ=255)

                                                                           cATGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCa              >  2:5170302/1‑93 (MQ=255)

                                                                            aTGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAg             >  2:3151051/1‑93 (MQ=255)

                                                                            aTGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAg             <  2:2578097/93‑1 (MQ=255)

                                                                             tGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGt            >  2:6263132/1‑93 (MQ=255)

                                                                             tGGAAATGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGt            >  1:211184/1‑93 (MQ=255)

                                                                                 aaTGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCt        <  2:8461914/93‑1 (MQ=255)

                                                                                 aaTGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCt        >  1:742064/1‑93 (MQ=255)

                                                                                  aTGGGTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTg       >  1:6845669/1‑93 (MQ=255)

                                                                                      gTCGCATTCGGCAAATTGCCACATAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCa   <  2:6931899/93‑1 (MQ=255)

                                                                                      gTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCa   <  1:10011690/93‑1 (MQ=255)

                                                                                      gTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCa   <  1:3079487/93‑1 (MQ=255)

                                                                                      gTCGCATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCa   <  2:3680643/93‑1 (MQ=255)

AAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCAT  >  NC_000913/3896517‑3896696

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