BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,388,932	A→G	A52A (GCT→GCC)	b4161 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,388,932	0	A	G	97.9%	438.2 / ‑0.5	143	A52A (GCT→GCC)	b4161	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (66/74); total (67/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

Only 100 of 147 total aligned reads displayed.

ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTT  >  NC_000913/4388842‑4389020

aCGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGc                                                                                        <  1:1364687/93‑1 (MQ=255)

aCGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGc                                                                                        <  1:4249025/93‑1 (MQ=255)

  gCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGt                                                                                      <  2:493053/93‑1 (MQ=255)

    cAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGc                                                                                    <  1:3844738/93‑1 (MQ=255)

    cAGACTACGCGGTCGCCGGTTACCAGCGAACGGACGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGc                                                                                    >  1:4422826/1‑93 (MQ=255)

     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                   <  2:2186078/93‑1 (MQ=255)

     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                   <  1:94374/93‑1 (MQ=255)

     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                   <  1:3526600/93‑1 (MQ=255)

     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                   <  1:1622582/93‑1 (MQ=255)

     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                   >  2:1361986/1‑93 (MQ=255)

       aCTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATa                                                                                 <  2:736575/93‑1 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                               >  1:1292405/1‑93 (MQ=255)

          aCGCGGTCGCCGGTTACGAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCTCATCGGCGTGCATGCCa                                                                              <  2:212368/93‑1 (MQ=255)

             cgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAg                                                                           <  1:1137043/93‑1 (MQ=255)

              ggTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGc                                                                          >  1:3648536/1‑93 (MQ=255)

                 tgcaggTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGc                                                                       <  2:4236622/89‑1 (MQ=255)

                  gCCGGTTACCAGCGAACGGATGGTACGGCAAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCt                                                                      <  2:1762055/93‑1 (MQ=255)

                   ccGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtg                                                                     >  2:1000356/1‑93 (MQ=255)

                    cGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtga                                                                    <  1:3613986/93‑1 (MQ=255)

                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                  <  1:2479576/93‑1 (MQ=255)

                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                  <  1:2605321/93‑1 (MQ=255)

                       gtACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatga                                                                 <  1:607602/92‑1 (MQ=255)

                        tACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGAc                                                                <  2:3361304/93‑1 (MQ=255)

                           cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGAt                                                             <  2:3098715/93‑1 (MQ=255)

                            aGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATa                                                            <  1:2048097/93‑1 (MQ=255)

                            aGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATCCCAAAGCGGCTGATGCCGATa                                                            <  2:2575741/93‑1 (MQ=255)

                              cGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAcc                                                          <  2:1771472/93‑1 (MQ=255)

                               gAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCt                                                         <  1:1713983/93‑1 (MQ=255)

                               gAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCt                                                         <  2:3673967/93‑1 (MQ=255)

                                aaCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCtt                                                        >  2:314989/1‑93 (MQ=255)

                                 aCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtc                                                       >  2:492150/1‑93 (MQ=255)

                                  cGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGTCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                      <  2:385220/93‑1 (MQ=255)

                                  cGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                      <  2:2379550/93‑1 (MQ=255)

                                  cGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                      <  1:1523103/93‑1 (MQ=255)

                                  cGGATGGTACGGAGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                      <  2:4499932/93‑1 (MQ=255)

                                  cGGATGGGACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                      <  2:2379551/93‑1 (MQ=255)

                                   ggATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcat                                                     >  1:2129155/1‑93 (MQ=255)

                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatc                                                    >  1:362939/1‑93 (MQ=255)

                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatc                                                    >  2:4347958/1‑93 (MQ=255)

                                     aTGGTACGGCTAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatca                                                   <  2:3953968/93‑1 (MQ=255)

                                      tGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                  >  2:2455626/1‑93 (MQ=255)

                                      tGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                  <  1:2561262/93‑1 (MQ=255)

                                       ggTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAgg                                                 >  2:782667/1‑93 (MQ=255)

                                       ggTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAgg                                                 >  1:4527387/1‑93 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                <  2:3500552/93‑1 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                <  1:3202610/93‑1 (MQ=255)

                                         tACGGCGAATATTGCAGTGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCt                                               <  1:2244657/93‑1 (MQ=255)

                                         tACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATCCCTTCATCCggtc                                               >  1:1187465/1‑91 (MQ=255)

                                         tACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatcaaggt                                               >  2:2772887/1‑89 (MQ=255)

                                         tACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCt                                               <  1:1479544/93‑1 (MQ=255)

                                         tACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCt                                               >  1:2949996/1‑93 (MQ=255)

                                          aCGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTc                                              <  1:2444824/93‑1 (MQ=255)

                                            ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGAAGATACCTTCATCAGGCTCGc                                            >  1:4688473/1‑93 (MQ=255)

                                             gCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGGATACCAAAGCGGCTGGTGACGATACCATCCTCAAGCTTGcc                                           >  2:1729815/1‑93 (MQ=255)

                                               gAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCCAAGCGGCTGATGACGATACCTTCATCATGCTCGCCaa                                         >  1:827070/1‑93 (MQ=255)

                                               gAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCaa                                         >  2:656949/1‑93 (MQ=255)

                                               gAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCaa                                         >  2:656950/1‑93 (MQ=255)

                                                  tatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACa                                      <  1:186118/93‑1 (MQ=255)

                                                  tatTGCAGCGGTGAACGTCCCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACa                                      <  2:3489967/93‑1 (MQ=255)

                                                   atTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAg                                     >  2:489979/1‑93 (MQ=255)

                                                   atTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAg                                     >  2:1717005/1‑93 (MQ=255)

                                                    ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                    <  1:4036264/93‑1 (MQ=255)

                                                    ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                    <  1:4036103/93‑1 (MQ=255)

                                                    ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                    >  2:1376595/1‑93 (MQ=255)

                                                    ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                    >  2:2781362/1‑93 (MQ=255)

                                                      tcAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAtt                                  <  2:1257833/92‑1 (MQ=255)

                                                         gCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtc                               <  1:2865305/93‑1 (MQ=255)

                                                           ggTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgt                             >  2:4580226/1‑93 (MQ=255)

                                                            gTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                            <  1:2639431/93‑1 (MQ=255)

                                                            gTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                            >  1:58964/1‑93 (MQ=255)

                                                            gTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                            <  1:2901490/93‑1 (MQ=255)

                                                              gAACGTCGCGATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTcgtcgt                          <  1:1785220/93‑1 (MQ=255)

                                                              gAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTcgtcgt                          <  2:1061848/93‑1 (MQ=255)

                                                               aaCGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTa                         >  2:4036658/1‑93 (MQ=255)

                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                       <  1:209714/93‑1 (MQ=255)

                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                       >  1:3413317/1‑93 (MQ=255)

                                                                 cGTCCCAATCTGCGGATTTCTCATCGGAGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGACAAACAGATTGTCGTCGTAGt                       <  2:1602246/93‑1 (MQ=255)

                                                                  gTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAgcc                      >  1:3268059/1‑91 (MQ=255)

                                                                  gTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTc                      >  1:3268060/1‑93 (MQ=255)

                                                                  gTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTc                      >  2:1042765/1‑93 (MQ=255)

                                                                   tCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCg                     >  2:4309486/1‑93 (MQ=255)

                                                                    cGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCgg                    <  1:4534584/93‑1 (MQ=255)

                                                                    cGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGAGACCTGCATCAGGCTCGCCAAACAGATTGTCGGCGTAGTCgg                    >  2:4222649/1‑93 (MQ=255)

                                                                      ccATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGc                  <  2:2793031/93‑1 (MQ=255)

                                                                       cATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCt                 >  1:2508814/1‑93 (MQ=255)

                                                                       cATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCt                 >  1:2744118/1‑93 (MQ=255)

                                                                       cATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCt                 <  2:555728/93‑1 (MQ=255)

                                                                       cATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTACTCGGGCt                 >  1:3751726/1‑93 (MQ=255)

                                                                        aTCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCtt                >  1:2445801/1‑93 (MQ=255)

                                                                        aTCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCtt                >  1:4644522/1‑93 (MQ=255)

                                                                        aTCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTAGCCAAACAGATTGTCGTCGTAGTCGGGCtt                >  2:2873120/1‑93 (MQ=255)

                                                                        aTCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGACCGCCAAACAGATTGTCGTCGTAGTCGGGGtt                >  2:3982226/1‑93 (MQ=255)

                                                                         tCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtc               <  1:1742168/93‑1 (MQ=255)

                                                                         tCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtc               >  1:1575185/1‑93 (MQ=255)

                                                                           ggcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtctc             >  1:1542984/1‑93 (MQ=255)

                                                                              ggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCtt          <  2:287864/93‑1 (MQ=255)

                                                                                  tCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGAc      <  1:535353/93‑1 (MQ=255)

                                                                                    cacaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGt    >  1:2166568/1‑93 (MQ=255)

                                                                                    cacaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCAGGCTTCTCCTTAGACGt    <  2:33166/93‑1 (MQ=255)

                                                                                      caTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGttt  >  2:1493769/1‑93 (MQ=255)

ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTT  >  NC_000913/4388842‑4389020

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