BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,013,177	+C	coding (543/999 nt)	ysgA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,013,177	1	.	C	99.6%	895.3 / ‑4.3	274	coding (543/999 nt)	ysgA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (137/136); total (137/137)
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 279 total aligned reads displayed.

TTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCCG‑GCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAC  >  NC_000913/4013092‑4013260

ttAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTa                                                                               >  2:6246541/1‑93 (MQ=255)

ttAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTa                                                                               >  2:2641634/1‑93 (MQ=255)

 tAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTaa                                                                              <  1:2938873/93‑1 (MQ=255)

   gcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGa                                                                            <  2:1099654/93‑1 (MQ=255)

   gcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGa                                                                            <  2:1696829/93‑1 (MQ=255)

   gcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGa                                                                            >  1:3942204/1‑93 (MQ=255)

   gcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGa                                                                            >  2:5255741/1‑93 (MQ=255)

   gcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGa                                                                            >  2:2944899/1‑93 (MQ=255)

   gcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGa                                                                            <  1:4216349/93‑1 (MQ=255)

    cagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGAt                                                                           <  2:3817023/93‑1 (MQ=255)

    cagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGAt                                                                           <  2:2092351/93‑1 (MQ=255)

     agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATc                                                                          <  2:3512928/93‑1 (MQ=255)

      gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCg                                                                         <  1:1899089/93‑1 (MQ=255)

      gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCg                                                                         >  1:869580/1‑93 (MQ=255)

       cccGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCgg                                                                        <  1:5603433/93‑1 (MQ=255)

       cccGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCgg                                                                        >  2:5531680/1‑93 (MQ=255)

        ccGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGt                                                                       <  2:1083955/93‑1 (MQ=255)

        ccGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGt                                                                       <  2:150475/93‑1 (MQ=255)

         cGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGtt                                                                      <  1:1408879/93‑1 (MQ=255)

          gCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTg                                                                     >  1:5841413/1‑93 (MQ=255)

          gCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTg                                                                     >  1:3188951/1‑93 (MQ=255)

            aGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCg                                                                   <  2:5984879/93‑1 (MQ=255)

             gcgcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGa                                                                  >  1:1544416/1‑93 (MQ=255)

             gcgcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGa                                                                  <  1:4518334/93‑1 (MQ=255)

              cgcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGat                                                                 <  1:1078137/93‑1 (MQ=255)

               gcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGata                                                                <  2:2424871/93‑1 (MQ=255)

                 cTGGGGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATc                                                              <  1:2259264/93‑1 (MQ=255)

                 cTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATc                                                              >  1:3618032/1‑93 (MQ=255)

                 cTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATc                                                              <  2:3699008/93‑1 (MQ=255)

                 cTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGACGAGAATCGGCGCGTTAAGATCGGTTGCGATATc                                                              >  1:279579/1‑93 (MQ=255)

                   ggCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCaa                                                            >  1:791502/1‑93 (MQ=255)

                     cgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACa                                                          >  2:6391146/1‑93 (MQ=255)

                       cATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAgg                                                        >  1:550225/1‑93 (MQ=255)

                        aTGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgt                                                       >  2:3771748/1‑93 (MQ=255)

                        aTGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgt                                                       >  2:1794868/1‑93 (MQ=255)

                         tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                      >  1:4348876/1‑93 (MQ=255)

                            tttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgttgtt                                                   <  1:4120620/93‑1 (MQ=255)

                            tttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgttgtt                                                   <  1:2365612/93‑1 (MQ=255)

                            tttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgttgtt                                                   <  1:3687721/93‑1 (MQ=255)

                             ttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCACCAGGTTGttt                                                  <  1:1934821/93‑1 (MQ=255)

                             ttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGttt                                                  <  1:5158579/93‑1 (MQ=255)

                                aaCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCgg                                               >  2:6605235/1‑93 (MQ=255)

                                aaCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCgc                                               <  2:4172941/93‑2 (MQ=255)

                                 aCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                              <  1:4767121/93‑1 (MQ=255)

                                 aCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                              >  2:5384198/1‑93 (MQ=255)

                                 aCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                              >  2:2585000/1‑93 (MQ=255)

                                 aCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                              <  1:2778029/93‑1 (MQ=255)

                                  cGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTg                                             >  2:5675895/1‑93 (MQ=255)

                                   gCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGa                                            >  2:2070628/1‑93 (MQ=255)

                                    ctctCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGaa                                           <  1:1694808/93‑1 (MQ=255)

                                    ctctCCTGCGGAATGCTGTTATCCTGACCACCATATAAGACGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGaa                                           >  1:6653209/1‑93 (MQ=255)

                                     tctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAt                                          >  1:5010335/1‑93 (MQ=255)

                                     tctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAt                                          >  2:4465053/1‑93 (MQ=255)

                                      ctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                         <  2:813655/93‑1 (MQ=255)

                                       tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTc                                        >  2:1221514/1‑93 (MQ=255)

                                       tcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTc                                        <  1:5280737/93‑1 (MQ=255)

                                        ccTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCa                                       >  1:3223949/1‑93 (MQ=255)

                                        ccTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCa                                       <  2:6309812/93‑1 (MQ=255)

                                         cTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAg                                      <  1:3277725/93‑1 (MQ=255)

                                          tGCGGAATGCTGTTATCCTGACCACCATATAAGACGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                                     >  2:4415341/1‑93 (MQ=255)

                                            cGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGTTGAATTCAGCGa                                   >  2:5816423/1‑93 (MQ=255)

                                            cGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGa                                   >  1:2563414/1‑93 (MQ=255)

                                               aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                                <  2:6176428/93‑1 (MQ=255)

                                                aTGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTg                               >  2:4522086/1‑93 (MQ=255)

                                                 tGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGt                              <  2:5130673/93‑1 (MQ=255)

                                                  tcTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc                             <  2:2046392/92‑1 (MQ=255)

                                                  gCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc                             >  2:2133590/1‑93 (MQ=255)

                                                  gCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc                             >  2:607708/1‑93 (MQ=255)

                                                  gCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc                             >  1:5537918/1‑93 (MQ=255)

                                                  gCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc                             >  1:5691482/1‑93 (MQ=255)

                                                    tGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGc                           <  1:2637188/93‑1 (MQ=255)

                                                    tGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGc                           >  1:6377530/1‑93 (MQ=255)

                                                      ttATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGccc                         >  1:3016353/1‑93 (MQ=255)

                                                          ccTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCa                     >  2:4053867/1‑93 (MQ=255)

                                                           cTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGATTCGGTGCATTCAGCGACTTGTCGCCCGTAAg                    >  1:2024120/1‑93 (MQ=255)

                                                             gACCACCATGTAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGtt                  <  1:5350826/93‑1 (MQ=255)

                                                             gACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGtt                  <  1:1824317/93‑1 (MQ=255)

                                                             gACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGtt                  <  1:1022471/93‑1 (MQ=255)

                                                             gACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGtt                  <  1:730333/93‑1 (MQ=255)

                                                                caccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGc               >  1:2149113/1‑93 (MQ=255)

                                                                caccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGc               >  2:340497/1‑93 (MQ=255)

                                                                caccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGc               <  2:5904738/93‑1 (MQ=255)

                                                                 accaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGcc              >  2:6682420/1‑93 (MQ=255)

                                                                 accaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGcc              >  2:728656/1‑93 (MQ=255)

                                                                  ccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCg             <  1:5826954/93‑1 (MQ=255)

                                                                  ccaTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCg             >  2:4791220/1‑93 (MQ=255)

                                                                   caTATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGt            >  1:6626310/1‑93 (MQ=255)

                                                                   caTAAAAGCCGAAAAACGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTTAATTCAGCGACTTGTCGCGCGTCAGTTTGCCGt            <  1:3400652/93‑1 (MQ=255)

                                                                    atatAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAACGACTTGTCGCCCGTCAGTTTGCCGTa           >  2:276638/1‑93 (MQ=255)

                                                                     tataAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAc          <  1:3640193/93‑1 (MQ=255)

                                                                      ataAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAcc         <  1:5950391/93‑1 (MQ=255)

                                                                      ataAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAcc         >  2:2252138/1‑93 (MQ=255)

                                                                      ataAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAcc         <  1:4959336/93‑1 (MQ=255)

                                                                       taAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCa        >  2:2498941/1‑93 (MQ=255)

                                                                         aGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACg      >  1:1777571/1‑93 (MQ=255)

                                                                           ccGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGcc    <  2:3269351/93‑1 (MQ=255)

                                                                           ccGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGcc    <  1:6254041/93‑1 (MQ=255)

                                                                           ccGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGcc    >  2:1446109/1‑93 (MQ=255)

                                                                            cGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCa   <  1:3842010/93‑1 (MQ=255)

                                                                             gagaATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGGATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAc  <  2:1394781/93‑1 (MQ=255)

TTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCCG‑GCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCAC  >  NC_000913/4013092‑4013260

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