BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,917,590	T→C	V688A (GTC→GCC)	b1834 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,917,590	0	T	C	99.1%	341.1 / NA	107	V688A (GTC→GCC)	b1834	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (54/52); total (54/53)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

TCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTC  >  NC_000913/1917500‑1917678

tcGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGccc                                                                                        <  1:150893/93‑1 (MQ=255)

 cGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCa                                                                                       >  2:1984210/1‑93 (MQ=255)

     aTTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGAcc                                                                                   >  2:2331905/1‑93 (MQ=255)

     aTTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGAcc                                                                                   >  1:2168735/1‑93 (MQ=255)

      ttGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCt                                                                                  >  2:4080007/1‑93 (MQ=255)

       tGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCtt                                                                                 <  1:786196/93‑1 (MQ=255)

       tGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCtt                                                                                 <  1:3787591/93‑1 (MQ=255)

       tGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCtt                                                                                 <  2:915996/93‑1 (MQ=255)

       tGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCtt                                                                                 <  2:692493/93‑1 (MQ=255)

       tGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCtt                                                                                 <  2:692494/93‑1 (MQ=255)

       tGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCtt                                                                                 >  1:2404153/1‑93 (MQ=255)

         atatCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTg                                                                               <  1:1034333/93‑1 (MQ=255)

         atatCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTg                                                                               >  1:338012/1‑93 (MQ=255)

             cGGGCAAATCCAGACGCTGTATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTcg                                                                           >  2:2199493/1‑93 (MQ=255)

             cGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTcg                                                                           <  2:1498638/93‑1 (MQ=255)

               ggCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTcgcg                                                                         >  1:2934711/1‑93 (MQ=255)

               ggCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTcgcg                                                                         >  2:3121659/1‑93 (MQ=255)

               ggCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTcgcg                                                                         <  1:1828057/93‑1 (MQ=255)

                 cAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATNNGCCCAGACCTTTGCTCGCngg                                                                       >  1:753188/1‑92 (MQ=255)

                  aaaTCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCggtg                                                                      <  2:1512132/93‑1 (MQ=255)

                  aaaTCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCggtg                                                                      >  1:4221134/1‑93 (MQ=255)

                  aaaTCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCggtg                                                                      >  2:3386557/1‑93 (MQ=255)

                    aTCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCggtggt                                                                    <  2:3558430/93‑1 (MQ=255)

                     tCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTa                                                                   <  2:2559521/93‑1 (MQ=255)

                     tCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTa                                                                   <  1:2766857/93‑1 (MQ=255)

                       cAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTAcg                                                                 >  2:3059382/1‑93 (MQ=255)

                        aGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTAcgc                                                                <  1:4564579/93‑1 (MQ=255)

                        aGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTAcgc                                                                <  2:194271/93‑1 (MQ=255)

                        aGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTAcgc                                                                <  2:4563281/93‑1 (MQ=255)

                            gCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCtt                                                            >  1:679492/1‑93 (MQ=255)

                            gCTGGATCTGATTACCGCACGCAATGAAGGACAGGCAAAGGCGGTACTTTACCAGGAATATGCGCAGACCTTTGCTCGCGGTGGTACGCGCtt                                                            >  1:679489/1‑93 (MQ=255)

                             cTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTtc                                                           <  2:575697/93‑1 (MQ=255)

                             cTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGCACGCGCTtc                                                           >  1:3766723/1‑93 (MQ=255)

                              tGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTtct                                                          <  2:4509834/93‑1 (MQ=255)

                              tGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTtct                                                          >  1:4576888/1‑93 (MQ=255)

                               tgATCTTATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTGCCCGGAATATGCCCAAACCTTTGCTCGCGGTGGTACTCGCTtctc                                                         <  2:3322777/92‑1 (MQ=255)

                               ggATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTtctc                                                         >  2:446076/1‑93 (MQ=255)

                               ggATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTtctc                                                         <  1:2558084/93‑1 (MQ=255)

                               ggATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGGGGTACGCGCTtctc                                                         >  2:2669303/1‑93 (MQ=255)

                                 aTCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAg                                                       <  1:4005657/93‑1 (MQ=255)

                                 aTCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAg                                                       <  2:4388706/93‑1 (MQ=255)

                                  tCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGt                                                      <  2:813619/93‑1 (MQ=255)

                                  tCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGt                                                      <  1:3112207/93‑1 (MQ=255)

                                   cTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTg                                                     <  1:2188193/93‑1 (MQ=255)

                                    tGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTgg                                                    <  2:991854/93‑1 (MQ=255)

                                    tGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTgg                                                    >  1:2430611/1‑93 (MQ=255)

                                     gATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGt                                                   <  2:777630/93‑1 (MQ=255)

                                       ttACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTca                                                 <  2:4372681/93‑1 (MQ=255)

                                       ttACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTca                                                 >  1:500019/1‑93 (MQ=255)

                                         aCCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTcaca                                               <  2:1061932/93‑1 (MQ=255)

                                          ccGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCacac                                              <  2:3621135/93‑1 (MQ=255)

                                           cgcacgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACAcc                                             <  2:1268219/93‑1 (MQ=255)

                                           cgcacgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACAcc                                             >  2:2981858/1‑93 (MQ=255)

                                              acgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCa                                          >  2:3874296/1‑93 (MQ=255)

                                              acgcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCa                                          >  1:1690461/1‑93 (MQ=255)

                                                gcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCaca                                        >  2:302990/1‑91 (MQ=255)

                                                gcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCaaa                                        >  2:4567328/1‑93 (MQ=255)

                                                gcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCaaa                                        <  2:3511589/93‑1 (MQ=255)

                                                gcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCaaa                                        >  1:257609/1‑93 (MQ=255)

                                                gcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCaaa                                        <  1:3607909/93‑1 (MQ=255)

                                                gcaATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCaaa                                        >  1:3809059/1‑93 (MQ=255)

                                                 caATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAAt                                       >  1:1750539/1‑93 (MQ=255)

                                                  aaTGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAAtt                                      >  2:2813786/1‑93 (MQ=255)

                                                  aaTGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAAtt                                      <  1:2016020/93‑1 (MQ=255)

                                                           cAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCt                             >  1:265180/1‑93 (MQ=255)

                                                           cAGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCt                             <  1:1950024/93‑1 (MQ=255)

                                                            aGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTg                            >  1:4394757/1‑93 (MQ=255)

                                                            aGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTg                            <  2:3758151/93‑1 (MQ=255)

                                                            aGGCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTg                            >  1:1307630/1‑93 (MQ=255)

                                                             ggCAAAGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTgg                           <  2:429209/93‑1 (MQ=255)

                                                                 aaGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGtt                       <  2:2769521/93‑1 (MQ=255)

                                                                 aaGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGtt                       >  2:2085744/1‑93 (MQ=255)

                                                                  aGGCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTg                      <  2:1839799/93‑1 (MQ=255)

                                                                   ggCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGa                     <  1:3312383/93‑1 (MQ=255)

                                                                    gCGGTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAg                    <  1:1109152/93‑1 (MQ=255)

                                                                      ggTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCa                  <  1:4054427/93‑1 (MQ=255)

                                                                       gTACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCAt                 >  2:3826391/1‑93 (MQ=255)

                                                                        tACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATc                >  1:1150096/1‑93 (MQ=255)

                                                                        tACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATc                >  1:1874468/1‑93 (MQ=255)

                                                                        tACTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATc                >  1:2559049/1‑93 (MQ=255)

                                                                         aCTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCt               >  1:3017026/1‑93 (MQ=255)

                                                                          cTTTACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCtt              <  2:1789116/93‑1 (MQ=255)

                                                                           tttACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTg             >  2:2273924/1‑93 (MQ=255)

                                                                            ttACCCGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGATGAGCATCTTGa            >  2:3517800/1‑93 (MQ=255)

                                                                               cccGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATAc         >  2:1618195/1‑93 (MQ=255)

                                                                               cccGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATAc         >  1:594714/1‑93 (MQ=255)

                                                                               cccGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCACAGTGGTCACACCGCAAATTTCGGCAGCTGGCATTGAGCATCTTGATAc         <  2:4217111/93‑1 (MQ=255)

                                                                                 cGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTa       <  1:1636141/93‑1 (MQ=255)

                                                                                 cGGAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTa       >  2:743964/1‑93 (MQ=255)

                                                                                  ggAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAt      >  2:4459496/1‑93 (MQ=255)

                                                                                   gAATATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtc     <  1:945828/93‑1 (MQ=255)

                                                                                    aaTATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcc    <  2:1907329/93‑1 (MQ=255)

                                                                                    aaTATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcc    >  1:3677779/1‑93 (MQ=255)

                                                                                    aaTATGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcc    <  2:1304239/93‑1 (MQ=255)

                                                                                     atatGCCCAGACCTTTGCTCGCTGAGGTACTCCCTTCTCCGTGGTCACACCGCAAATTTCGGCAGGTTGCGTTGAGCCTCTTGATACTAtgct   >  2:3549948/1‑93 (MQ=255)

                                                                                     atatGCCCAGACCTTTGCTCGCGGTGGTCCGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcct   >  2:4480628/1‑93 (MQ=255)

                                                                                     atatGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcct   >  2:4480629/1‑93 (MQ=255)

                                                                                     atatGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcct   >  2:915421/1‑93 (MQ=255)

                                                                                      tatGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcctc  >  2:864259/1‑93 (MQ=255)

                                                                                      tatGCCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcctc  >  2:2231043/1‑93 (MQ=255)

TCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTC  >  NC_000913/1917500‑1917678

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