BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,895,138	T→C	S2P (TCT→CCT)	yieI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,895,138	0	T	C	99.7%	934.9 / ‑5.1	293	S2P (TCT→CCT)	yieI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (129/163); total (129/164)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 297 total aligned reads displayed.

GGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTC  >  NC_000913/3895048‑3895227

ggTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCt                                                                                         >  1:3992739/1‑93 (MQ=255)

ggTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCt                                                                                         <  1:1142184/93‑1 (MQ=255)

ggTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCt                                                                                         <  2:2597457/93‑1 (MQ=255)

ggTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTCAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCt                                                                                         >  2:3941364/1‑93 (MQ=255)

 gTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTg                                                                                        >  1:1448738/1‑93 (MQ=255)

 gTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTg                                                                                        <  2:956724/93‑1 (MQ=255)

   ttggATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGtt                                                                                      <  1:1036993/91‑1 (MQ=255)

   tGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGtt                                                                                      <  2:3012101/93‑1 (MQ=255)

   tGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGtt                                                                                      <  1:925162/93‑1 (MQ=255)

    gggATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGttt                                                                                     >  1:2642585/1‑93 (MQ=255)

     ggATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTc                                                                                    >  2:2056954/1‑93 (MQ=255)

     ggATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTc                                                                                    <  2:4524881/93‑1 (MQ=255)

      gATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCa                                                                                   <  2:4250482/93‑1 (MQ=255)

      gATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCa                                                                                   <  2:814998/93‑1 (MQ=255)

       atatTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCAc                                                                                  <  2:3338394/93‑1 (MQ=255)

          ttACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTc                                                                               >  1:3662499/1‑93 (MQ=255)

            aCGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCgg                                                                             >  2:1221204/1‑93 (MQ=255)

             cGGCATAGTTCTTCACACTCCCTTCACTTCCCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCggg                                                                            <  2:2240520/93‑1 (MQ=255)

               gCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTa                                                                          >  2:2180486/1‑93 (MQ=255)

               gCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTa                                                                          >  1:2784117/1‑93 (MQ=255)

                cATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTaa                                                                         >  1:4067157/1‑93 (MQ=255)

                cATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTaa                                                                         >  2:3325864/1‑93 (MQ=255)

                 aTAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAAt                                                                        >  1:1324848/1‑93 (MQ=255)

                  tAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATa                                                                       >  1:2425004/1‑93 (MQ=255)

                  tAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATa                                                                       >  1:2727896/1‑93 (MQ=255)

                  tAGTTCTGCGCGCTCCCTTCACTTCCCCCGCTTAGATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATa                                                                       <  2:3690578/93‑1 (MQ=255)

                    gTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACa                                                                     <  1:3720318/93‑1 (MQ=255)

                      tcttcCCCCTCCCTTCACTTCCCCCGCTTAAATGGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATc                                                                   <  1:1636942/93‑1 (MQ=255)

                      tcttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATc                                                                   <  1:2061480/93‑1 (MQ=255)

                       cttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCa                                                                  <  2:3464725/93‑1 (MQ=255)

                        ttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCAc                                                                 >  1:3076449/1‑93 (MQ=255)

                        ttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCAc                                                                 <  2:2938143/93‑1 (MQ=255)

                           acacTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGa                                                              >  2:3262729/1‑93 (MQ=255)

                            cacTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGAc                                                             <  1:1096558/93‑1 (MQ=255)

                             ccTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACt                                                            <  1:4260473/92‑1 (MQ=255)

                             acTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACt                                                            >  2:1090638/1‑93 (MQ=255)

                             acTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACt                                                            >  2:3775249/1‑93 (MQ=255)

                             acTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGGAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACt                                                            >  2:3549019/1‑93 (MQ=255)

                               tCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACttt                                                          >  1:1576496/1‑93 (MQ=255)

                                 ccTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAt                                                        <  2:778327/93‑1 (MQ=255)

                                  cTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                                                       <  1:3417821/93‑1 (MQ=255)

                                  cTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                                                       <  2:1899967/93‑1 (MQ=255)

                                  cTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                                                       <  1:4456384/93‑1 (MQ=255)

                                   ttCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAttt                                                      <  1:2250085/93‑1 (MQ=255)

                                   ttCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAttt                                                      <  1:294492/93‑1 (MQ=255)

                                    tCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAtttt                                                     >  1:3971547/1‑93 (MQ=255)

                                    tCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAtttt                                                     <  2:1117782/93‑1 (MQ=255)

                                     cACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGCGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTg                                                    >  1:2514850/1‑93 (MQ=255)

                                       cTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCa                                                  <  1:3278534/93‑1 (MQ=255)

                                       cTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCa                                                  >  2:1445529/1‑93 (MQ=255)

                                       cTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCa                                                  <  2:1826476/93‑1 (MQ=255)

                                        ttACCCCGCTTAAATTGGCGCTCAAAGGGAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACCTCAAGGATTTTATTTTGCAg                                                 >  2:3102091/1‑93 (MQ=255)

                                         tACCCCGCTTTAATTGGCGCTCAAAGGTAAATAAAGGTAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGt                                                <  2:3744206/93‑1 (MQ=255)

                                         tACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGt                                                <  2:2155027/93‑1 (MQ=255)

                                          aCCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGtt                                               <  2:1667491/93‑1 (MQ=255)

                                           ccccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGttt                                              <  2:840900/93‑1 (MQ=255)

                                              cGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAg                                           >  1:352066/1‑93 (MQ=255)

                                              cGCTTAAATTCGCACTCAAAGGTAAGTAAAGCGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAg                                           <  1:4190409/93‑1 (MQ=255)

                                               gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAgg                                          <  2:238119/93‑1 (MQ=255)

                                               gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAgg                                          >  2:972937/1‑93 (MQ=255)

                                               gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAgg                                          >  2:2535331/1‑93 (MQ=255)

                                               gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAgg                                          >  2:3396850/1‑93 (MQ=255)

                                               gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAgg                                          >  1:801914/1‑93 (MQ=255)

                                                cTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGa                                         <  1:3546938/93‑1 (MQ=255)

                                                 ttAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGAc                                        >  1:1255024/1‑93 (MQ=255)

                                                 ttAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGAc                                        <  1:461144/93‑1 (MQ=255)

                                                  tAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGAcc                                       <  1:3861690/93‑1 (MQ=255)

                                                  tAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGAcc                                       >  1:2286187/1‑93 (MQ=255)

                                                    aaTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGt                                     >  2:4418307/1‑93 (MQ=255)

                                                    aaTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGt                                     >  1:2123360/1‑93 (MQ=255)

                                                     aTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGtt                                    >  1:3449571/1‑93 (MQ=255)

                                                       tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTaa                                  >  2:1606111/1‑93 (MQ=255)

                                                          cgcTCAAAGGTAAGTAAAGAGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTg                               <  1:681583/93‑1 (MQ=255)

                                                              cAAAGGTAAATAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGgtgt                           <  1:2142433/93‑1 (MQ=255)

                                                               aaaGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGtt                          >  1:619029/1‑93 (MQ=255)

                                                                aaGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTc                         <  2:2453448/93‑1 (MQ=255)

                                                                aaGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTc                         >  2:2826739/1‑93 (MQ=255)

                                                                  ggTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTg                       <  1:18372/93‑1 (MQ=255)

                                                                  ggTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTg                       >  1:3778851/1‑93 (MQ=255)

                                                                  ggTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGATTTAGGACCGTTAATTGGTGTTCTg                       >  2:1275770/1‑93 (MQ=255)

                                                                   gtaagtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGt                      >  2:1561939/1‑93 (MQ=255)

                                                                   gtaagtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGt                      <  1:1967772/93‑1 (MQ=255)

                                                                   gtaagtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGt                      >  1:4087660/1‑93 (MQ=255)

                                                                    taagtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGtt                     <  1:1485100/93‑1 (MQ=255)

                                                                      agtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTc                   <  1:3146135/93‑1 (MQ=255)

                                                                       gtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCt                  >  1:842704/1‑93 (MQ=255)

                                                                        taaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCtt                 <  1:4591682/93‑1 (MQ=255)

                                                                        taaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCtt                 <  2:4586498/93‑1 (MQ=255)

                                                                         aaaGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTg                >  1:499976/1‑93 (MQ=255)

                                                                         aaaGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTg                >  2:4294238/1‑93 (MQ=255)

                                                                          aaGGTAGTTTGATATGCTTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGt               <  2:456021/93‑1 (MQ=255)

                                                                           aGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTc              <  1:4537391/93‑1 (MQ=255)

                                                                            gggAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTcc             <  2:179596/93‑1 (MQ=255)

                                                                               aGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCctct          >  1:1940558/1‑93 (MQ=255)

                                                                                gTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTa         >  1:4107692/1‑93 (MQ=255)

                                                                                   tGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACat      >  1:3633640/1‑93 (MQ=255)

                                                                                      tatGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATAtt   <  1:1208531/93‑1 (MQ=255)

                                                                                      tatGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATAtt   >  2:4058721/1‑93 (MQ=255)

                                                                                      tatGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATAtt   >  1:4650326/1‑93 (MQ=255)

                                                                                       atGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATAttc  >  2:2310059/1‑93 (MQ=255)

GGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTC  >  NC_000913/3895048‑3895227

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: `‑`