BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,589,600	G→A	G134G (GGC→GGT)	ydeT ←	fimbrial usher domain‑containing protein YdeT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,589,600	0	G	A	99.7%	1214.7 / NA	285	G134G (GGC→GGT)	ydeT	fimbrial usher domain‑containing protein YdeT
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (130/154); total (130/155)
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 293 total aligned reads displayed.

GACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATGCCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACATTAA  >  NC_000913/1589502‑1589697

 aCGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATaccg                                                                                                >  2:2591926/1‑101 (MQ=255)

  cGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATaccgc                                                                                               <  2:2591567/101‑1 (MQ=255)

  cGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATaccgc                                                                                               >  1:2591929/1‑101 (MQ=255)

    aTGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCt                                                                                             <  1:2591415/101‑1 (MQ=255)

    aTGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCt                                                                                             >  2:2591934/1‑101 (MQ=255)

     tGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCtt                                                                                            <  1:2591465/101‑1 (MQ=255)

      gTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCttt                                                                                           >  1:2591937/1‑101 (MQ=255)

       tATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTc                                                                                          <  2:2591469/101‑1 (MQ=255)

       tATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTc                                                                                          <  1:2591357/101‑1 (MQ=255)

        aTTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCa                                                                                         <  2:2591409/101‑1 (MQ=255)

        aTTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCa                                                                                         <  2:2591329/101‑1 (MQ=255)

        aTTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCa                                                                                         >  1:2591940/1‑101 (MQ=255)

            ccTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAAt                                                                                      <  1:2591528/100‑1 (MQ=255)

              tCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGt                                                                                   <  1:2591202/101‑1 (MQ=255)

              tCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGt                                                                                   >  1:2591944/1‑101 (MQ=255)

              tCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGt                                                                                   <  1:2591331/101‑1 (MQ=255)

                cGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGttt                                                                                 >  2:2591948/1‑101 (MQ=255)

                cGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGttt                                                                                 <  2:2591353/101‑1 (MQ=255)

                 gTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTg                                                                                <  1:2591429/101‑1 (MQ=255)

                 gTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTg                                                                                <  1:2591248/101‑1 (MQ=255)

                 gTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTACTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTg                                                                                <  1:2591407/101‑1 (MQ=255)

                   gggTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGAc                                                                              <  1:2591511/101‑1 (MQ=255)

                     gtgtTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACAt                                                                            <  2:2591601/101‑1 (MQ=255)

                      tgtTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATa                                                                           >  2:2591951/1‑101 (MQ=255)

                        ttACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACt                                                                         <  2:2591698/101‑1 (MQ=255)

                         tACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTg                                                                        <  2:2591519/101‑1 (MQ=255)

                          aCCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGt                                                                       >  2:2591957/1‑101 (MQ=255)

                          aCCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGt                                                                       <  1:2591619/101‑1 (MQ=255)

                          aCCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGt                                                                       <  1:2591745/101‑1 (MQ=255)

                            cGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTaa                                                                     <  2:2591577/101‑1 (MQ=255)

                             gACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAAc                                                                    >  2:2591960/1‑101 (MQ=255)

                               ccTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTg                                                                  >  1:2591962/1‑101 (MQ=255)

                               ccTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTg                                                                  <  2:2591589/101‑1 (MQ=255)

                                  gAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCg                                                               >  1:2591966/1‑101 (MQ=255)

                                    aCGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGtt                                                             >  2:2591969/1‑101 (MQ=255)

                                     cGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGttt                                                            >  2:2591972/1‑101 (MQ=255)

                                      gCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTg                                                           >  1:2591974/1‑101 (MQ=255)

                                        tataATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGt                                                              <  1:2591503/96‑1 (MQ=255)

                                        tataATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGaa                                                         >  1:2591978/1‑101 (MQ=255)

                                         ataATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAAt                                                        >  2:2591979/1‑101 (MQ=255)

                                          taATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAAtt                                                       <  2:2591368/101‑1 (MQ=255)

                                          taATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAAtt                                                       <  1:2591145/101‑1 (MQ=255)

                                          taATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAAtt                                                       <  2:2591547/101‑1 (MQ=255)

                                            aTTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTAc                                                     <  2:2591419/101‑1 (MQ=255)

                                             ttCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACg                                                    <  2:2591305/101‑1 (MQ=255)

                                             ttCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACg                                                    >  1:2591984/1‑101 (MQ=255)

                                              tCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGa                                                   >  2:2591985/1‑101 (MQ=255)

                                              tCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGa                                                   <  2:2591458/101‑1 (MQ=255)

                                                aGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGaaa                                                 <  1:2591717/101‑1 (MQ=255)

                                                 ggTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAAt                                                >  1:2591989/1‑101 (MQ=255)

                                                    tattatCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGcc                                             <  1:2591557/101‑1 (MQ=255)

                                                      ttatCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCga                                           >  1:2591992/1‑101 (MQ=255)

                                                       tatCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCgac                                          <  2:2591629/101‑1 (MQ=255)

                                                       tatCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCgac                                          >  2:2591993/1‑101 (MQ=255)

                                                        atCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCgact                                         >  2:2591995/1‑101 (MQ=255)

                                                        atCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCgact                                         <  1:2591370/101‑1 (MQ=255)

                                                        atCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCgact                                         <  2:2591438/101‑1 (MQ=255)

                                                         tCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCgactg                                        <  2:2591494/101‑1 (MQ=255)

                                                         tCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCgactg                                        <  1:2591291/101‑1 (MQ=255)

                                                           cGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCg                                                                                    <  2:2591232/55‑1 (MQ=255)

                                                              cagcagAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGActgt                                   <  1:2591664/101‑1 (MQ=255)

                                                              cagcagAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGActgt                                   >  1:2592001/1‑101 (MQ=255)

                                                              cagcagAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGActgt                                   >  2:2592000/1‑101 (MQ=255)

                                                                gcagAGTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTg                                                                        >  1:2592005/1‑62 (MQ=255)

                                                                  agagTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGActgtctgt                               <  1:2591479/101‑1 (MQ=255)

                                                                  agagTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGActgtctgt                               >  2:2592008/1‑101 (MQ=255)

                                                                   gagTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTa                              <  2:2591833/101‑1 (MQ=255)

                                                                    agTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTAc                             >  2:2592011/1‑101 (MQ=255)

                                                                     gTGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACg                            <  2:2591615/101‑1 (MQ=255)

                                                                      tGCCATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGc                           >  1:2592013/1‑101 (MQ=255)

                                                                        ccATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCAt                         >  2:2592017/1‑101 (MQ=255)

                                                                        ccATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCAt                         <  1:2591645/101‑1 (MQ=255)

                                                                         cATAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATc                        >  1:2592018/1‑101 (MQ=255)

                                                                           tAAACCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCa                      <  1:2591439/101‑1 (MQ=255)

                                                                            aaaCCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCa                          <  1:2591487/96‑1 (MQ=255)

                                                                            aaaCCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCaa                     >  1:2592022/1‑101 (MQ=255)

                                                                            aaaCCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCaa                     <  1:2591318/101‑1 (MQ=255)

                                                                            aaaCCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCaa                     >  2:2592025/1‑101 (MQ=255)

                                                                             aaCCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAAt                    <  2:2591525/101‑1 (MQ=255)

                                                                             aaCCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAAt                    >  2:2592029/1‑101 (MQ=255)

                                                                              aCCCCCGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATg                   >  1:2592031/1‑101 (MQ=255)

                                                                                ccccGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGAc                 >  2:2592033/1‑101 (MQ=255)

                                                                                  ccGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACAt                                                                            >  2:2592037/1‑40 (MQ=255)

                                                                                  ccGATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTg               <  2:2591648/101‑1 (MQ=255)

                                                                                    gATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGaa             >  2:2592041/1‑101 (MQ=255)

                                                                                    gATAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGaa             <  1:2591581/101‑1 (MQ=255)

                                                                                     aTAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAg            >  1:2592043/1‑101 (MQ=255)

                                                                                      tAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAgg           >  2:2592048/1‑101 (MQ=255)

                                                                                      tAGATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAgg           >  1:2592046/1‑101 (MQ=255)

                                                                                        gATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGa         <  1:2591565/101‑1 (MQ=255)

                                                                                        gATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGa         >  1:2592052/1‑101 (MQ=255)

                                                                                        gATTGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGa         >  2:2592054/1‑101 (MQ=255)

                                                                                          ttGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAAc       >  2:2592057/1‑101 (MQ=255)

                                                                                          ttGGTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAAc       <  2:2591299/101‑1 (MQ=255)

                                                                                            ggTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACAt     <  2:2591400/101‑1 (MQ=255)

                                                                                             gTCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACAtt    >  2:2592061/1‑101 (MQ=255)

                                                                                              tCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGaa        >  1:2592065/1‑96 (MQ=255)

                                                                                              tCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACATTa   <  1:2591460/101‑1 (MQ=255)

                                                                                              tCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACATTa   <  1:2591548/101‑1 (MQ=255)

                                                                                              tCATACCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACATTa   <  1:2591913/101‑1 (MQ=255)

GACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATGCCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACATTAA  >  NC_000913/1589502‑1589697

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`