BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,206,988	A→G	T89A (ACT→GCT)	ycfK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,206,988	0	A	G	98.4%	384.2 / ‑2.9	128	T89A (ACT→GCT)	ycfK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (59/67); total (61/67)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 130 total aligned reads displayed.

GAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1206900‑1207077

gaAGGGGCGGCACGACCAGTTACTTCCCGCACTTCGCGCGCTGCTGTAAAGCGGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCtgt                                                                                       <  2:2808058/93‑1 (MQ=18)

 aaGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCtgtg                                                                                      <  2:1298422/93‑1 (MQ=255)

  aGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGc                                                                                     >  2:3120547/1‑93 (MQ=255)

  aGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGc                                                                                     >  1:925474/1‑93 (MQ=255)

   ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCa                                                                                    >  2:1845399/1‑93 (MQ=255)

   ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCa                                                                                    >  1:266598/1‑93 (MQ=255)

     ggcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCaaa                                                                                  <  2:3750495/93‑1 (MQ=255)

     ggcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCaaa                                                                                  <  2:900338/93‑1 (MQ=255)

     ggcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCaaa                                                                                  >  2:84423/1‑93 (MQ=255)

     ggcggcACGACCAGTTACGTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCaaa                                                                                  <  2:2620434/93‑1 (MQ=255)

      gcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAAc                                                                                 >  2:1303259/1‑93 (MQ=255)

        ggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACgg                                                                               >  1:822338/1‑93 (MQ=255)

           aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGctc                                                                            >  1:3438661/1‑93 (MQ=255)

           aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGctc                                                                            <  1:4188110/93‑1 (MQ=255)

           aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGctc                                                                            <  2:3695564/93‑1 (MQ=255)

           aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGctc                                                                            <  2:2534592/93‑1 (MQ=255)

            cGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGctct                                                                           >  2:1117787/1‑93 (MQ=255)

             gACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCtctc                                                                          >  1:1740936/1‑93 (MQ=255)

             gACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCtctc                                                                          <  2:2067123/93‑1 (MQ=255)

             gACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCtctc                                                                          <  1:1049023/93‑1 (MQ=255)

               ccAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGa                                                                        >  1:1360163/1‑93 (MQ=255)

               ccAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGa                                                                        >  1:2200699/1‑93 (MQ=255)

                  gTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGaaaa                                                                     >  2:1563379/1‑93 (MQ=255)

                  gTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGaaaa                                                                     >  1:3232939/1‑93 (MQ=255)

                  gTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGaaaa                                                                     >  1:1644674/1‑93 (MQ=255)

                   ttacttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAAc                                                                    >  1:2096911/1‑93 (MQ=255)

                    tacttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGAGGGCGCTGTGCAAACGGCTCTCGAAAAcc                                                                   <  1:4084374/93‑1 (MQ=255)

                     aCTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCT                                                                  <  2:3430385‑M1/93‑1 (MQ=255)

                     aCTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCT                                                                  <  2:2897566‑M1/93‑1 (MQ=255)

                      cTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTT                                                                 <  2:3046902‑M1/93‑1 (MQ=255)

                      cTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTT                                                                 <  1:4188400‑M1/93‑1 (MQ=255)

                      cTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTT                                                                 >  2:3780163‑M1/1‑93 (MQ=255)

                       ttACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTG                                                                <  1:1453889‑M1/93‑1 (MQ=255)

                        tACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGG                                                               <  1:833493‑M1/93‑1 (MQ=255)

                            gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTG                                                           <  1:361350‑M1/93‑1 (MQ=255)

                            gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTG                                                           >  1:1151257‑M1/1‑93 (MQ=255)

                            gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTG                                                           <  1:2537777‑M1/93‑1 (MQ=255)

                              aCTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGG                                                         <  1:1020554‑M1/93‑1 (MQ=255)

                              aCTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGG                                                         >  2:4080228‑M1/1‑93 (MQ=255)

                                ttCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAG                                                       <  2:77936‑M1/93‑1 (MQ=255)

                                 tCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                      <  1:2012958‑M1/93‑1 (MQ=255)

                                  cgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAA                                                     <  2:202190‑M1/93‑1 (MQ=255)

                                  cgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAA                                                     <  1:2566264‑M1/93‑1 (MQ=255)

                                  cgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAA                                                     <  1:3929164‑M1/93‑1 (MQ=255)

                                    cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                   <  1:586281‑M1/93‑1 (MQ=255)

                                      cgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGct                                                 <  2:3048670‑M1/93‑3 (MQ=255)

                                       gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctc                                                <  1:1574492‑M1/93‑4 (MQ=255)

                                       gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                >  2:1075167/1‑93 (MQ=255)

                                        ctgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                               >  2:2953186‑M1/1‑89 (MQ=255)

                                          gctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                             <  1:4199817‑M1/93‑7 (MQ=255)

                                           ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                            <  1:3014744‑M1/93‑8 (MQ=255)

                                            tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcat                                           >  2:830133‑M1/1‑85 (MQ=255)

                                            tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcat                                           >  1:2660079‑M1/1‑85 (MQ=255)

                                            tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                           <  2:1235014/93‑1 (MQ=255)

                                               tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                        <  1:775945‑M1/93‑12 (MQ=255)

                                                aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacc                                       <  1:3886086‑M1/93‑13 (MQ=255)

                                                aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacc                                       <  2:304635‑M1/93‑13 (MQ=255)

                                                aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCaa                                       >  1:3112969/1‑93 (MQ=255)

                                                 aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                      <  2:4450682‑M1/93‑14 (MQ=255)

                                                  gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctg                                     >  2:2931373‑M1/1‑79 (MQ=255)

                                                    cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                   <  2:4543751‑M1/93‑17 (MQ=255)

                                                    cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                   <  2:2777280‑M1/93‑17 (MQ=255)

                                                     gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                  >  2:4142735‑M1/1‑76 (MQ=255)

                                                     gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                  <  2:709877‑M1/93‑18 (MQ=255)

                                                     gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                  >  2:522405/1‑93 (MQ=255)

                                                     gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                  >  1:4275467/1‑93 (MQ=255)

                                                     gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                  >  2:1608100/1‑93 (MQ=255)

                                                     gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                  >  1:539488‑M1/1‑76 (MQ=255)

                                                      cAAGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                 >  1:3009344/1‑93 (MQ=255)

                                                       aaGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                <  1:1423756‑M1/93‑20 (MQ=255)

                                                       aaGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                <  1:2877073‑M1/93‑20 (MQ=255)

                                                        aGAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCaa                               <  1:2803997/93‑1 (MQ=255)

                                                         gAATCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgt                              <  1:4356048‑M1/93‑22 (MQ=255)

                                                           aTCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcc                            >  2:822364‑M1/1‑70 (MQ=255)

                                                            tCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                           <  2:3751211‑M1/93‑25 (MQ=255)

                                                               gTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        <  2:4263601/93‑1 (MQ=255)

                                                               gTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        >  2:690717/1‑93 (MQ=255)

                                                                tttGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggtt                       <  2:2885983‑M1/93‑29 (MQ=255)

                                                                 ttGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                      <  1:3990304/93‑1 (MQ=255)

                                                                  tGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttcc                     <  2:1527254‑M1/93‑31 (MQ=255)

                                                                  tGGCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttcc                     >  2:2807305‑M1/1‑63 (MQ=255)

                                                                   ggCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                    <  2:3937661/93‑1 (MQ=255)

                                                                    gCGATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                   <  1:779980‑M1/93‑33 (MQ=255)

                                                                      gATATCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtgg                 >  1:1194656‑M1/1‑59 (MQ=255)

                                                                       atatCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                >  2:1227773/1‑93 (MQ=255)

                                                                        tatCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcc               >  2:3478122‑M1/1‑57 (MQ=255)

                                                                        tatCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcc               >  2:1015518‑M1/1‑57 (MQ=255)

                                                                         atCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct              >  1:2452138‑M1/1‑56 (MQ=255)

                                                                          tCAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcctt             >  2:95832‑M1/1‑55 (MQ=255)

                                                                           cAAATCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttc            <  1:4519954‑M1/93‑40 (MQ=255)

                                                                            aaaTCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttca           <  1:1713443‑M1/93‑41 (MQ=255)

                                                                             aaTCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcag          >  2:2053215‑M1/1‑52 (MQ=255)

                                                                              aTCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagc         >  1:2036054‑M1/1‑51 (MQ=255)

                                                                               tCGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccntcagcc        >  1:1424098‑M1/1‑50 (MQ=255)

                                                                                cGGATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt       <  2:787880/93‑1 (MQ=255)

                                                                                  gATGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc     >  2:1407187/1‑93 (MQ=255)

                                                                                    tGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactcc   >  2:2097090‑M1/1‑45 (MQ=255)

                                                                                    tGGCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactcc   >  2:1651301‑M1/1‑45 (MQ=255)

                                                                                     ggCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccg  >  2:2187813‑M1/1‑44 (MQ=255)

                                                                                     ggCGCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa  <  1:1471583/93‑1 (MQ=255)

GAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1206900‑1207077

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