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.4%	543.9 / ‑4.7	169	V42A (GTC→GCC)	yhjK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (105/63); total (106/63)
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.17e-01

Only 100 of 170 total aligned reads displayed.

CTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAA  >  NC_000913/3683001‑3683180

  gACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCaga                                                                                       >  1:1583212/1‑93 (MQ=255)

   aCCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCagag                                                                                      <  1:168736/93‑1 (MQ=255)

    ccACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagc                                                                                     <  1:663139/93‑1 (MQ=255)

      aCGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcga                                                                                   <  1:1259950/93‑1 (MQ=255)

       cGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcgag                                                                                  >  2:507198/1‑93 (MQ=255)

       cGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcgag                                                                                  >  2:507197/1‑93 (MQ=255)

       cGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAgagcgag                                                                                  <  1:2956274/93‑1 (MQ=255)

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

          cTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCaa                                                                               <  2:799020/93‑1 (MQ=255)

          cTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCaa                                                                               <  2:3662272/93‑1 (MQ=255)

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

             ccGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGc                                                                            >  1:2330021/1‑93 (MQ=255)

               ggTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCtt                                                                          >  2:1167617/1‑93 (MQ=255)

               ggTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCtt                                                                          >  1:4442578/1‑93 (MQ=255)

               cgTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCtt                                                                          <  2:1209797/92‑1 (MQ=255)

                gTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCttt                                                                         <  1:2693655/93‑1 (MQ=255)

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

                gTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCttt                                                                         <  2:3276090/93‑1 (MQ=255)

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

                 tttAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTc                                                                        <  2:3433377/93‑1 (MQ=255)

                  ttAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTcc                                                                       >  1:3651107/1‑93 (MQ=255)

                  ttAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTcc                                                                       >  1:1107319/1‑93 (MQ=255)

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

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

                       gCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTg                                                                  <  1:1182967/93‑1 (MQ=255)

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

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

                             aaGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGc                                                            <  2:174104/93‑1 (MQ=255)

                             aaGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGc                                                            >  2:2135677/1‑93 (MQ=255)

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

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

                              aGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGcc                                                           <  2:2297291/93‑1 (MQ=255)

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

                              aGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGcc                                                           <  2:3852616/93‑1 (MQ=255)

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

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

                                 aTAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCgtg                                                        <  2:4470099/93‑1 (MQ=255)

                                  tAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCgtgt                                                       <  1:1450844/93‑1 (MQ=255)

                                   aGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGtt                                                      <  2:134052/93‑1 (MQ=255)

                                   aGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGtt                                                      >  2:642070/1‑93 (MQ=255)

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

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

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

                                       tCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAg                                                  <  1:928529/93‑1 (MQ=255)

                                         cGCTTCCAGAATATCTGCTTGCAAAACAGCTGAATATAATGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGGGTTATAGCg                                                <  2:3510636/93‑1 (MQ=255)

                                          gCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGa                                               <  1:12160/93‑1 (MQ=255)

                                            ttagggAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAtt                                             >  2:2909221/4‑93 (MQ=255)

                                            ttCGGGAATATCGCCTTTCAAAATGGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAtt                                             >  2:274214/1‑93 (MQ=255)

                                            ttCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAtt                                             <  1:751587/93‑1 (MQ=255)

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

                                              cGGGAATATCGCCTTTCAAAATAGCTGAAGATACGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTct                                           <  1:4575097/93‑1 (MQ=255)

                                              cGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTct                                           <  1:2065636/93‑1 (MQ=255)

                                              cGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTct                                           <  2:273322/93‑1 (MQ=255)

                                              cGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTct                                           <  2:3102078/93‑1 (MQ=255)

                                                ggAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTctgc                                         >  1:2635129/1‑93 (MQ=255)

                                                 gAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTctgct                                        >  2:2118284/1‑93 (MQ=255)

                                                 gAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTctgct                                        <  2:206009/93‑1 (MQ=255)

                                                   atatCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTATTATAGCGATTCTCCTgg                                      >  1:571819/1‑93 (MQ=255)

                                                      tCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGAcc                                   >  2:2006181/1‑93 (MQ=255)

                                                      tCGCCTTTCAAAATAGCTGAAGACAAGGGTTCGCGGGCAGAGCGAGCAATCCTTTCCGGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGAca                                   >  2:1667200/1‑92 (MQ=255)

                                                        gCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAg                                 >  2:4048663/1‑93 (MQ=255)

                                                        gCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAg                                 >  2:1910016/1‑93 (MQ=255)

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

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

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

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

                                                           tttCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATg                              >  2:2075777/1‑93 (MQ=255)

                                                           tttCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATg                              >  2:1310399/1‑93 (MQ=255)

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

                                                             tCAAAATAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCATAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGa                            <  2:4645351/93‑1 (MQ=255)

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

                                                               aaaaTAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAAc                          <  1:3468337/93‑1 (MQ=255)

                                                               aaaaTAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAAc                          <  1:1916447/93‑1 (MQ=255)

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

                                                                aaaTAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTATGCTGGACCAGATGGAACa                         >  2:4185666/1‑93 (MQ=255)

                                                                aaaTAGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAAAGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACa                         >  1:2459211/1‑93 (MQ=255)

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

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

                                                                    aGCTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCaa                     >  1:3191012/1‑93 (MQ=255)

                                                                      cTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCaaaa                   >  2:186576/1‑93 (MQ=255)

                                                                      cTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCaaaa                   <  1:402736/93‑1 (MQ=255)

                                                                      cTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCaaaa                   >  2:1100481/1‑93 (MQ=255)

                                                                      cTGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCaaaa                   <  1:1758064/93‑1 (MQ=255)

                                                                       tGAAGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAAc                  <  1:2213551/93‑1 (MQ=255)

                                                                          aGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGt               >  1:1473541/1‑93 (MQ=255)

                                                                          aGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGt               <  1:1606296/93‑1 (MQ=255)

                                                                          aGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGt               >  2:2057743/1‑93 (MQ=255)

                                                                          aGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGt               <  2:3746445/93‑1 (MQ=255)

                                                                          aGATAAGGGTTCGCGGGCAGAGCGAGCAATGCTGTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGt               <  1:4578636/93‑1 (MQ=255)

                                                                           gATAAGGGTTCGCGGGCAGAGCGTGCAATGCTTTCCAGTTGCGTTGCCGTGCTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTg              >  2:3790856/1‑93 (MQ=255)

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

                                                                            aTAAGGGTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGc             >  2:1407999/1‑93 (MQ=255)

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

                                                                                  gTTCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCaaaaaa       >  1:4024430/1‑93 (MQ=255)

                                                                                   ttCGCGGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTATGCTCGACCAGATGGAACAGCAAAACGGAGCAAAAAAt      >  2:3285846/1‑93 (MQ=255)

                                                                                   ttCGCGAGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAAt      <  1:3679312/93‑1 (MQ=255)

                                                                                      gcgGGCAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATaaa   <  1:2366998/93‑1 (MQ=255)

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

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

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

CTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAA  >  NC_000913/3683001‑3683180

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