BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,573,190	Δ1 bp	coding (1487/1542 nt)	yjiU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,189	0	G	.	99.7%	1487.4 / ‑3.2	328	coding (1486/1542 nt)	yjiU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base . (168/158); total (170/158)
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 353 total aligned reads displayed.

TTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAACCC  >  NC_000913/4573097‑4573277

ttatCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAg                                                                                          <  2:583999/93‑1 (MQ=255)

ttatCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAg                                                                                          >  2:3647281/1‑93 (MQ=255)

 tatCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAga                                                                                         <  2:3724204/93‑2 (MQ=255)

   tCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAgagc                                                                                       >  1:5787259/1‑90 (MQ=255)

   tCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAgagc                                                                                       <  1:7029128/93‑4 (MQ=255)

     aaTACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCat                                                                                    >  2:4406171/1‑93 (MQ=255)

      aTACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCata                                                                                   <  2:2110510/93‑1 (MQ=255)

      aTACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCata                                                                                   <  1:272754/93‑1 (MQ=255)

      aTACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCata                                                                                   <  1:6539148/93‑1 (MQ=255)

       tACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAGACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCatat                                                                                  >  1:5902687/1‑93 (MQ=255)

       tACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCatat                                                                                  >  1:4968303/1‑93 (MQ=255)

       tACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCatat                                                                                  >  2:756049/1‑93 (MQ=255)

        aCATGGCTGAAAGAAGCATGGGCCGAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATg                                                                                 <  1:1672326/93‑1 (MQ=255)

        aCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATg                                                                                 <  2:6106643/93‑1 (MQ=255)

        aCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATg                                                                                 <  1:3751626/93‑1 (MQ=255)

         cATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGa                                                                                <  1:4796888/93‑1 (MQ=255)

          aTGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAGCCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGAc                                                                               <  1:2115681/93‑1 (MQ=255)

          aTGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGAc                                                                               <  2:648562/93‑1 (MQ=255)

          aTGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGAc                                                                               <  1:3379064/93‑1 (MQ=255)

           tGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACt                                                                              <  2:4675370/93‑1 (MQ=255)

           tGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACt                                                                              >  1:6908385/1‑93 (MQ=255)

           tGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACt                                                                              <  1:7162614/93‑1 (MQ=255)

           tGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACt                                                                              >  1:1296565/1‑93 (MQ=255)

            ggCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACtt                                                                             <  1:2386193/93‑1 (MQ=255)

             gCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACttt                                                                            >  1:3929872/1‑93 (MQ=255)

             gCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACttt                                                                            >  1:4587389/1‑93 (MQ=255)

              cTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACtttt                                                                           <  2:4198777/93‑1 (MQ=255)

               tGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACttttt                                                                          >  1:2214170/1‑93 (MQ=255)

                gAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTc                                                                         >  2:5971395/1‑93 (MQ=255)

                  aagaagCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCtt                                                                       >  1:1001789/1‑93 (MQ=255)

                  aagaagCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCtt                                                                       <  1:7715868/93‑1 (MQ=255)

                   agaagCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCttt                                                                      <  1:6808591/93‑1 (MQ=255)

                    gaagCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTa                                                                     <  1:7405849/93‑1 (MQ=255)

                      agCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAAt                                                                   <  1:5139258/93‑1 (MQ=255)

                       gCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATc                                                                  >  1:5560179/1‑93 (MQ=255)

                       gCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATc                                                                  >  1:4097576/1‑93 (MQ=255)

                       gCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATc                                                                  >  1:1487650/1‑93 (MQ=255)

                        cATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCa                                                                 >  2:3973617/1‑93 (MQ=255)

                        cATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCa                                                                 >  2:6600002/1‑93 (MQ=255)

                           gggCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACgg                                                              >  2:2040884/1‑93 (MQ=255)

                            ggCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGa                                                             <  1:7896034/93‑1 (MQ=255)

                             gCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGAt                                                            <  1:204574/93‑1 (MQ=255)

                                cAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGcc                                                         >  2:3295328/1‑93 (MQ=255)

                                 aaCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCt                                                        <  2:2169924/93‑1 (MQ=255)

                                  aCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTg                                                       >  2:6254702/1‑93 (MQ=255)

                                  aCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTg                                                       >  2:1708250/1‑93 (MQ=255)

                                    ttACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGa                                                     >  1:1095000/1‑93 (MQ=255)

                                    ttACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGa                                                     <  2:2772418/93‑1 (MQ=255)

                                      tccGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAtt                                                   <  2:3927383/92‑1 (MQ=255)

                                      aCCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAtt                                                   <  2:744794/93‑1 (MQ=255)

                                         ggACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGc                                                <  1:4300040/93‑1 (MQ=255)

                                          gACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGcc                                               >  2:3178812/1‑93 (MQ=255)

                                           aCCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGccc                                              <  2:4725718/93‑1 (MQ=255)

                                              gCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGt                                           >  2:6970338/1‑93 (MQ=255)

                                              gCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGt                                           >  2:2677802/1‑93 (MQ=255)

                                               cTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTa                                          <  2:6423103/93‑1 (MQ=255)

                                                ttGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTaa                                         >  2:1217331/1‑93 (MQ=255)

                                                ttGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTaa                                         >  2:7930067/1‑93 (MQ=255)

                                                  gCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAcc                                       >  2:7343285/1‑93 (MQ=255)

                                                   cAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCa                                      <  2:4911371/93‑1 (MQ=255)

                                                    aGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCaa                                     >  1:1813385/1‑93 (MQ=255)

                                                     gTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAAt                                    <  2:2547772/93‑1 (MQ=255)

                                                        cTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaaa                                 <  2:1487521/93‑1 (MQ=255)

                                                        cTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaaa                                 >  2:3551492/1‑93 (MQ=255)

                                                         tGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaaaa                                >  1:2351180/1‑93 (MQ=255)

                                                             aCCCGATGGCAACCGTTTTTATTTACCGAAA‑GTGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCt                            >  1:5712971/1‑93 (MQ=255)

                                                              cccGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTg                           <  1:3716106/93‑1 (MQ=255)

                                                               ccGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTgg                          <  2:533931/93‑1 (MQ=255)

                                                                cGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGa                         >  1:3593751/1‑93 (MQ=255)

                                                                cGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGa                         <  1:87582/93‑1 (MQ=255)

                                                                  aTGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATa                       <  2:1776992/93‑1 (MQ=255)

                                                                  aTGCCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATa                       <  1:2973197/93‑1 (MQ=255)

                                                                   tGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATAc                      <  2:5699282/93‑1 (MQ=255)

                                                                   tGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATAc                      >  1:3100228/1‑93 (MQ=255)

                                                                   tGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATAc                      >  1:7594562/1‑93 (MQ=255)

                                                                   tGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATAc                      <  2:898499/93‑1 (MQ=255)

                                                                       aaCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCt                  >  1:504744/1‑93 (MQ=255)

                                                                         ccGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGTCCGGTAACCAATAAAATCCTGGATACGGCttt                >  2:3086427/1‑93 (MQ=255)

                                                                         ccGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCttt                >  1:6662881/1‑93 (MQ=255)

                                                                         ccGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCttt                >  1:67960/1‑93 (MQ=255)

                                                                          cGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCtttt               >  1:4408225/1‑93 (MQ=255)

                                                                          cGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCtttt               <  1:4736386/93‑1 (MQ=255)

                                                                          cGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCtttt               <  1:1637094/93‑1 (MQ=255)

                                                                           gTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTa              >  2:4541390/1‑93 (MQ=255)

                                                                           gTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTa              >  2:3007592/1‑93 (MQ=255)

                                                                           gTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTa              >  2:7779858/1‑93 (MQ=255)

                                                                           gTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTa              >  1:7474368/1‑93 (MQ=255)

                                                                            tttttATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTaa             <  2:5583134/93‑1 (MQ=255)

                                                                            tttttATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTaa             <  2:8152/93‑1 (MQ=255)

                                                                             ttttATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTaaa            >  1:5393666/1‑93 (MQ=255)

                                                                             ttttATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTaaa            <  1:544128/93‑1 (MQ=255)

                                                                             ttttATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTaaa            <  1:6217748/93‑1 (MQ=255)

                                                                                 atttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCt        <  1:7849584/93‑1 (MQ=255)

                                                                                 atttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCt        <  1:3487109/93‑1 (MQ=255)

                                                                                  tttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCtt       >  2:4116829/1‑93 (MQ=255)

                                                                                   ttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTa      >  2:6801878/1‑93 (MQ=255)

                                                                                   ttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTa      >  2:5154804/1‑93 (MQ=255)

                                                                                    taCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTaa     <  1:2510717/93‑1 (MQ=255)

                                                                                     aCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAAc    <  2:561939/93‑1 (MQ=255)

                                                                                       cGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAAccc  <  2:2441809/93‑1 (MQ=255)

TTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGCTTTTAAAGGCTTAACCC  >  NC_000913/4573097‑4573277

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