BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,945,882	A→G	K58K (AAA→AAG)	yifE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,945,882	0	A	G	82.5%	761.7 / 144.0	320	K58K (AAA→AAG)	yifE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (34/22); new base G (127/137); total (161/159)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 326 total aligned reads displayed.

ACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCC  >  NC_000913/3945790‑3945971

 ccATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAAc                                                                                          >  1:5740574/1‑93 (MQ=255)

  cATCAAAGAGGCACAATTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGct                                                                                         <  2:2236829/93‑1 (MQ=255)

  cATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGct                                                                                         <  2:5940942/93‑1 (MQ=255)

  cATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAAct                                                                                         <  2:1679464/93‑1 (MQ=255)

  cATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAAct                                                                                         >  1:2170198/1‑93 (MQ=255)

   aTCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAActc                                                                                        <  2:5078781/93‑1 (MQ=255)

     cAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCtt                                                                                      <  1:4481659/93‑1 (MQ=255)

      aaaGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTc                                                                                     <  2:3865563/93‑1 (MQ=255)

      aaaGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTc                                                                                     <  1:364429/93‑1 (MQ=255)

      aaaGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTc                                                                                     <  2:4592287/93‑1 (MQ=255)

       aaGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCg                                                                                    >  1:3311862/1‑93 (MQ=255)

        agagGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGt                                                                                   >  1:2823797/1‑93 (MQ=255)

        agagGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGt                                                                                   >  1:1307575/1‑93 (MQ=255)

         gagGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTa                                                                                  >  1:2566561/1‑93 (MQ=255)

          agGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAg                                                                                 >  1:1629180/1‑93 (MQ=255)

           ggCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGc                                                                                <  2:1912238/93‑1 (MQ=255)

           ggCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGc                                                                                <  1:5537554/93‑1 (MQ=255)

              acaACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGt                                                                             <  1:303936/93‑1 (MQ=255)

               caACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTa                                                                            <  1:4653195/93‑1 (MQ=255)

               caACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTa                                                                            <  1:5421947/93‑1 (MQ=255)

                atcTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTAt                                                                           <  2:5606450/91‑1 (MQ=255)

                 aCTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATg                                                                          <  2:6336140/93‑1 (MQ=255)

                  ctgctgGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGc                                                                         <  1:2969642/93‑1 (MQ=255)

                  cTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGc                                                                         <  2:4778168/93‑1 (MQ=255)

                     gttGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGt                                                                      <  1:244859/92‑1 (MQ=255)

                     cTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGt                                                                      >  2:523944/1‑93 (MQ=255)

                      ttGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTg                                                                     >  2:5721862/1‑93 (MQ=255)

                      ttGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTg                                                                     <  2:1206240/93‑1 (MQ=255)

                       tGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTgg                                                                    >  1:140856/1‑93 (MQ=255)

                       tGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTgg                                                                    >  1:2473961/1‑93 (MQ=255)

                       tGAGCGTCATGGTGATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCTGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTgg                                                                    <  2:4666313/93‑1 (MQ=255)

                        gAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGc                                                                   <  1:1082998/93‑1 (MQ=255)

                        gAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGc                                                                   <  1:5680371/93‑1 (MQ=255)

                          gCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGa                                                                 <  2:3235679/93‑1 (MQ=255)

                            gTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAAc                                                               >  1:479116/1‑93 (MQ=255)

                             tCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACg                                                              >  1:1838324/1‑93 (MQ=255)

                              cATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGt                                                             <  2:4021007/93‑1 (MQ=255)

                               aTGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTg                                                            >  2:2499355/1‑93 (MQ=255)

                               aTGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTg                                                            >  1:63153/1‑93 (MQ=255)

                                  gTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGc                                                         >  2:1229474/1‑93 (MQ=255)

                                    tATGCTTTCAATGAGTTGGATCTTGGCAACCGGGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCa                                                       <  2:4427057/93‑1 (MQ=255)

                                     aTGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAg                                                      <  2:879833/93‑1 (MQ=255)

                                     aTGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAg                                                      <  2:4879650/93‑1 (MQ=255)

                                     aTGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAg                                                      >  1:3944128/1‑93 (MQ=255)

                                        ccTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGa                                                   >  2:3357286/1‑93 (MQ=255)

                                           tCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAg                                                >  1:3129850/1‑93 (MQ=255)

                                            cAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGa                                               >  2:4959747/1‑93 (MQ=255)

                                            cAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGa                                               <  2:392935/93‑1 (MQ=255)

                                             aaTGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGaa                                              <  1:3248340/93‑1 (MQ=255)

                                             aaTGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGaa                                              <  2:3380863/93‑1 (MQ=255)

                                              aTGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAg                                             >  2:4093980/1‑93 (MQ=255)

                                               tGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGc                                            <  2:2623653/93‑1 (MQ=255)

                                                gAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCa                                           >  1:3400626/1‑93 (MQ=255)

                                                gAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCa                                           <  2:2367553/93‑1 (MQ=255)

                                                gAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCa                                           <  2:371540/93‑1 (MQ=255)

                                                  gTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGa                                         >  2:5288508/1‑93 (MQ=255)

                                                    tGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAc                                       >  2:6070310/1‑93 (MQ=255)

                                                    tGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAc                                       <  1:4190288/93‑1 (MQ=255)

                                                     ggATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcg                                      >  1:199147/1‑93 (MQ=255)

                                                      gATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgc                                     <  2:380277/93‑1 (MQ=255)

                                                         cTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgtg                                  >  1:6357145/1‑93 (MQ=255)

                                                         cTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgtg                                  <  1:1708159/93‑1 (MQ=255)

                                                          ttGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgtgt                                 <  2:2040956/93‑1 (MQ=255)

                                                            ggCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTgg                               >  1:4721747/1‑93 (MQ=255)

                                                             gCAAACGCGAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGt                              <  2:1782211/93‑1 (MQ=255)

                                                               aatcgcgAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTcc                            <  2:2338803/90‑1 (MQ=255)

                                                               aaaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTcc                            >  1:4904567/1‑93 (MQ=255)

                                                                  cgcgAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAg                         <  1:1282420/93‑1 (MQ=255)

                                                                  cgcgAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAg                         >  1:6580784/1‑93 (MQ=255)

                                                                   gcgAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGt                        <  1:4372827/93‑1 (MQ=255)

                                                                   gcgAGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGt                        >  1:968411/1‑93 (MQ=255)

                                                                   gcgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGt                        <  1:2421164/93‑1 (MQ=255)

                                                                    cgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGta                       <  2:1377289/93‑1 (MQ=255)

                                                                      aGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                     <  1:1863672/93‑1 (MQ=255)

                                                                      aGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtata                     <  1:4573381/93‑1 (MQ=255)

                                                                      aGCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGGGACAGAACCAGAACGCGTGTGGTCCAAGtata                     >  2:6006040/1‑93 (MQ=255)

                                                                       gCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTatat                    <  2:1508276/93‑1 (MQ=255)

                                                                       gCCGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTatat                    <  1:5125937/93‑1 (MQ=255)

                                                                        ccGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATg                   <  2:3078541/93‑1 (MQ=255)

                                                                        ccGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATg                   >  1:3814739/1‑93 (MQ=255)

                                                                         cGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGa                  <  1:5894236/93‑1 (MQ=255)

                                                                         cGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGa                  <  1:2318790/93‑1 (MQ=255)

                                                                         cGGTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGa                  <  2:4535862/93‑1 (MQ=255)

                                                                         cGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGa                  >  2:6655206/1‑93 (MQ=255)

                                                                           gTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcg                >  2:1283241/1‑93 (MQ=255)

                                                                           gTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcg                <  2:714219/93‑1 (MQ=255)

                                                                           gTTACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcg                <  1:5308310/93‑1 (MQ=255)

                                                                             tACCGAGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcgcg              >  2:6504831/1‑93 (MQ=255)

                                                                               ccGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTa            >  1:3840095/1‑93 (MQ=255)

                                                                                  aGGAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTa         <  1:6122396/93‑1 (MQ=255)

                                                                                   ggAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTaa        <  2:4284293/93‑1 (MQ=255)

                                                                                    gAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAg       <  1:2752394/93‑1 (MQ=255)

                                                                                    gAAGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAg       <  2:5483158/93‑1 (MQ=255)

                                                                                     aaGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGc      <  1:2374398/93‑1 (MQ=255)

                                                                                     aaGAGAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGc      <  2:587245/93‑1 (MQ=255)

                                                                                      agagAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCg     >  2:3202794/1‑93 (MQ=255)

                                                                                      agagAAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCg     >  2:6187816/1‑93 (MQ=255)

                                                                                       gagaAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTGTATGACGCGTATTAAGCGt    <  1:6484249/93‑1 (MQ=255)

                                                                                         gaAGCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTcc  >  1:5011748/1‑93 (MQ=255)

                                                                                         gaAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTcc  >  2:289940/1‑93 (MQ=255)

ACCATCAAAGAGGCACAACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCC  >  NC_000913/3945790‑3945971

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