BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	271,281	G→A	G76G (GGG→GGA)	b0257 →	putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	271,281	0	G	A	99.3%	975.6 / ‑4.6	278	G76G (GGG→GGA)	b0257	putative transposase
Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (122/154); total (124/154)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 297 total aligned reads displayed.

TGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGT  >  NC_000913/271191‑271373

tGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACt                                                                                            >  1:613087/1‑93 (MQ=255)

 gCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTg                                                                                           >  2:4009475/1‑93 (MQ=255)

 gCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTg                                                                                           >  2:360701/1‑93 (MQ=255)

    ggAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTc                                                                                        >  2:1721486/1‑93 (MQ=255)

     gaccAAGAAGCATCGCAACAAGGGAAACCCAGAGAGGATACCAGATGGGGCGCTGCCTTGCTGGCAGCCTCATGAAAGAGCCGGGACTGGTCa                                                                                       <  2:2300171/90‑1 (MQ=255)

        cAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGtt                                                                                    <  1:4295264/93‑1 (MQ=255)

         aagaagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTAGGACCGGGCAGTTg                                                                                   >  1:2091697/1‑93 (MQ=17)

          agaagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGATCTTGGTCTGGTCAGTTGc                                                                                  >  1:5068029/1‑93 (MQ=255)

          agaagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGc                                                                                  >  2:7221835/1‑93 (MQ=255)

          agaagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGc                                                                                  >  1:526226/1‑93 (MQ=255)

            aagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCca                                                                                <  2:482028/93‑1 (MQ=255)

            aagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCca                                                                                <  2:5185594/93‑1 (MQ=255)

            aagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCca                                                                                <  1:5685988/93‑1 (MQ=255)

              gCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCcagc                                                                              <  2:567919/93‑1 (MQ=255)

              gCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCcagc                                                                              >  2:1833393/1‑93 (MQ=255)

               cATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCcagca                                                                             <  1:199826/93‑1 (MQ=255)

               cATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCcagca                                                                             <  1:4767973/93‑1 (MQ=255)

                aTCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCcagcag                                                                            >  1:1145008/1‑93 (MQ=255)

                aTCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCcagcag                                                                            >  2:4540343/1‑93 (MQ=255)

                 tCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCagcagc                                                                           <  1:1007447/93‑1 (MQ=255)

                   gCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGTCg                                                                         >  1:2788545/1‑93 (MQ=255)

                   gCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCg                                                                         >  1:7000548/1‑93 (MQ=255)

                   cccACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCg                                                                         >  2:1168137/2‑93 (MQ=255)

                     cacaATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGAc                                                                       <  1:3259444/93‑1 (MQ=255)

                        aaTGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTcc                                                                    >  1:6683938/1‑92 (MQ=255)

                        aaTGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCa                                                                    >  1:2838156/1‑93 (MQ=255)

                              aaCCCAGAGAGGATACCCGATGGGGCGCTGGCTTGCTGGCCGACTCCTGAAAGAGCTGGGACAGGGCAGTTGCCCGCACCCGACTCCCCggat                                                              >  2:6028870/1‑91 (MQ=255)

                              aaCCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGta                                                              <  1:4524940/93‑1 (MQ=255)

                               aCCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGtat                                                             >  1:5938574/1‑93 (MQ=255)

                               aCCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGtat                                                             >  2:2023756/1‑93 (MQ=255)

                                cccAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGtata                                                            <  2:4952883/93‑1 (MQ=255)

                                 ccAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATaa                                                           <  1:2597522/93‑1 (MQ=255)

                                   agagagTATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGc                                                         <  2:4311766/93‑1 (MQ=255)

                                   agagagGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGc                                                         <  1:3163198/93‑1 (MQ=255)

                                    gagagGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCg                                                        >  1:354477/1‑93 (MQ=255)

                                    gagagGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCg                                                        >  2:1327945/1‑93 (MQ=255)

                                    gagagGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCg                                                        >  2:677235/1‑93 (MQ=255)

                                     agagGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGt                                                       <  2:2715845/93‑1 (MQ=255)

                                     agagGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGt                                                       <  1:4422147/93‑1 (MQ=255)

                                         gATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCg                                                   >  1:1523395/1‑93 (MQ=255)

                                          aTACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCgg                                                  >  1:1897026/1‑93 (MQ=255)

                                           tACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGt                                                 >  1:1342814/1‑93 (MQ=255)

                                            aCCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTc                                                <  2:7381588/93‑1 (MQ=255)

                                             ccAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCa                                               >  1:3406808/1‑93 (MQ=255)

                                             ccAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCa                                               <  1:1056513/93‑1 (MQ=255)

                                             ccAGATGGGGCGCTGGCTTCCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCa                                               <  1:7693083/93‑1 (MQ=255)

                                              cAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCAt                                              <  2:2123734/93‑1 (MQ=255)

                                                gATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCCGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGa                                            >  1:5350826/1‑93 (MQ=255)

                                                gATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGa                                            <  1:7297307/93‑1 (MQ=255)

                                                gATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGa                                            >  1:5866362/1‑93 (MQ=255)

                                                  tGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGc                                          <  2:6697790/93‑1 (MQ=255)

                                                   ggggCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCa                                         <  1:4201781/93‑1 (MQ=255)

                                                    gggCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCAc                                        >  1:7215724/1‑93 (MQ=255)

                                                      gcgcTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGt                                      <  2:1666556/93‑1 (MQ=255)

                                                       ggcaggCTTGCTGCCAGACACATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGCCTCACCGGTATAAGCGTGGCGGTCATGAGCACGtt                                     <  1:3019659/89‑1 (MQ=255)

                                                       cgcTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGtt                                     <  1:2288565/93‑1 (MQ=255)

                                                        gcTGGCTTGCTGTCAGACTCATGAAAGAGCTGAGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTg                                    <  1:1423981/93‑1 (MQ=255)

                                                        gcTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTg                                    >  2:3746407/1‑93 (MQ=255)

                                                           ggCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTa                                 <  1:4612989/93‑1 (MQ=255)

                                                            gCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCCGTCATGAGCACGTTGCTAt                                >  1:5427558/1‑93 (MQ=255)

                                                             cTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATc                               <  1:4443103/93‑1 (MQ=255)

                                                              ttGCTGGCAGTCTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATcc                              <  1:7396883/93‑1 (MQ=255)

                                                              ttGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATcc                              >  2:3019016/1‑93 (MQ=255)

                                                               tGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATccc                             <  2:925806/93‑1 (MQ=255)

                                                               tGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATccc                             <  2:750364/93‑1 (MQ=255)

                                                               tGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATccc                             <  1:86710/93‑1 (MQ=255)

                                                                   ggCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAat                         >  2:2404977/1‑93 (MQ=255)

                                                                   ggCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAat                         <  1:4122273/93‑1 (MQ=255)

                                                                     cAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatca                       <  2:7336452/93‑1 (MQ=255)

                                                                      aggcTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCAGGTGCACGTTGCGATCCCGAatcat                      >  2:3208508/4‑93 (MQ=18)

                                                                      aGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatcat                      <  2:6153059/93‑1 (MQ=255)

                                                                       gACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatcatc                     >  1:6410594/1‑93 (MQ=255)

                                                                       gACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatcatc                     >  1:3619460/1‑93 (MQ=255)

                                                                       gACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatcatc                     <  1:6825489/93‑1 (MQ=255)

                                                                        gcTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCt                    <  2:6847428/92‑1 (MQ=255)

                                                                        aCTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCt                    <  2:1225631/93‑1 (MQ=255)

                                                                        aCTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCt                    <  2:6987525/93‑1 (MQ=255)

                                                                        aCTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCt                    <  2:4622461/93‑1 (MQ=255)

                                                                         cTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCtt                   <  1:2672213/93‑1 (MQ=255)

                                                                           cATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGa                 >  2:390952/1‑93 (MQ=255)

                                                                           cATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGa                 <  2:5481064/93‑1 (MQ=255)

                                                                           cATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGa                 <  2:5129456/93‑1 (MQ=255)

                                                                            aTGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAg                >  2:4972971/1‑93 (MQ=255)

                                                                             tGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGc               <  1:5398334/93‑1 (MQ=255)

                                                                             tGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGc               <  2:593031/93‑1 (MQ=255)

                                                                             tGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGc               >  2:4491725/1‑93 (MQ=255)

                                                                             tGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTAGCCCGAATCATCTTGAGc               >  1:7047311/1‑93 (MQ=255)

                                                                               aaaGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGa             <  2:1534693/93‑1 (MQ=255)

                                                                                aaGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGAc            <  1:2903600/93‑1 (MQ=255)

                                                                                aaGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGAc            <  1:6089573/93‑1 (MQ=255)

                                                                                 agagCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACa           <  2:2568682/93‑1 (MQ=255)

                                                                                   agCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGt         >  1:2404724/1‑93 (MQ=255)

                                                                                   agCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGt         <  2:7699099/93‑1 (MQ=255)

                                                                                    gCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATTCCGAATCATCTTGAGCGACAGtt        >  2:782237/1‑93 (MQ=255)

                                                                                    gCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGtt        <  2:6219033/93‑1 (MQ=255)

                                                                                    gCTGGGACTGGTCAGTTGACAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGtt        <  1:3751520/93‑1 (MQ=255)

                                                                                     cTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTc       <  1:5111461/93‑1 (MQ=255)

                                                                                     cTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTc       >  2:159415/1‑93 (MQ=255)

                                                                                      tGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTcc      >  2:2689485/1‑92 (MQ=255)

                                                                                      tGGGACTGGTAAGTTGCCAGCAGCCGACTCACCGGAATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATTAACTTGAGTGACAGTTCg      >  2:1401821/1‑93 (MQ=18)

TGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGT  >  NC_000913/271191‑271373

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