BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,789,635	T→C	Y137C (TAT→TGT)	b1707 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,789,635	0	T	C	99.1%	1108.6 / ‑2.0	344	Y137C (TAT→TGT)	b1707	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (170/171); total (172/172)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 363 total aligned reads displayed.

CTTTTATCCAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCAATAAACGGA  >  NC_000913/1789545‑1789724

   ttATCCAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACt                                                                                      <  1:2505963/93‑1 (MQ=255)

   ttATCCAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACt                                                                                      <  1:5089885/93‑1 (MQ=255)

      tCCAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGc                                                                                   >  1:7042039/1‑93 (MQ=255)

      tCCAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGc                                                                                   >  2:6919351/1‑93 (MQ=255)

        cAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGa                                                                                 >  1:9171019/1‑93 (MQ=255)

          gCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATa                                                                               <  2:6149702/93‑1 (MQ=255)

           cATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATaa                                                                              <  2:5766588/93‑1 (MQ=255)

            aTGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAg                                                                             <  1:10036831/93‑1 (MQ=255)

            aTGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAg                                                                             >  1:6236861/1‑93 (MQ=255)

             tGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAgg                                                                            <  2:9732162/93‑1 (MQ=255)

             tGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAgg                                                                            >  2:4620655/1‑93 (MQ=255)

             tGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAgg                                                                            <  1:6193034/93‑1 (MQ=255)

             tGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTCCTCCCAATTTACCTAACACAAGGGGACAAATTTGCGATAAga                                                                            >  1:2824343/1‑92 (MQ=255)

              gACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGa                                                                           >  1:5735393/1‑93 (MQ=255)

               aCGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGAt                                                                          <  2:6043277/93‑1 (MQ=255)

                   ggggGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATaaaa                                                                      >  2:2366946/1‑93 (MQ=255)

                   gggTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATaaaa                                                                      >  1:9365639/1‑93 (MQ=255)

                    ggTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAg                                                                     >  2:2056236/1‑93 (MQ=255)

                     gTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACGTAACACAATGGGGCACACTTGCGATAAGGACAAAAGa                                                                    >  2:7557787/1‑93 (MQ=255)

                     gTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGa                                                                    <  1:7133089/93‑1 (MQ=255)

                      tGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGAc                                                                   >  1:8018826/1‑93 (MQ=255)

                       gAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGAcc                                                                  >  2:140362/1‑93 (MQ=255)

                       gAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGAcc                                                                  <  2:7634060/93‑1 (MQ=255)

                         aaaGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACcgc                                                                >  2:876088/1‑93 (MQ=255)

                           aGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGATACACTTGCGATAAGGATAAAAGACCGCGa                                                              <  1:9972579/93‑1 (MQ=255)

                           aGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGa                                                              <  1:2665586/93‑1 (MQ=255)

                            gACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGAt                                                             <  1:6627158/93‑1 (MQ=255)

                               cATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTAt                                                          >  2:7526748/1‑93 (MQ=255)

                               cATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTAt                                                          >  1:4684848/1‑93 (MQ=255)

                                aTCAAAAACCTCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATc                                                         <  1:5869704/93‑1 (MQ=255)

                                aTCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTA‑CTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCt                                                        <  1:1889449/93‑1 (MQ=255)

                                 tCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACCCACTTGCGATAAGGCTAAAAGACCGCGATTATCt                                                        >  1:2187233/1‑93 (MQ=255)

                                 tCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCt                                                        >  2:6177267/1‑93 (MQ=255)

                                  cAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCtt                                                       <  1:1227889/93‑1 (MQ=255)

                                   aaaaaCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCttt                                                      >  1:9593485/1‑93 (MQ=255)

                                    aaaaCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCtttt                                                     <  1:2063381/93‑1 (MQ=255)

                                     aaaCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTc                                                    <  2:3751443/93‑1 (MQ=255)

                                     aaaCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTc                                                    <  2:3041179/93‑1 (MQ=255)

                                     aaaCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTc                                                    <  2:7849203/93‑1 (MQ=255)

                                     aaaCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTc                                                    >  2:8539497/1‑93 (MQ=255)

                                      aaCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTcc                                                   <  1:8728310/93‑1 (MQ=255)

                                      aaCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTcc                                                   >  1:3681918/1‑93 (MQ=255)

                                       aCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCt                                                  <  2:474791/93‑1 (MQ=255)

                                       aCCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCt                                                  <  2:9603021/93‑1 (MQ=255)

                                        ccGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCtt                                                 >  1:3120005/1‑93 (MQ=255)

                                          gCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCa                                               <  1:6791761/93‑1 (MQ=255)

                                          gCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCa                                               >  2:3821572/1‑93 (MQ=255)

                                            tttCATTGTGCTATTGCCTGCTCTCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCAtt                                             <  2:5499034/93‑1 (MQ=255)

                                            tttCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCAtt                                             >  2:1506102/1‑93 (MQ=255)

                                             ttCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTg                                            <  1:7394820/93‑1 (MQ=255)

                                              tCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTga                                           >  1:641447/1‑93 (MQ=255)

                                               cATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTgag                                          <  2:532778/93‑1 (MQ=255)

                                               cATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTgag                                          <  1:3562189/93‑1 (MQ=255)

                                                aTTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTgaga                                         <  1:7549717/93‑1 (MQ=255)

                                                  tgtgCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATg                                       <  2:8074982/93‑1 (MQ=255)

                                                    tgCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCg                                     <  2:3587313/93‑1 (MQ=255)

                                                     gCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCgg                                    >  2:1717563/1‑93 (MQ=255)

                                                      cTATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCggg                                   <  2:9158085/93‑1 (MQ=255)

                                                       tATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGt                                  >  2:1131068/1‑93 (MQ=255)

                                                       tATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGt                                  <  2:6479497/93‑1 (MQ=255)

                                                       tATTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGt                                  >  2:4287126/1‑93 (MQ=255)

                                                        aTTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTa                                 >  1:8111094/1‑93 (MQ=255)

                                                        aTTGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTa                                 >  2:5035493/1‑93 (MQ=255)

                                                         ttGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTaa                                >  2:3909685/1‑93 (MQ=255)

                                                          tGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAAt                               >  1:1392596/1‑93 (MQ=255)

                                                          tGCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGATAAt                               >  1:4926573/1‑93 (MQ=255)

                                                           gCCTGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAAtt                              >  2:8322301/1‑93 (MQ=255)

                                                              tGCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATtctc                           >  1:5300977/1‑93 (MQ=255)

                                                              tGCTCCCAAATTACCTAAAACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATtctc                           >  2:4092430/1‑93 (MQ=255)

                                                               gCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCa                          <  2:9508083/93‑1 (MQ=255)

                                                               gCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCa                          <  2:2645406/93‑1 (MQ=255)

                                                               gCTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCa                          <  2:2547802/93‑1 (MQ=255)

                                                                cTCCCAAATTACCTAACACAAGGTGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCAt                         <  2:7236281/93‑1 (MQ=255)

                                                                cTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCAt                         <  1:6345227/93‑1 (MQ=255)

                                                                cTCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCAt                         >  2:6322180/1‑93 (MQ=255)

                                                                 tCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCAtt                        <  2:8848140/93‑1 (MQ=255)

                                                                 tCCCAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCAtt                        >  1:5452907/1‑93 (MQ=255)

                                                                  cccAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTa                       <  2:6767355/93‑1 (MQ=255)

                                                                    cAAATTACCTAACACAAGGGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAGTTCTCATTAAt                     >  1:9103917/1‑93 (MQ=255)

                                                                    cAAATTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAAt                     >  1:4264143/1‑93 (MQ=255)

                                                                     aaaTTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAAta                    >  2:7198866/1‑92 (MQ=255)

                                                                     aaaTTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATc                    <  1:8311078/93‑1 (MQ=255)

                                                                     aaaTTACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATc                    >  2:4482951/1‑93 (MQ=255)

                                                                        ttACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAAc                 <  1:669612/93‑1 (MQ=255)

                                                                         tACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACa                >  1:3826731/1‑93 (MQ=255)

                                                                         tACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACa                >  2:3337723/1‑93 (MQ=255)

                                                                         tACCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGAGTATCTTTTCCTTCATTGAGATGCGGGTAAATCTCATTAATCAACa                >  1:8932638/1‑93 (MQ=255)

                                                                          aCCTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAg               <  1:4860742/93‑1 (MQ=255)

                                                                            cTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGtt             >  2:6341868/1‑93 (MQ=255)

                                                                            cTAACACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGtt             >  1:1063695/1‑93 (MQ=255)

                                                                              aaCACAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCa           <  2:8303010/93‑1 (MQ=255)

                                                                               acacAAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCaa          >  1:4161616/1‑93 (MQ=255)

                                                                                cacaAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCAAt         >  1:5178242/1‑93 (MQ=255)

                                                                                cacaAGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCAAt         <  1:1585985/93‑1 (MQ=255)

                                                                                    aGGGGACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCAATAAAc     <  1:689807/93‑1 (MQ=255)

                                                                                     ggggACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCAATAAACg    >  2:8676966/1‑93 (MQ=255)

                                                                                       ggACACACTTGCGGTCAGGATACGAGACCGCGATTATCTTTTCCTACATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCAATAGACGGa  <  1:8976785/93‑1 (MQ=255)

                                                                                       ggACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGTTAATTCTCATTAATCAACAGTTCAATAAACGGa  <  2:4938693/93‑1 (MQ=255)

                                                                                       ggACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCAATAAACGGa  <  1:2757697/93‑1 (MQ=255)

                                                                                       ggACACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCAATAAACGGa  <  2:3387106/93‑1 (MQ=255)

CTTTTATCCAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCTCATTAATCAACAGTTCAATAAACGGA  >  NC_000913/1789545‑1789724

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