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.6%	1153.2 / ‑3.7	245	coding (1388/1458 nt)	b3876	putative permease
Reads supporting (aligned to +/- strand): ref base C (1/0); new base . (146/98); total (147/98)
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 267 total aligned reads displayed.

TTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTA  >  NC_000913/4061867‑4062047

tttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGc                                                                                          >  1:6546308/1‑93 (MQ=255)

tttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGc                                                                                          >  1:5502453/1‑93 (MQ=255)

tttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGc                                                                                          >  2:4302683/1‑93 (MQ=255)

 ttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCGGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGcc                                                                                         <  1:1815313/93‑1 (MQ=255)

  tttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccg                                                                                        >  2:2659121/1‑92 (MQ=255)

   ttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgt                                                                                       <  2:3672203/93‑3 (MQ=255)

     cATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgttt                                                                                     >  1:753918/1‑89 (MQ=255)

     cATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgttt                                                                                     <  1:6097759/93‑5 (MQ=255)

     cATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccgttt                                                                                     >  1:2938211/1‑89 (MQ=255)

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

        aataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGtttttt                                                                                 >  2:309971/1‑93 (MQ=255)

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

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

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

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

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

          aaaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttct                                                                               <  1:6321317/92‑1 (MQ=255)

             gTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                            >  1:6603277/1‑93 (MQ=255)

              tcgaTTTTTTAAACGTAGGGAGGCCGCCACCGACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                           >  1:4202749/4‑93 (MQ=255)

              tCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttctt                                                                           >  2:2479216/1‑92 (MQ=255)

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

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

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

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

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

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

               ccATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGt                                                                          >  1:3491605/1‑93 (MQ=255)

               ccATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGt                                                                          >  2:5761882/1‑93 (MQ=255)

                 aTTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCg                                                                        >  1:5641999/1‑93 (MQ=255)

                 aTTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCg                                                                        >  1:4563550/1‑93 (MQ=255)

                 aTTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCg                                                                        >  2:1107626/1‑93 (MQ=255)

                     tttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAt                                                                    >  2:3762811/1‑93 (MQ=255)

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

                        aaaCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtctc                                                                 <  1:4266942/93‑1 (MQ=255)

                         aaCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATctct                                                                <  1:859868/93‑1 (MQ=255)

                          aCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCtt                                                               <  2:2079558/93‑1 (MQ=255)

                          aCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCtt                                                               >  2:1852130/1‑93 (MQ=255)

                          aCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCtt                                                               <  1:338022/93‑1 (MQ=255)

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

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

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

                             tAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAt                                                            >  2:2248816/1‑93 (MQ=255)

                             tAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAt                                                            >  1:6381815/1‑93 (MQ=255)

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

                                ggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtttt                                                         >  2:3843773/1‑93 (MQ=255)

                                 gAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttttt                                                        >  1:3656038/1‑93 (MQ=255)

                                 gAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttttt                                                        >  2:898992/1‑93 (MQ=255)

                                   ggCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGc                                                      <  2:5062291/93‑1 (MQ=255)

                                   ggCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGc                                                      >  1:2379531/1‑93 (MQ=255)

                                    gccgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCg                                                     <  1:1724078/93‑1 (MQ=255)

                                      cgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgtt                                                   <  2:2517765/93‑1 (MQ=255)

                                      cgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgtt                                                   >  2:2124549/1‑93 (MQ=255)

                                      cgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgtt                                                   >  1:2268542/1‑93 (MQ=255)

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

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

                                         cacccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggt                                                >  1:1283345/1‑93 (MQ=255)

                                         cacccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggt                                                >  2:293864/1‑93 (MQ=255)

                                          acccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggtt                                               >  2:4414012/1‑93 (MQ=255)

                                          acccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggtt                                               >  2:4732976/1‑93 (MQ=255)

                                             caccaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAt                                            <  2:3069918/93‑1 (MQ=255)

                                              accaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAtt                                           >  1:2480169/1‑93 (MQ=255)

                                              accaGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGAtt                                           <  2:522236/93‑1 (MQ=255)

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

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

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

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

                                                      ggAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaat                                   <  2:5188939/93‑1 (MQ=255)

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

                                                      ggAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaat                                   <  2:6285331/93‑1 (MQ=255)

                                                       gAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaata                                  <  1:4597248/93‑1 (MQ=255)

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

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

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

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

                                                               ggCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGc                          <  2:5946011/93‑1 (MQ=255)

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

                                                                  tATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCaaa                       >  2:359168/1‑93 (MQ=255)

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

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

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

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

                                                                       aCATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtt                  >  2:401997/1‑93 (MQ=255)

                                                                       aCATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtt                  >  1:1624639/1‑93 (MQ=255)

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

                                                                       aCATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtt                  >  2:1423877/1‑93 (MQ=255)

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

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

                                                                          tttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttttt               >  1:5992228/1‑93 (MQ=255)

                                                                             ataACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAt            >  2:971987/1‑93 (MQ=255)

                                                                             ataACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAt            >  2:5471584/1‑93 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCCATAATCAATGCAAACCGTTTTTCAtt           >  2:1199554/1‑93 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           >  2:1562946/1‑93 (MQ=255)

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

                                                                               aaCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTg          >  1:4037491/1‑93 (MQ=255)

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

                                                                               aaCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTg          >  1:4806982/1‑93 (MQ=255)

                                                                                   ttACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCCATGCAAACCGTTTTTCATTGAgtg      <  1:6404471/93‑1 (MQ=255)

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

                                                                                    tACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgt     <  1:3733100/93‑1 (MQ=255)

                                                                                     aCGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgtg    >  1:4700569/1‑93 (MQ=255)

TTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTA  >  NC_000913/4061867‑4062047

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