BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,064,160	T→C	16.2%	L42P (CTG→CCG) ‡	yccE →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,064,160	0	T	C	16.2%	334.9 / 56.9	157	L42P (CTG→CCG)	yccE	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base T (57/73); new base C (17/8); total (74/83)
Fisher's exact test for biased strand distribution p-value = 3.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-02

Only 100 of 160 total aligned reads displayed.

GTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGA  >  NC_000913/1064061‑1064259

  tGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTgg                                                                                                   >  2:2308030/1‑100 (MQ=255)

  tGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTgg                                                                                                   <  1:1739960/100‑1 (MQ=255)

   gCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACTAAAAAAATAAAACATATTCAATCACGCTTTTTGAAAAAACAATg                                                                                                    >  2:2451645/1‑98 (MQ=255)

    cACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtg                                                                                                 >  1:519816/1‑100 (MQ=255)

    cACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtg                                                                                                 <  2:4062263/100‑1 (MQ=255)

      cTGCAAATAATTATTTGAAACAGGCTTGCAACGATATAAAAAATGAAAAAGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtgtg                                                                                               <  2:3697362/100‑1 (MQ=255)

      cTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtgtg                                                                                               >  1:145570/1‑100 (MQ=255)

      cTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtgtg                                                                                               >  1:4095756/1‑100 (MQ=255)

        gCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtgtgt                                                                                              >  1:2409023/1‑99 (MQ=255)

         cAAAAAATTATTTGAAACAGGCCTGGAACGATATAAAAAATAAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGttt                                                                                            <  2:1008910/100‑1 (MQ=255)

          aaaTAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGtttt                                                                                           >  1:209689/1‑100 (MQ=255)

           aataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTa                                                                                          <  1:4087710/100‑1 (MQ=255)

            ataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttca                                                                                         <  1:1086911/100‑3 (MQ=255)

            ataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttca                                                                                         <  1:2371187/100‑3 (MQ=255)

            ataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttca                                                                                         <  1:2371104/100‑3 (MQ=255)

            ataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTAAATCACGCTTTTTGAAAACACACCAGTGTGttttca                                                                                         <  2:1947325/100‑3 (MQ=255)

              gattattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGc                                                                                       >  1:55577/2‑100 (MQ=255)

               attattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCg                                                                                      >  2:3085500/1‑100 (MQ=255)

                ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                     <  2:519811/100‑1 (MQ=255)

                ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                     >  1:1884194/1‑100 (MQ=255)

                 tattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGt                                                                                    <  1:3298007/100‑1 (MQ=255)

                 tattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttcaaaaag                                                                                    <  2:1802105/100‑8 (MQ=38)

                  attTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGtt                                                                                   >  1:3481029/1‑100 (MQ=255)

                  attTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGtt                                                                                   >  1:3026292/1‑100 (MQ=255)

                  attTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGtt                                                                                   >  2:3704990/1‑100 (MQ=255)

                   tttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTta                                                                                  >  2:2714203/1‑100 (MQ=255)

                     tGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTtata                                                                                <  1:3360495/100‑1 (MQ=255)

                      gAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATAc                                                                               <  1:3970892/100‑1 (MQ=255)

                      gAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATAc                                                                               <  2:1598604/100‑1 (MQ=255)

                      gAAACACGCCTTGAACGAGATAAAAAATGCGTACTAAAAAAATGAAACATATTCAATCACTCTTTGTTAAAACACTCTGGTGTGTTTTATGCTGTtatat                                                                               <  2:556246/100‑2 (MQ=255)

                       aaaCAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACa                                                                              <  1:729835/100‑1 (MQ=255)

                           aGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGa                                                                          <  2:607271/100‑1 (MQ=255)

                              ccTGGAACGATATAAAAAATGAGTACGAAAAAAATCGAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACt                                                                       >  1:2915749/1‑100 (MQ=255)

                              ccTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACt                                                                       >  1:872084/1‑100 (MQ=255)

                               cTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTc                                                                      <  1:3972271/100‑1 (MQ=255)

                                 ggAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAg                                                                    >  1:3165629/1‑100 (MQ=255)

                                 ggAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAg                                                                    <  2:2821993/100‑1 (MQ=255)

                                   aaCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGAc                                                                  <  2:633760/100‑1 (MQ=255)

                                   aaCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGAc                                                                  >  1:2636267/1‑100 (MQ=255)

                                    aCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACg                                                                 <  1:1298774/100‑1 (MQ=255)

                                    aCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACg                                                                 >  1:3829128/1‑100 (MQ=255)

                                       atttaAAAAAAGAGTAAAAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgta                                                               <  2:1304541/96‑1 (MQ=255)

                                       atatAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaa                                                              <  2:169410/100‑1 (MQ=255)

                                        tataAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaa                                                             <  1:1743607/100‑1 (MQ=255)

                                         ataAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaag                                                            >  1:2108655/1‑100 (MQ=255)

                                          taAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagt                                                           >  1:1647870/1‑100 (MQ=255)

                                          taAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagt                                                           >  2:743967/1‑100 (MQ=255)

                                           aaaaaaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagta                                                          >  1:1915407/1‑100 (MQ=255)

                                           aaaaaaTGAGTAAGAAAAAAATCAAACATGTTCAATCACTCTGTTTGCAAACACACTGGTGTGTTTTGTGTGGTTATACAATGAACGCATAGgtaaagta                                                          <  2:2910545/100‑1 (MQ=255)

                                               aaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTTTTTTTTTTTTGGTGTGTTTTATGCGTTTATTCTTTTAACTTAGACGTAGAGTAAATg                                                      >  2:1320937/1‑100 (MQ=255)

                                                 tGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGaa                                                    >  1:1813450/1‑100 (MQ=255)

                                                   aGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAaga                                                  >  2:3195180/1‑100 (MQ=255)

                                                    gTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAagag                                                 <  1:2226557/100‑1 (MQ=255)

                                                    gTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAagag                                                 <  1:1531150/100‑1 (MQ=255)

                                                      aCGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGa                                               <  2:3887028/100‑1 (MQ=255)

                                                      aCGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGa                                               <  1:2598406/100‑1 (MQ=255)

                                                         aaaaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGTTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATAAAGAgg                                                >  2:2253750/1‑96 (MQ=255)

                                                          aaaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGCTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATa                                             >  2:1820492/1‑98 (MQ=255)

                                                             aaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCa                                        <  2:3901813/100‑1 (MQ=255)

                                                              aaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCa                                        <  2:3572111/99‑1 (MQ=255)

                                                               aTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATg                                      >  2:3998362/1‑100 (MQ=255)

                                                                 cAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTc                                    >  2:2355608/1‑100 (MQ=255)

                                                                  aaaCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCt                                   <  1:959377/100‑1 (MQ=255)

                                                                   aaCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTg                                  <  1:1485144/100‑1 (MQ=255)

                                                                   aaCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTg                                  <  1:2275195/100‑1 (MQ=255)

                                                                         ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                            <  1:359787/100‑1 (MQ=255)

                                                                         ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                            <  1:2763301/100‑1 (MQ=255)

                                                                          tCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgtg                           <  1:2153175/100‑1 (MQ=255)

                                                                          tCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgtg                           <  1:264853/100‑1 (MQ=255)

                                                                           cAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgtg                           <  1:555471/99‑1 (MQ=255)

                                                                           cAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGa                          <  2:3517201/100‑1 (MQ=255)

                                                                            aaTCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGaa                         <  2:3411725/100‑1 (MQ=255)

                                                                             aTCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAAt                        >  1:386385/1‑100 (MQ=255)

                                                                             aTCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAAt                        >  1:741310/1‑100 (MQ=255)

                                                                              tCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATc                       <  1:3129555/100‑1 (MQ=255)

                                                                              tCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATc                       <  2:2750028/100‑1 (MQ=255)

                                                                                 cGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACt                    >  1:3910995/1‑100 (MQ=255)

                                                                                 cGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACt                    >  2:1618896/1‑100 (MQ=255)

                                                                                 cGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACt                    <  2:754514/100‑1 (MQ=255)

                                                                                 cGCTTTTTGAAAACACACTGGTGTGTTTTAGGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACt                    <  2:3312635/100‑1 (MQ=255)

                                                                                  gCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTa                   <  1:3376618/100‑1 (MQ=255)

                                                                                     cgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaa                >  1:2475513/6‑100 (MQ=255)

                                                                                     cgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaa                >  2:3466346/6‑100 (MQ=255)

                                                                                     cgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATAACTAGaa                >  2:368267/6‑100 (MQ=255)

                                                                                     cgcataaaaCACACCAGTGTGTTGTATGCGGTTATACAATGAAATCAGACGTAAAGTAAATGAAGAGGCTACTCCATATCTGGAATATGAATCAc                     >  2:1539167/6‑95 (MQ=255)

                                                                                      gcataaaaCACACTGGTGTGTTTTATGCGGTTATACAACGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa               >  2:2540002/5‑100 (MQ=255)

                                                                                      gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa               >  2:2525109/5‑100 (MQ=255)

                                                                                      gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACATAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa               >  2:2032108/5‑100 (MQ=255)

                                                                                       cataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa              >  1:3081766/4‑100 (MQ=255)

                                                                                       cataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa              >  1:1014291/4‑100 (MQ=255)

                                                                                        tGAAAATACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAAACGTAAAGTAAATGAAGTGGATACTCCATGTCTGGAATGTTAAACACTAGaaaaa             <  2:1772977/100‑1 (MQ=255)

                                                                                         gAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAg            <  1:2895536/100‑1 (MQ=255)

                                                                                           aaaCACACTGGTATGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTAAAAAAATAAAAAAAg            <  2:2140425/98‑1 (MQ=255)

                                                                                            aaCACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAAt         <  2:678229/100‑1 (MQ=255)

                                                                                             acacacTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAAtt        <  2:338246/100‑1 (MQ=255)

                                                                                             acacacATGTGAGTTTTATGCAGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTTGAATGTGAATCACTAGAAAAAGAAt         <  2:1968305/99‑1 (MQ=255)

                                                                                               acacTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTg      <  1:36303/100‑1 (MQ=255)

                                                                                                 acTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAg    <  2:3842562/100‑1 (MQ=255)

                                                                                                 acTGGTGTGTTTTATGCGGTTATACAATGAACGCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAg    <  2:887195/100‑1 (MQ=255)

                                                                                                   tGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGa  <  2:166362/100‑1 (MQ=255)

GTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGA  >  NC_000913/1064061‑1064259

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