BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,683,091	A→G	V42A (GTC→GCC)	yhjK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,683,091	0	A	G	99.8%	1262.3 / ‑5.7	400	V42A (GTC→GCC)	yhjK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (227/171); total (229/171)
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.94e-01

Only 100 of 409 total aligned reads displayed.

CTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAA  >  NC_000913/3683001‑3683180

cTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCa                                                                                         <  2:1126280/93‑1 (MQ=255)

 tgcccACGCCCGCCGGGTTAATGCTGGCAAGAATAGCTGCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCag                                                                                        <  1:4926202/90‑1 (MQ=255)

 tGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCag                                                                                        <  2:3638811/93‑1 (MQ=255)

  gACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATATCTGAAGATAAGGGTTCGCGAGCaga                                                                                       >  1:6536690/1‑93 (MQ=255)

  gACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCaga                                                                                       <  2:5965166/93‑1 (MQ=255)

      cgcccTGCCGTTTTAAGGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcga                                                                                   <  1:1105870/90‑1 (MQ=255)

      ccgcgTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcga                                                                                   <  2:2852526/92‑1 (MQ=255)

      aCGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcga                                                                                   >  2:2737791/1‑93 (MQ=255)

        gCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcgagc                                                                                 >  1:5000526/1‑93 (MQ=255)

        gCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcgagc                                                                                 <  1:3476381/93‑1 (MQ=255)

        gCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcgagc                                                                                 >  2:615454/1‑93 (MQ=255)

           tGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAAt                                                                              >  1:7854539/1‑93 (MQ=255)

           tGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAAt                                                                              >  2:1720843/1‑93 (MQ=255)

           tGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAAt                                                                              >  2:4254812/1‑93 (MQ=255)

           tGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAAt                                                                              >  1:1290755/1‑93 (MQ=255)

           tGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAAt                                                                              >  1:6251886/1‑93 (MQ=255)

           tGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAAt                                                                              >  1:6353607/1‑93 (MQ=255)

               ggTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCtt                                                                          <  1:3822858/93‑1 (MQ=255)

               ggTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCtt                                                                          <  2:757404/93‑1 (MQ=255)

               ggTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCtt                                                                          <  2:2215775/93‑1 (MQ=255)

                gTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCttt                                                                         >  2:681671/1‑93 (MQ=255)

                 tttAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTc                                                                        >  1:6690318/1‑93 (MQ=255)

                 tttAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGATCGCTGGCAGAGCGAGCAATGCTTTc                                                                        >  1:7009692/1‑93 (MQ=255)

                    aaTGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAg                                                                     >  1:2589288/1‑93 (MQ=255)

                    aaTGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAg                                                                     >  2:7215586/1‑93 (MQ=255)

                     aTGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGt                                                                    <  1:4327254/93‑1 (MQ=255)

                     aTGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGt                                                                    >  1:3783773/1‑93 (MQ=255)

                     aTGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGt                                                                    <  1:2433556/93‑1 (MQ=255)

                      tGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGtt                                                                   >  1:7459328/1‑93 (MQ=255)

                         tGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCg                                                                >  2:1906961/1‑93 (MQ=255)

                         tGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCg                                                                >  2:6421081/1‑93 (MQ=255)

                         tGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCg                                                                >  1:1884650/1‑93 (MQ=255)

                         tGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCGTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCg                                                                >  2:4954996/1‑93 (MQ=255)

                          ggCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGt                                                               >  1:4531074/1‑93 (MQ=255)

                          ggCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGt                                                               >  2:1297294/1‑93 (MQ=255)

                           gCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTa                                                              <  2:3725368/93‑1 (MQ=255)

                             aaGAATAGCTTCCGCTTCGGGAATATCGCATTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGc                                                            <  1:5390382/93‑1 (MQ=255)

                              aGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGcc                                                           >  2:3832015/1‑93 (MQ=255)

                               gAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCg                                                          >  1:5842957/1‑93 (MQ=255)

                               gAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCg                                                          >  1:1527599/1‑93 (MQ=255)

                                aaTAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCgt                                                         >  2:862841/1‑93 (MQ=255)

                                aaTAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCgt                                                         >  2:6566152/1‑93 (MQ=255)

                                 aTAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCgtg                                                        >  1:1170789/1‑93 (MQ=255)

                                  tAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGCAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTTCAGTTGCGTAGTCgtgt                                                       >  2:7079082/1‑93 (MQ=255)

                                    gCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTta                                                     <  2:73030/93‑1 (MQ=255)

                                    gCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTta                                                     >  2:7917512/1‑93 (MQ=255)

                                    gCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTta                                                     >  1:5186056/1‑93 (MQ=255)

                                     cTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTtat                                                    <  1:4139748/93‑1 (MQ=255)

                                      ttCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTtata                                                   <  1:7691752/93‑1 (MQ=255)

                                      ttCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTtata                                                   >  2:1176728/1‑93 (MQ=255)

                                       tCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAg                                                  >  2:610422/1‑93 (MQ=255)

                                        ccGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGc                                                 <  2:4514205/93‑1 (MQ=255)

                                         cGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCg                                                <  2:5301158/93‑1 (MQ=255)

                                            ttCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAtt                                             >  1:3195613/1‑93 (MQ=255)

                                            ttCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAtt                                             >  1:4730138/1‑93 (MQ=255)

                                            ttCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAtt                                             >  2:548097/1‑93 (MQ=255)

                                             tCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTc                                            <  1:791179/93‑1 (MQ=255)

                                                ggAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTctgc                                         <  1:6747964/93‑1 (MQ=255)

                                                ggAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTctgc                                         <  2:4769813/93‑1 (MQ=255)

                                                  aaTATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTctgctg                                       <  1:3907940/93‑1 (MQ=255)

                                                   atatCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTgg                                      <  1:6008661/93‑1 (MQ=255)

                                                    tatCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGa                                     >  1:2716532/1‑93 (MQ=255)

                                                    tatCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGa                                     >  2:6691257/1‑93 (MQ=255)

                                                    tatCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGa                                     >  2:3035531/1‑93 (MQ=255)

                                                     atCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGAc                                    >  2:1040614/1‑93 (MQ=255)

                                                       cGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCa                                  <  2:3975643/93‑1 (MQ=255)

                                                          cTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGAt                               >  2:2924772/1‑93 (MQ=255)

                                                          cTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGAt                               <  2:7453070/93‑1 (MQ=255)

                                                          cTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGAt                               <  1:104417/93‑1 (MQ=255)

                                                          cTTTCAAAATAGCTGAAGATAAGCGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGAt                               <  2:4304589/93‑1 (MQ=255)

                                                            ttCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATgg                             >  1:6129859/1‑93 (MQ=255)

                                                              cAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGaa                           <  2:7680641/93‑1 (MQ=255)

                                                              cAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGaa                           <  1:2329177/93‑1 (MQ=255)

                                                               aaaaTAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAAc                          >  2:2393898/1‑93 (MQ=255)

                                                               aaaaTAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAAc                          >  2:5586399/1‑93 (MQ=255)

                                                                aaaTAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACa                         <  1:1660772/93‑1 (MQ=255)

                                                                aaaTAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACa                         <  1:7341636/93‑1 (MQ=255)

                                                                 aaTAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAg                        >  2:5150381/1‑93 (MQ=255)

                                                                   tAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCa                      >  1:4045330/1‑93 (MQ=255)

                                                                   tAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCa                      >  2:5879992/1‑93 (MQ=255)

                                                                     gCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCaaa                    >  2:1644819/1‑93 (MQ=255)

                                                                       tGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAAc                  >  2:1540408/1‑93 (MQ=255)

                                                                       tGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAAc                  >  1:2850973/1‑93 (MQ=255)

                                                                        gAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACg                 <  1:4283903/93‑1 (MQ=255)

                                                                        gAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACg                 >  2:4593385/1‑93 (MQ=255)

                                                                        gAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACg                 >  1:863511/1‑93 (MQ=255)

                                                                        gAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACg                 >  1:2981847/1‑93 (MQ=255)

                                                                         aaGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACgg                >  1:2198374/1‑93 (MQ=255)

                                                                           gATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTg              >  2:2455454/1‑93 (MQ=255)

                                                                           gATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTg              >  1:5552602/1‑93 (MQ=255)

                                                                             tAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCa            >  2:7247944/1‑93 (MQ=255)

                                                                             tAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCa            >  1:1775361/1‑93 (MQ=255)

                                                                             tAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCa            >  1:3241240/1‑93 (MQ=255)

                                                                             tAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCa            <  2:2558420/93‑1 (MQ=255)

                                                                              aaGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCaa           >  2:438810/1‑93 (MQ=255)

                                                                                gggTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCaaaa         >  2:3423182/1‑93 (MQ=255)

                                                                                gggTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCaaaa         <  1:4492475/93‑1 (MQ=255)

                                                                                  gTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCaaaaaa       <  2:3299762/93‑1 (MQ=255)

                                                                                       cgGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATaaaa  <  2:4155705/93‑1 (MQ=255)

                                                                                       cgGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATaaaa  <  1:549621/93‑1 (MQ=255)

CTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAA  >  NC_000913/3683001‑3683180

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