BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,061,961	Δ1 bp	coding (1386/1458 nt)	b3876 ←	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,061,959	0	C	.	99.5%	2555.9 / 1.0	567	coding (1388/1458 nt)	b3876	putative permease
Reads supporting (aligned to +/- strand): ref base C (2/1); new base . (315/249); total (317/250)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 615 total aligned reads displayed.

TTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTA  >  NC_000913/4061867‑4062047

   ttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgt                                                                                       >  1:2245521/1‑91 (MQ=255)

    tCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgtt                                                                                      <  1:8067965/93‑4 (MQ=255)

    tCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgtt                                                                                      >  2:4228611/1‑90 (MQ=255)

      ataataAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGtttt                                                                                   >  1:7566443/1‑93 (MQ=255)

        aataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGtttttt                                                                                 <  2:1476617/93‑1 (MQ=255)

        aataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGtttttt                                                                                 <  2:8034734/93‑1 (MQ=255)

         ataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttc                                                                                <  2:4663991/93‑1 (MQ=255)

          taaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttct                                                                               >  1:488227/1‑93 (MQ=255)

           aaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttctt                                                                              >  2:7043174/1‑93 (MQ=255)

            aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttcttc                                                                             <  1:6092153/93‑1 (MQ=255)

            aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttcttc                                                                             <  1:632385/93‑1 (MQ=255)

            aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttcttc                                                                             <  1:2436400/93‑1 (MQ=255)

             gTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                            <  2:6488366/93‑1 (MQ=255)

              tCCATTTTTTAAACGTAGTTAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGCAG‑CCGTTTTTTCTTCTg                                                                           >  2:3753723/1‑93 (MQ=255)

              tCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                           >  2:3929908/1‑93 (MQ=255)

              tCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                           <  1:9059592/93‑1 (MQ=255)

              tCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGAAATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                           >  2:2337389/1‑93 (MQ=255)

                cATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTc                                                                         >  2:3480375/1‑93 (MQ=255)

                cATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACACTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTc                                                                         <  2:6378568/93‑1 (MQ=255)

                 aTTTTTTAAACGTAGGGAGGCCGCCACCCCCCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCg                                                                        >  2:7706010/1‑93 (MQ=255)

                 aTTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCg                                                                        <  2:9229209/93‑1 (MQ=255)

                 aTTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGACGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTgtgtag                                                                        >  1:2856049/1‑91 (MQ=255)

                  ttttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGCAG‑CCGTTTTTTCTTCTGTGGc                                                                       >  1:978246/1‑93 (MQ=255)

                  ttttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGc                                                                       >  2:9553822/1‑93 (MQ=255)

                  ttttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGc                                                                       >  2:8284115/1‑93 (MQ=255)

                  ttttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGc                                                                       >  1:3324834/1‑93 (MQ=255)

                   tttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGcc                                                                      >  2:7316145/1‑93 (MQ=255)

                    ttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCa                                                                     <  1:3711060/93‑1 (MQ=255)

                       tAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtct                                                                  >  1:3907068/1‑93 (MQ=255)

                           cGTAGGGAGGCCGCCACCCAACAGGATGGAAATAGCGGCTATTGACATTTATCACGTTACGGAAG‑CCGTCTTTTCTTCTGTCGCCATGTCttt                                                              >  1:5036229/1‑93 (MQ=255)

                            gTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTa                                                             >  2:6787101/1‑93 (MQ=255)

                            gTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTa                                                             <  1:4774831/93‑1 (MQ=255)

                            gTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTa                                                             >  2:1616904/1‑93 (MQ=255)

                              aGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtt                                                           <  1:6040156/93‑1 (MQ=255)

                              aGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtt                                                           >  1:7674289/1‑93 (MQ=255)

                               gggAGGCCGCCACCCACCAGGATGGAAATAGAGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGACATCTCTTTAttt                                                          >  2:1702868/1‑93 (MQ=255)

                                ggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtttt                                                         <  2:3321659/93‑1 (MQ=255)

                                 gAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttttt                                                        <  2:5476835/93‑1 (MQ=255)

                                    gccgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCg                                                     >  2:2238824/1‑93 (MQ=255)

                                     ccgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgt                                                    <  2:5730161/93‑1 (MQ=255)

                                       gccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttg                                                  >  2:862174/1‑93 (MQ=255)

                                          acccaccAGGATGGAAATAGCGGCTATTGACATTTTTAACGTTACGGAAG‑CCGATTTTTCTTCTGGCGCCATCTCTTTATTTTTCCgttggtt                                               >  2:6242858/1‑93 (MQ=255)

                                             caccaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAt                                            >  1:4164653/1‑93 (MQ=255)

                                                caGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTc                                         <  1:5385022/93‑1 (MQ=255)

                                                caGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTc                                         >  1:8701984/1‑93 (MQ=255)

                                                 aGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTcc                                        >  1:9208664/1‑93 (MQ=255)

                                                  ggATGGAAATAGCGTCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCt                                       <  2:7860579/93‑1 (MQ=255)

                                                  ggATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCt                                       >  2:2085014/1‑93 (MQ=255)

                                                  ggATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCt                                       >  1:7757697/1‑93 (MQ=255)

                                                    aTGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                     >  2:2989887/1‑93 (MQ=255)

                                                    aTGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                     >  1:4937320/1‑93 (MQ=255)

                                                    aTGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                     <  1:7492592/93‑1 (MQ=255)

                                                    aTGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                     >  1:4516751/1‑93 (MQ=255)

                                                    aTGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                     <  1:9450341/93‑1 (MQ=255)

                                                     tGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaa                                    <  1:3815166/93‑1 (MQ=255)

                                                      ggAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaat                                   >  2:4028733/1‑93 (MQ=255)

                                                       gAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaata                                  >  1:3552428/1‑93 (MQ=255)

                                                       gAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaata                                  >  1:8398545/1‑93 (MQ=255)

                                                        aaaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataa                                 <  2:5869888/93‑1 (MQ=255)

                                                        aaaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataa                                 <  1:3107706/93‑1 (MQ=255)

                                                          aTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATc                               >  2:2491834/1‑93 (MQ=255)

                                                           tAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCa                              <  2:4542065/93‑1 (MQ=255)

                                                            aGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCaa                             <  1:5196360/93‑1 (MQ=255)

                                                             gCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTTTCGCCATCTCTTTATTTTTGCGTTGGTTTATTTCCTCAATCATCAAt                            >  2:7504058/1‑93 (MQ=255)

                                                             gCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAAt                            <  2:4768857/93‑1 (MQ=255)

                                                              gggCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATg                           <  1:695906/92‑1 (MQ=255)

                                                               ggCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGc                          >  1:2030996/1‑93 (MQ=255)

                                                                gCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCa                         <  2:2807196/93‑1 (MQ=255)

                                                                gCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCa                         <  1:8338602/93‑1 (MQ=255)

                                                                   aTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAc                      <  2:4462285/93‑1 (MQ=255)

                                                                   aTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAc                      >  1:182615/1‑93 (MQ=255)

                                                                   aTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAc                      <  1:7344269/93‑1 (MQ=255)

                                                                   aTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAc                      >  2:5075703/1‑93 (MQ=255)

                                                                   aTTGACATTTATAACGTTACGGAAG‑CCGCTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATCATCAATGCAAAc                      >  1:1387204/1‑93 (MQ=255)

                                                                    ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                     >  1:4038765/1‑93 (MQ=255)

                                                                    ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                     >  2:9009675/1‑93 (MQ=255)

                                                                    ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                     >  2:887209/1‑93 (MQ=255)

                                                                    ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                     >  1:66964/1‑93 (MQ=255)

                                                                     tGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCg                    <  2:9636239/93‑1 (MQ=255)

                                                                     tGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCg                    <  2:3131191/93‑1 (MQ=255)

                                                                      gACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGt                   <  2:2700328/93‑1 (MQ=255)

                                                                       aCATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtt                  <  1:3475237/93‑1 (MQ=255)

                                                                        cATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttt                 >  2:5554240/1‑93 (MQ=255)

                                                                         aTTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                <  1:5711112/93‑1 (MQ=255)

                                                                          tttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttttt               <  1:4268300/93‑1 (MQ=255)

                                                                           ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc              <  1:5507495/93‑1 (MQ=255)

                                                                           ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc              >  2:8417208/1‑93 (MQ=255)

                                                                            tataACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCa             <  2:1211720/93‑1 (MQ=255)

                                                                             ataACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAt            <  1:1103716/93‑1 (MQ=255)

                                                                             ataACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAt            <  1:5076697/93‑1 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           <  2:7404224/93‑1 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           <  1:9612016/93‑1 (MQ=255)

                                                                               aaCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTg          <  1:9942941/93‑1 (MQ=255)

                                                                                aCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGa         <  2:664207/93‑1 (MQ=255)

                                                                                aCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGa         <  1:1167383/93‑1 (MQ=255)

                                                                                aCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGa         >  2:5232037/1‑93 (MQ=255)

                                                                                  gTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgt       >  2:1800495/1‑93 (MQ=255)

                                                                                   ttACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtg      >  2:355026/1‑93 (MQ=255)

                                                                                   ttACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtg      <  2:489297/93‑1 (MQ=255)

                                                                                      cGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgtgt   <  1:6513680/93‑1 (MQ=255)

TTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTA  >  NC_000913/4061867‑4062047

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