BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,751,393	G→A	T154T (ACC→ACT)	b1673 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,751,393	0	G	A	99.6%	857.3 / NA	229	T154T (ACC→ACT)	b1673	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (105/123); minor base C (1/0); total (106/123)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

GCCGCTGTGGTTACGGCTATTCAACATGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTGGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCTTTAA  >  NC_000913/1751303‑1751481

 ccgcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTc                                                                                       >  2:6069583/1‑93 (MQ=255)

 ccgcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTc                                                                                       >  1:4540759/1‑93 (MQ=255)

   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                                                                                     >  2:3657180/1‑93 (MQ=255)

     tgtgGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATc                                                                                   >  1:1219032/1‑93 (MQ=255)

     tgtgGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATc                                                                                   <  1:5587476/93‑1 (MQ=255)

     tgtgGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGGTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATc                                                                                   >  1:1659556/1‑93 (MQ=255)

      gtgGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCAGCCACACAGGTTTCCGGACTAGTCAATCg                                                                                  >  2:4630529/1‑93 (MQ=255)

      gtgGTTACGGCTATTCAACGTGACAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCg                                                                                  >  2:5788375/1‑93 (MQ=255)

         gTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACa                                                                               <  2:3211773/93‑1 (MQ=255)

         gTTACGGCTATTCAACGTCCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACa                                                                               <  1:5308536/93‑1 (MQ=255)

          ttACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACaa                                                                              <  2:613491/93‑1 (MQ=255)

           tACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACaaa                                                                             >  1:1600240/1‑93 (MQ=255)

            aCGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAt                                                                            <  2:6034703/93‑1 (MQ=255)

             cGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAtt                                                                           >  1:6028611/1‑93 (MQ=255)

             cGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAtt                                                                           <  1:2207146/93‑1 (MQ=255)

               gCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttc                                                                         <  2:5825569/93‑1 (MQ=255)

                cTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttct                                                                        <  1:1576097/93‑1 (MQ=255)

                  aTTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttcttc                                                                      <  1:6230711/93‑1 (MQ=255)

                  aTTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttcttc                                                                      <  1:46867/93‑1 (MQ=255)

                  aTTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttcttc                                                                      <  1:1782102/93‑1 (MQ=255)

                   ttCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTcc                                                                     <  2:2106646/93‑1 (MQ=255)

                   ttCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTcc                                                                     <  2:2371359/93‑1 (MQ=255)

                   ttCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTcc                                                                     <  1:4394132/93‑1 (MQ=255)

                   ttCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTcc                                                                     <  2:6257473/93‑1 (MQ=255)

                    tCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTATTCcg                                                                    <  2:3355121/93‑1 (MQ=255)

                     cAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCcgt                                                                   <  2:3013120/93‑1 (MQ=255)

                      aaCGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCcgtc                                                                  <  2:2908471/93‑1 (MQ=255)

                        cgtGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCcgtcgt                                                                <  1:2368808/91‑1 (MQ=255)

                         gtGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGtt                                                               >  1:5855099/2‑93 (MQ=255)

                          tGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTg                                                              >  2:5016291/1‑93 (MQ=255)

                           gCCAGAGATAGGAACAAGGATTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTgc                                                             >  1:6560415/1‑93 (MQ=255)

                            ccAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTgcg                                                            <  1:2778555/93‑1 (MQ=255)

                              agagagagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTgcgcg                                                          >  2:5605998/1‑93 (MQ=255)

                               gagagagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTgcgcgc                                                         >  1:93795/1‑93 (MQ=255)

                                agtgagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGcgcgcg                                                        >  1:5376167/4‑93 (MQ=255)

                                 gagagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGt                                                       >  2:3530212/1‑93 (MQ=255)

                                  agagGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGTTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTc                                                      <  1:2887342/93‑1 (MQ=255)

                                  agagGAACAAGGTTTTCCCCAGCATGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTc                                                      <  2:3174113/93‑1 (MQ=255)

                                    gggAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTccc                                                    <  1:2588803/92‑1 (MQ=255)

                                    agGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTccc                                                    >  2:4217127/1‑93 (MQ=255)

                                    agGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTccc                                                    <  2:2707631/93‑1 (MQ=255)

                                    agGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTccc                                                    <  1:3328542/93‑1 (MQ=255)

                                     ggAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTTCGCGCGTCGCt                                                   >  1:3977537/1‑93 (MQ=255)

                                     ggAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCt                                                   <  1:2579041/93‑1 (MQ=255)

                                     gcaaCAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCt                                                   <  1:1016702/91‑1 (MQ=255)

                                      gAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCtt                                                  <  1:862068/93‑1 (MQ=255)

                                       aaCAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCttt                                                 <  2:896259/93‑1 (MQ=255)

                                       aaCAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCttt                                                 <  1:4296239/93‑1 (MQ=255)

                                         cAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTc                                               <  2:370523/93‑1 (MQ=255)

                                          aaGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTcc                                              >  2:116996/1‑93 (MQ=255)

                                          aaGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTcc                                              >  2:5434494/1‑93 (MQ=255)

                                             gTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCa                                           <  2:1474165/93‑1 (MQ=255)

                                             gTTTTCACCAGCCTGACCAATAGCAGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCa                                           >  1:5559385/1‑93 (MQ=255)

                                              ttttCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCAt                                          <  2:1565524/93‑1 (MQ=255)

                                               tttCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATa                                         <  1:2439784/93‑1 (MQ=255)

                                                  ccccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGa                                      >  2:3785800/1‑93 (MQ=255)

                                                   cccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGaa                                     <  1:1412911/93‑1 (MQ=255)

                                                   cccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGaa                                     <  2:3292283/93‑1 (MQ=255)

                                                   cccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGaa                                     <  2:1826259/93‑1 (MQ=255)

                                                    ccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGaaa                                    >  2:4909355/1‑93 (MQ=255)

                                                    ccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGaaa                                    >  1:5021092/1‑93 (MQ=255)

                                                    ccAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGaaa                                    <  2:4854732/93‑1 (MQ=255)

                                                      aGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATc                                  <  2:4604681/93‑1 (MQ=255)

                                                      aGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATc                                  >  1:2651129/1‑93 (MQ=255)

                                                       gCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCg                                 >  2:4387696/1‑93 (MQ=255)

                                                       gCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCg                                 >  1:2273019/1‑93 (MQ=255)

                                                        ccTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCgg                                >  1:1721985/1‑93 (MQ=255)

                                                         cTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGc                               <  1:1088926/93‑1 (MQ=255)

                                                         cTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGc                               >  2:3519311/1‑93 (MQ=255)

                                                         cTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGc                               <  2:1913691/93‑1 (MQ=255)

                                                          tGACCAATAGCCGCCACACAGGTTTCCGGACTCGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATACGAAATCGGCt                              >  1:647468/1‑93 (MQ=255)

                                                          tGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCt                              >  1:3045438/1‑93 (MQ=255)

                                                           gACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtt                             <  2:4080649/93‑1 (MQ=255)

                                                           gACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtt                             >  1:5277467/1‑93 (MQ=255)

                                                            aCCAATAGCCGCCACACATGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCttt                            <  2:833587/93‑1 (MQ=255)

                                                            aCCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCttt                            <  1:4818762/93‑1 (MQ=255)

                                                             ccAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtttt                           <  1:4796973/93‑1 (MQ=255)

                                                             ccAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCtttt                           <  1:714777/93‑1 (MQ=255)

                                                               aaTAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTc                         >  2:3320585/1‑93 (MQ=255)

                                                                aTAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCACTTTTCCCCATAAGAAATCGGCTTTTTcc                        >  1:4970200/1‑93 (MQ=255)

                                                                 tAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCa                       <  1:3601584/93‑1 (MQ=255)

                                                                 tAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCa                       <  1:2760867/93‑1 (MQ=255)

                                                                  aGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAg                      <  2:19967/93‑1 (MQ=255)

                                                                   gccgccACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAgg                     <  2:4702337/93‑1 (MQ=255)

                                                                     cgccACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCt                   <  2:4727512/93‑1 (MQ=255)

                                                                      gccACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCCTCCCTTTTCCCACTAAGAAATCGGCTTTTTCCAGGCTa                  >  1:6677543/1‑93 (MQ=255)

                                                                        cacacaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAAc                <  1:6165282/93‑1 (MQ=255)

                                                                         acacaGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACt               >  1:41507/1‑93 (MQ=255)

                                                                            caGGTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTg            <  2:2472928/93‑1 (MQ=255)

                                                                              ggTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTTCAGGCTAACTTTGTc          >  2:5729741/1‑93 (MQ=255)

                                                                              ggTTTCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTc          <  2:5525206/93‑1 (MQ=255)

                                                                               gctgccGGACTAGTCAATCGACAAATTGCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCa         <  2:6436885/89‑1 (MQ=255)

                                                                                 ttCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATc       <  2:6558851/93‑1 (MQ=255)

                                                                                  tCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCt      <  1:6115189/93‑1 (MQ=255)

                                                                                  tCCGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCt      <  2:4945252/93‑1 (MQ=255)

                                                                                   ccGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCtt     <  1:4047138/93‑1 (MQ=255)

                                                                                    cGGACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCttt    <  1:2044519/93‑1 (MQ=255)

                                                                                      gACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTTCAGGCTAACTTTGTCATCTTTaa  >  1:5058031/1‑93 (MQ=255)

                                                                                      gACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCTTTaa  <  2:1629191/93‑1 (MQ=255)

                                                                                      gACTAGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCTTTaa  <  1:5758140/93‑1 (MQ=255)

GCCGCTGTGGTTACGGCTATTCAACATGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTGGTCAATCGACAAATTTCTTCCGTCGTTGCGCGCGTCCCTTTTCCCCATAAGAAATCGGCTTTTTCCAGGCTAACTTTGTCATCTTTAA  >  NC_000913/1751303‑1751481

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