BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,064,161	G→A	16.1%	L42L (CTG→CTA) ‡	yccE →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,064,161	0	G	A	16.1%	345.2 / 53.2	155	L42L (CTG→CTA)	yccE	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (54/75); new base A (17/8); total (71/84)
Fisher's exact test for biased strand distribution p-value = 2.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

Only 100 of 161 total aligned reads displayed.

TGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAA  >  NC_000913/1064063‑1064260

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

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

 gCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgt                                                                                                   >  1:2212083/1‑100 (MQ=255)

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

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

   aCTGCAAATAATTATTTGAAACAGGCCTGGAACTATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACAATGgtg                                                                                                  >  2:2482733/1‑99 (MQ=255)

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

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

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

      gCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGtt                                                                                              >  1:2393318/1‑100 (MQ=255)

      gCAAAGAGTTGTCTGAAAAGGGCCTGGAAGGATAGAAAAAATGAGTACCAAAAAAATCAAACATATTCAAGCACGCTCTGTGAAAACACACTGGTGTGtt                                                                                              <  2:1218040/100‑1 (MQ=255)

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

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

          ccgatcttttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttca                                                                                          <  1:2859015/93‑3 (MQ=255)

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

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

           taatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttcaa                                                                                         <  1:2472650/100‑4 (MQ=255)

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

            aatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGc                                                                                        <  1:2735615/100‑1 (MQ=255)

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

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

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

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

               tattTGAAACAGGCGTGGAACGATATAAAAAATGAGGACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGt                                                                                     <  2:1418655/100‑1 (MQ=255)

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

                attTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACAAACTGATGAGTTTTATGCgg                                                                                      >  2:3496538/1‑98 (MQ=255)

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

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

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

                       aCAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAAt                                                                             <  1:3118232/100‑1 (MQ=255)

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

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

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

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

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

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

                                   cGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgt                                                                 >  1:2603284/1‑100 (MQ=255)

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

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

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

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

                                            aaaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAAt                                                        <  1:3346436/100‑1 (MQ=255)

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

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

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

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

                                                    accaaaaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGa                                                <  2:3603625/97‑1 (MQ=255)

                                                    aaaaaaaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGa                                                <  2:2820551/97‑1 (MQ=255)

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

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

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

                                                           aaaTCAAACATATTCAATCACGCTTCTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCa                                         >  1:4105059/1‑100 (MQ=255)

                                                            aaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCAt                                        >  1:178175/1‑100 (MQ=255)

                                                             aTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATg                                       <  1:757584/100‑1 (MQ=255)

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

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

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

                                                                  aCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTgg                                  <  2:1582095/100‑1 (MQ=255)

                                                                   cATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTgg                                  <  2:3380203/99‑1 (MQ=255)

                                                                       ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                             >  1:1816369/1‑100 (MQ=255)

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

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

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

                                                                         cAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGa                           >  2:2048851/1‑100 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                                                    tttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa                >  2:1822371/1‑100 (MQ=255)

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

                                                                                    gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGCACCCAGACGTAAAGTAAATGAAGAGGGTACTGCATGTCTGGAATGTGAATAAcaaaaa                 >  2:3802992/5‑94 (MQ=255)

                                                                                    gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa                >  1:3147419/5‑100 (MQ=255)

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

                                                                                     ttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa               >  1:3901150/1‑100 (MQ=255)

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

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

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

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

                                                                                     cataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa               >  2:1649644/4‑100 (MQ=255)

                                                                                      tGAAAACACACGGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaaa              <  1:398500/100‑1 (MQ=255)

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

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

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

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

                                                                                            cacacTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAAttt        <  1:3682922/100‑1 (MQ=255)

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

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

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

                                                                                                cTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAgg    >  2:3897867/1‑100 (MQ=255)

                                                                                                cTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAgg    >  1:2317178/1‑100 (MQ=255)

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

                                                                                                  ggTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGaa  >  2:2559535/1‑100 (MQ=255)

                                                                                                  ggTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGaa  >  1:1579013/1‑100 (MQ=255)

TGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAA  >  NC_000913/1064063‑1064260

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