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.8%	2210.6 / ‑1.1	462	coding (1388/1458 nt)	b3876	putative permease
Reads supporting (aligned to +/- strand): ref base C (0/0); major base . (267/193); minor base T (1/0); total (269/193)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 510 total aligned reads displayed.

ATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAA  >  NC_000913/4061864‑4062050

  cTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccg                                                                                           <  2:9296471/96‑2 (MQ=255)

  cTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccg                                                                                           >  1:736150/1‑95 (MQ=255)

  cTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGCCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccg                                                                                           <  2:3271078/96‑2 (MQ=255)

   tttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTCTTGACATTTATAACGTTACGGAAGccgt                                                                                          >  2:9006107/1‑94 (MQ=255)

    ttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgtt                                                                                         >  2:10834955/1‑93 (MQ=255)

    ttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgtt                                                                                         >  1:1676153/1‑93 (MQ=255)

       tCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGttttt                                                                                     >  1:2438337/1‑96 (MQ=255)

        cATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGtttttt                                                                                    <  2:10344317/96‑1 (MQ=255)

        aataataAGTCCATTTTTTACACGTCGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGGTATTGACATTTATAACGTTACGGAAG‑CCGtttttt                                                                                    <  1:4658278/95‑1 (MQ=255)

         ataataAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttc                                                                                   >  2:5863617/1‑96 (MQ=255)

          taataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttct                                                                                  >  2:11576260/1‑96 (MQ=255)

          taataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttct                                                                                  >  2:6430862/1‑96 (MQ=255)

           aataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttctt                                                                                 >  1:3297220/1‑96 (MQ=255)

           aataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttctt                                                                                 >  1:12640323/1‑96 (MQ=255)

             taaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                               <  2:11310550/96‑1 (MQ=255)

             taaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                               <  2:13916289/96‑1 (MQ=255)

             taaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                               <  1:7623223/96‑1 (MQ=255)

              aaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGTATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                              <  1:7924473/96‑1 (MQ=255)

              aaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                              >  2:6937522/1‑96 (MQ=255)

                gTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTc                                                                            >  1:10750290/1‑96 (MQ=255)

                 tCCATTTTTTAAACTTAGGGAGGCCGCCCCCCCCCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCg                                                                           <  1:4400189/96‑1 (MQ=255)

                   cATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTAGcc                                                                         <  1:6242476/96‑1 (MQ=255)

                      tttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtc                                                                      >  2:13988339/1‑96 (MQ=255)

                      tttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑ACGTTTTTTCTTCTCTCGCCAtc                                                                      >  2:6609084/1‑96 (MQ=255)

                      tgtttAACCGTAGGGAGGCCGCCACCCGCCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtc                                                                      <  1:2884608/94‑1 (MQ=255)

                        tttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtctc                                                                    <  1:11358401/96‑1 (MQ=255)

                         ttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATctct                                                                   >  1:12183031/1‑96 (MQ=255)

                          tAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCtt                                                                  >  2:8356190/1‑96 (MQ=255)

                             aCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAt                                                               >  1:8854019/1‑96 (MQ=255)

                              cGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtt                                                              >  2:7173834/1‑96 (MQ=255)

                              cGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtt                                                              >  1:11110411/1‑96 (MQ=255)

                                tAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtttt                                                            >  1:6729408/1‑96 (MQ=255)

                                 aGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttttt                                                           >  2:2147337/1‑96 (MQ=255)

                                 aGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttttt                                                           <  2:4263793/96‑1 (MQ=255)

                                  gggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTg                                                          <  1:5863617/96‑1 (MQ=255)

                                   ggAGGCCGCCACCCACCAGGATGGAAATAGCGGGTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGc                                                         >  2:5516611/1‑96 (MQ=255)

                                    gAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCg                                                        >  1:13392378/1‑96 (MQ=255)

                                     aGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGATTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgt                                                       >  2:113728/1‑96 (MQ=255)

                                       gccgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttg                                                     >  2:3482787/1‑96 (MQ=255)

                                       gccgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgtcg                                                     <  1:11887979/96‑3 (MQ=255)

                                         cgccACCCACCAGGATGGAAATAACGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggt                                                   >  2:13669561/1‑96 (MQ=255)

                                          gccCCCCTCCAGTATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggct                                                  <  1:11837128/96‑3 (MQ=255)

                                          gccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggtt                                                  >  1:9808548/1‑96 (MQ=255)

                                          gccACCAACCAGGATGGAAATAGCGGCTATTGACATATATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggtt                                                  <  2:9417952/96‑1 (MQ=255)

                                           ccacccacCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggttg                                                 <  2:8159818/96‑1 (MQ=255)

                                               ccaccaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAttt                                             >  2:11165818/1‑96 (MQ=255)

                                                 accaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTcc                                           >  2:5693096/1‑96 (MQ=255)

                                                    gggATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                        >  1:11684880/2‑96 (MQ=255)

                                                    aagATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                        >  1:3565634/3‑96 (MQ=255)

                                                    aGGATGGAAATAGCGGCTATTGACATTTATACCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                        <  2:1285363/96‑1 (MQ=255)

                                                    aGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCa                                        >  2:13451048/1‑96 (MQ=255)

                                                     ggATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaa                                       <  2:1590705/96‑1 (MQ=255)

                                                     ggATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaa                                       >  1:8709769/1‑96 (MQ=255)

                                                      gATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaat                                      <  2:11900634/96‑1 (MQ=255)

                                                       aTGGAAATAGCGGCTTTTGACATTTATAACGGTACGGAAG‑CCGTTTGTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaata                                     >  2:1690449/1‑96 (MQ=255)

                                                        tGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataa                                    >  2:12153921/1‑96 (MQ=255)

                                                         ggAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataat                                   <  2:8773597/96‑1 (MQ=255)

                                                          gAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATc                                  <  2:14154544/96‑1 (MQ=255)

                                                          caaaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATc                                  <  1:2147337/95‑1 (MQ=255)

                                                           aaaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCa                                 <  2:27858/96‑1 (MQ=255)

                                                             aTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAAt                               <  1:4129860/96‑1 (MQ=255)

                                                             aTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAAt                               >  2:3689619/1‑96 (MQ=255)

                                                               aGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGc                             <  1:10013913/96‑1 (MQ=255)

                                                               aGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGc                             >  1:263857/1‑96 (MQ=255)

                                                                gCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCa                            >  2:5165943/1‑96 (MQ=255)

                                                                 cGGCTATTGACATTTATAATGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaa                           <  1:14368560/96‑1 (MQ=255)

                                                                 cGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaa                           >  1:1851549/1‑96 (MQ=255)

                                                                  ggCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaaa                          >  1:9513071/1‑96 (MQ=255)

                                                                    caaTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                        >  1:9195445/3‑96 (MQ=255)

                                                                    cTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                        >  2:381272/1‑96 (MQ=255)

                                                                    cTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                        <  2:4855706/96‑1 (MQ=255)

                                                                     tATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCg                       <  2:10084624/96‑1 (MQ=255)

                                                                     tATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCg                       <  2:13058718/96‑1 (MQ=255)

                                                                     tATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTAAATAATCAATGCAAACCg                       >  1:13161056/1‑96 (MQ=255)

                                                                      aTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGt                      <  1:8131521/96‑1 (MQ=255)

                                                                      aTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGt                      <  2:6157106/96‑1 (MQ=255)

                                                                       ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtt                     >  2:14267005/1‑96 (MQ=255)

                                                                       ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtt                     <  1:14060242/96‑1 (MQ=255)

                                                                       ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtt                     >  2:4486755/1‑96 (MQ=255)

                                                                         gACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                   >  1:8475017/1‑96 (MQ=255)

                                                                         gACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                   >  1:3780534/1‑96 (MQ=255)

                                                                         gACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                   >  1:932087/1‑96 (MQ=255)

                                                                            agttATAAGGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCa                >  1:10502053/3‑96 (MQ=255)

                                                                            aTTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCa                >  2:9649332/1‑96 (MQ=255)

                                                                            aTTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCa                >  1:7118251/1‑96 (MQ=255)

                                                                              ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCACACCGTTTTTCAtt              <  1:13464904/96‑1 (MQ=255)

                                                                              ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt              <  2:10300172/96‑1 (MQ=255)

                                                                              ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt              <  1:1383274/96‑1 (MQ=255)

                                                                              ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt              >  1:6005295/1‑96 (MQ=255)

                                                                              ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt              >  1:13985468/1‑96 (MQ=255)

                                                                                ataACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCGTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCCAACCGTTTTTCATTGa            >  2:7958121/1‑96 (MQ=255)

                                                                                 taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAg           >  1:5015910/1‑96 (MQ=255)

                                                                                  aaCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgt          <  1:6430862/96‑1 (MQ=255)

                                                                                  aaCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgt          <  2:2466096/96‑1 (MQ=255)

                                                                                  aaCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGGTTTTCATTGAgt          >  2:7608982/1‑96 (MQ=255)

                                                                                   tcgatACGTAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtg         <  1:5440251/92‑1 (MQ=255)

                                                                                     gTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgtg       <  1:5701745/96‑1 (MQ=255)

                                                                                     gTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgtg       <  2:12411857/96‑1 (MQ=255)

                                                                                       tACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTa     <  2:11801633/96‑1 (MQ=255)

                                                                                        aCGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTaa    <  1:6942100/96‑1 (MQ=255)

ATCTTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTAAAA  >  NC_000913/4061864‑4062050

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