BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,130,334	A→G	V30V (GTT→GTC)	b2983 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,130,334	0	A	G	99.7%	937.9 / ‑5.2	286	V30V (GTT→GTC)	b2983	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (130/155); total (131/155)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 327 total aligned reads displayed.

TACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACAACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGCCAA  >  NC_000913/3130242‑3130423

  cGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGAACCAGCAAGGCGCTCacgac                                                                                         >  1:3764906/1‑93 (MQ=255)

   gATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCacgaca                                                                                        <  2:2193824/93‑1 (MQ=255)

   gATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCACCAAGGCGCTCacgaca                                                                                        <  2:6708500/93‑1 (MQ=255)

   gATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCCTTAACGGCACACGAGAGCGCC‑CACAGACCCAGCCAAGCGCTCacgaca                                                                                        >  2:816391/1‑93 (MQ=255)

    aTTACACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACAt                                                                                       <  2:3306361/93‑1 (MQ=255)

     ttGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACAtt                                                                                      >  1:6228102/1‑93 (MQ=255)

     ttGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACAtt                                                                                      <  2:6972399/93‑1 (MQ=255)

     ttGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACAta                                                                                      <  1:4480031/93‑2 (MQ=255)

         acaATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTCCCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTg                                                                                  <  1:9393224/93‑1 (MQ=255)

          caATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTgg                                                                                 >  2:6848306/1‑93 (MQ=255)

            aTCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAg                                                                               >  2:5905083/1‑93 (MQ=255)

              cACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGcc                                                                             >  1:7478614/1‑93 (MQ=255)

                  aGTACGCCAACCATGGCTGCCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAAc                                                                         <  2:1195473/93‑1 (MQ=255)

                  aGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAAc                                                                         <  1:1024760/93‑1 (MQ=255)

                  aGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAAc                                                                         >  1:7952170/1‑93 (MQ=255)

                   gTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACCTTACTGGAGCCTAAcc                                                                        >  1:9184359/1‑93 (MQ=255)

                    tACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCa                                                                       >  2:9650506/1‑93 (MQ=255)

                     acaccAAACATGGCGGGCGTCATCCCTTTCCCGGCACAGGACAGCCCA‑AACCGCCCCAGAAAGGCCCGCACGACATTACTGGAGCCTAACCAg                                                                      <  2:7567092/90‑1 (MQ=255)

                     aCGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAg                                                                      <  2:8775682/93‑1 (MQ=255)

                       gCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGCCCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGgc                                                                    <  2:3357657/93‑1 (MQ=255)

                         cAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGgcgc                                                                  >  1:9555684/1‑93 (MQ=255)

                          aaaCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCt                                                                 >  2:3695495/1‑93 (MQ=255)

                            aCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAt                                                               <  1:6000465/93‑1 (MQ=255)

                             cATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAtt                                                              <  2:1054810/93‑1 (MQ=255)

                              aTGGCGGGCGTCCTCCCCTTACCGGCACACGAGCGCGCC‑AACAGAACCAGCAAGGCGCTCACGCCATTACTGGAGCCTAACAAGGCGCTATTg                                                             >  2:7112332/1‑93 (MQ=255)

                               tGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGc                                                            <  1:3590792/93‑1 (MQ=255)

                                 gcgggcgTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGa                                                          <  2:1434121/93‑1 (MQ=255)

                                    ggcgTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaata                                                       >  2:2484787/1‑93 (MQ=255)

                                     gcgTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataa                                                      <  2:1706269/93‑1 (MQ=255)

                                     gcgTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCCCTATTGCGaataa                                                      <  1:6721026/93‑1 (MQ=255)

                                      cttCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataat                                                     <  2:8095837/91‑1 (MQ=255)

                                      cgTCATCCCTTTCCCGGCAGACGAGAGCTCC‑AACAGACCCAGCAAGGGGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataat                                                     <  2:2659168/93‑1 (MQ=255)

                                      cgTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataat                                                     <  2:4957076/93‑1 (MQ=255)

                                      cgTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataat                                                     <  1:3027567/93‑1 (MQ=255)

                                       gTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATCATg                                                    >  1:2408600/1‑93 (MQ=255)

                                       gTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATg                                                    <  1:3422486/93‑1 (MQ=255)

                                       gTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATg                                                    <  1:1539176/93‑1 (MQ=255)

                                       cccttCCTTTTCCCGGCACAGAAGACCCAA‑AACAGCCCAAGCAGGGCGCCCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATg                                                    <  1:3858089/89‑1 (MQ=255)

                                        tCATCCCTTTCCCGGCACACGAGAGCGCC‑AACAGACCCAGCAGGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATgg                                                   <  1:8834870/93‑1 (MQ=255)

                                        tCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATgg                                                   >  2:5468465/1‑93 (MQ=255)

                                         cATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGt                                                  >  2:9329750/1‑93 (MQ=255)

                                         cATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGt                                                  <  2:9903128/93‑1 (MQ=255)

                                           tCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGc                                                >  2:4050232/1‑93 (MQ=255)

                                           tCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGc                                                >  1:5704027/1‑93 (MQ=255)

                                           ccccTTTCCCGGCCCACG‑GAGCCCCAAACAGCCCCAGCAAGGCCCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGc                                                <  2:7758347/92‑1 (MQ=255)

                                            cccTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCg                                               >  2:1762832/1‑93 (MQ=255)

                                             ccTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGa                                              >  1:9703231/1‑93 (MQ=255)

                                             ccTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGa                                              >  2:9741636/1‑93 (MQ=255)

                                               tttACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGaaa                                            >  1:3731005/1‑93 (MQ=255)

                                                ttACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAAc                                           >  2:1364304/1‑93 (MQ=255)

                                                  aCCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAAcgc                                         <  2:9528886/93‑1 (MQ=255)

                                                    cGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACgcgc                                       <  2:6201966/93‑1 (MQ=255)

                                                     ggCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGcc                                      >  1:8401398/1‑93 (MQ=255)

                                                     ggCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGcc                                      <  2:7896307/93‑1 (MQ=255)

                                                       cacaCGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCAt                                    >  2:6504409/1‑93 (MQ=255)

                                                         cacGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCt                                  >  2:7500083/1‑93 (MQ=255)

                                                           cGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGc                                >  1:8224164/1‑93 (MQ=255)

                                                            gagaGCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCa                               >  1:1092182/1‑93 (MQ=255)

                                                             agagCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaa                              >  2:8721848/1‑93 (MQ=255)

                                                             agagCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaa                              >  2:5175830/1‑93 (MQ=255)

                                                             agagCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaa                              <  1:3501828/93‑1 (MQ=255)

                                                              gagCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaaa                             >  2:2823975/1‑93 (MQ=255)

                                                              gagCGCC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaaa                             >  1:1984056/1‑93 (MQ=255)

                                                                 cgcC‑AACAGCCCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCa                          <  2:553882/93‑1 (MQ=255)

                                                                  gcC‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCaa                         <  1:7717340/93‑1 (MQ=255)

                                                                   cc‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAAt                        <  2:3099589/93‑1 (MQ=255)

                                                                   cc‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAAt                        <  1:7836834/93‑1 (MQ=255)

                                                                   cc‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAAt                        >  1:7138388/1‑93 (MQ=255)

                                                                    c‑AACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATg                       <  2:7042300/93‑1 (MQ=255)

                                                                       aCAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGtt                     <  2:8503220/93‑1 (MQ=255)

                                                                        cAGCCCCAGCAAGCCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGttt                    <  2:434530/93‑1 (MQ=255)

                                                                        cAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGttt                    <  1:1192511/93‑1 (MQ=255)

                                                                        cAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGttt                    <  1:1715265/93‑1 (MQ=255)

                                                                         agccccAGCACGGCGCCCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTg                   <  2:5787737/90‑1 (MQ=255)

                                                                         aGACCCAGCAAGTCCCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTg                   <  1:2837601/93‑1 (MQ=255)

                                                                         aGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTg                   >  2:3461957/1‑93 (MQ=255)

                                                                          gACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGa                  >  1:6983387/1‑93 (MQ=255)

                                                                          gACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGa                  >  2:8238796/1‑93 (MQ=255)

                                                                          agcccGGCAAGGCCCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGa                  <  1:2155416/91‑1 (MQ=255)

                                                                                gCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaa            >  1:8671337/1‑93 (MQ=255)

                                                                                gCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaa            <  2:8796973/93‑1 (MQ=255)

                                                                                 cgaGGCCCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaa           <  2:1896739/91‑1 (MQ=255)

                                                                                 cAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaa           >  1:9524588/1‑93 (MQ=255)

                                                                                 cAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaa           >  1:5307313/1‑93 (MQ=255)

                                                                                 cAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaa           <  1:4921196/93‑1 (MQ=255)

                                                                                 cAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaa           >  1:2963029/1‑93 (MQ=255)

                                                                                 cAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaa           <  2:6402208/93‑1 (MQ=255)

                                                                                 aaaGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaa           <  1:4092408/92‑1 (MQ=255)

                                                                                  aaGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaaa          >  2:4819734/1‑93 (MQ=255)

                                                                                  aaGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTaaaa          >  1:4588519/1‑93 (MQ=255)

                                                                                   agccgcTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAAc         <  2:7281087/90‑1 (MQ=255)

                                                                                   aGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAAc         >  2:6109468/1‑93 (MQ=255)

                                                                                    ggCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAAcc        <  1:2240145/93‑1 (MQ=255)

                                                                                    ggCGCCCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAAcc        <  1:5778673/93‑1 (MQ=255)

                                                                                      cgcTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCg      <  1:5922228/93‑1 (MQ=255)

                                                                                        cTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGcc    >  2:4752922/1‑93 (MQ=255)

                                                                                        cTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGcc    >  2:9061272/1‑93 (MQ=255)

                                                                                        cTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGcc    <  2:943236/93‑1 (MQ=255)

                                                                                         tCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGCCa   >  1:4720627/1‑93 (MQ=255)

                                                                                         tCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGCCa   >  2:6308688/1‑93 (MQ=255)

TACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCC‑AACAGACCCAGCAAGGCGCTCACAACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTTGATGTTAAAACCCGCCAA  >  NC_000913/3130242‑3130423

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