BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,168,777	A→G	L274P (CTG→CCG)	ycfS ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,168,777	0	A	G	99.3%	471.4 / ‑4.4	146	L274P (CTG→CCG)	ycfS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (81/64); total (81/65)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

Only 100 of 153 total aligned reads displayed.

ACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATC  >  NC_000913/1168689‑1168866

aCGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGcggcag                                                                                       <  2:4235677/93‑1 (MQ=255)

   gACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTg                                                                                    >  1:1991371/1‑93 (MQ=255)

    aCATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGc                                                                                   >  2:2185750/1‑93 (MQ=255)

    aCATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGc                                                                                   <  2:340467/93‑1 (MQ=255)

     catccatcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCg                                                                                  <  2:2155906/93‑1 (MQ=255)

     catccatcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCg                                                                                  >  2:4551351/1‑93 (MQ=255)

       tccatcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCggc                                                                                >  1:3655462/1‑92 (MQ=255)

       tccatcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGa                                                                                <  1:1071018/93‑1 (MQ=255)

       tccatcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGa                                                                                >  1:3655461/1‑93 (MQ=255)

       tccatcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGa                                                                                >  1:2384798/1‑93 (MQ=255)

         catcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATc                                                                              <  1:135975/93‑1 (MQ=255)

         catcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATc                                                                              >  1:402510/1‑93 (MQ=255)

         catcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATc                                                                              <  1:1302407/93‑1 (MQ=255)

         catcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATc                                                                              >  2:4320334/1‑93 (MQ=255)

          atcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCg                                                                             >  1:3496632/1‑93 (MQ=255)

           tcACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGt                                                                            <  1:3112324/93‑1 (MQ=255)

            cacaTGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGtc                                                                           <  2:2794296/93‑1 (MQ=255)

             acaTGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGtca                                                                          >  1:3814308/1‑93 (MQ=255)

             acaTGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGtca                                                                          <  2:4600497/93‑1 (MQ=255)

               aTGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGtcatc                                                                        <  2:2063385/93‑1 (MQ=255)

                tGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGtcatca                                                                       >  1:4571580/1‑93 (MQ=255)

                tGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGtcatca                                                                       <  2:2788716/93‑1 (MQ=255)

                tGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGtcatca                                                                       >  2:172014/1‑93 (MQ=255)

                tGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGtcatca                                                                       >  1:2551225/1‑93 (MQ=255)

                 gTTGCATCACTTCAGCGTCTGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCaa                                                                      >  2:2163010/1‑93 (MQ=255)

                 gTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCaa                                                                      <  1:1901421/93‑1 (MQ=255)

                 gTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCaa                                                                      >  1:2850381/1‑93 (MQ=255)

                  ttGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAAt                                                                     >  1:268856/1‑93 (MQ=255)

                    gCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAAtct                                                                   <  1:2035716/93‑1 (MQ=255)

                     cATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAAtctt                                                                  <  1:3227034/93‑1 (MQ=255)

                     cATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAAtctt                                                                  >  1:2430307/1‑93 (MQ=255)

                     cATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAAtctt                                                                  >  2:4729190/1‑93 (MQ=255)

                      aTCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAAtcttc                                                                 >  2:3198330/1‑93 (MQ=255)

                      aTCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAAtcttc                                                                 <  1:2092532/93‑1 (MQ=255)

                        cACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTtctc                                                               <  2:2481976/93‑1 (MQ=255)

                        cACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTtctc                                                               <  1:4155877/93‑1 (MQ=255)

                        cACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTtctc                                                               <  2:2803811/93‑1 (MQ=255)

                         aCTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTctct                                                              >  2:3194080/1‑93 (MQ=255)

                         aCTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTctct                                                              <  2:667883/93‑1 (MQ=255)

                          cTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCCATCTTCTCTg                                                             >  1:2756950/1‑93 (MQ=255)

                          cTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTg                                                             >  1:1255002/1‑93 (MQ=255)

                           ttCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGa                                                            >  1:1890210/1‑93 (MQ=255)

                           ttCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGa                                                            >  1:4707732/1‑93 (MQ=255)

                           ttCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGa                                                            >  1:4160733/1‑93 (MQ=255)

                            tCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGAc                                                           >  2:1761753/1‑93 (MQ=255)

                              aGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAg                                                         <  1:804681/93‑1 (MQ=255)

                               gCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGc                                                        <  1:2258008/93‑1 (MQ=255)

                                 gTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCgg                                                      >  2:3167394/1‑93 (MQ=255)

                                 gTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCgg                                                      >  2:1338267/1‑93 (MQ=255)

                                  tCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGc                                                     >  2:3435693/1‑93 (MQ=255)

                                  tCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGc                                                     >  1:4414777/1‑93 (MQ=255)

                                  tCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGc                                                     >  1:3036086/1‑93 (MQ=255)

                                   cAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCt                                                    >  2:1915257/1‑93 (MQ=255)

                                       ttGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATg                                                <  2:3582369/93‑1 (MQ=255)

                                        tgtgCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGt                                               >  1:713507/1‑93 (MQ=255)

                                        tgtgCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGt                                               <  1:307188/93‑1 (MQ=255)

                                             tgcATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACtta                                          >  1:4050619/1‑92 (MQ=255)

                                             tgcATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTc                                          <  2:3933291/93‑1 (MQ=255)

                                             tgcATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTc                                          >  1:291519/1‑93 (MQ=255)

                                                aTCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAAc                                       >  1:1752073/1‑93 (MQ=255)

                                                aTCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAAc                                       >  1:4288542/1‑93 (MQ=255)

                                                aTCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAAc                                       <  1:2235140/93‑1 (MQ=255)

                                                    ttAAATGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGGTGTACTTCAACCAGa                                   >  2:4468021/1‑93 (MQ=255)

                                                      aaaTGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGAcg                                 >  2:4419335/1‑93 (MQ=255)

                                                       aaTGATTGCATTGCGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGAcgc                                <  1:3548867/93‑1 (MQ=255)

                                                          gATTGCATTGCGCTAGTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAAACGTCACTGACAGCGGCTGAAGTACTACAACCAGACGCGCa                             >  1:4201824/1‑93 (MQ=255)

                                                           attgcattgcGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCAc                            <  2:2624674/93‑1 (MQ=255)

                                                            ttgcattgcGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCAcc                           >  2:241938/1‑93 (MQ=255)

                                                              gcattgcGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTAGTTCAACCAGACGCGGACCGt                         >  1:3380407/1‑93 (MQ=255)

                                                               cattgcGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGtt                        <  2:4045976/93‑1 (MQ=255)

                                                               cattgcGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGtt                        <  1:1967738/93‑1 (MQ=255)

                                                                attgcGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGttt                       >  2:2779158/1‑93 (MQ=255)

                                                                attgcGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGttt                       <  1:776676/93‑1 (MQ=255)

                                                                 ttgcGCTATTCAGCGTAATTGGCGGCAGCTGCGGACCGTCATCATTCTTCTCCGACAGCGGCTGATGTACTTCAACCAGACGCGCCCCGTTTg                      <  2:3574720/93‑1 (MQ=255)

                                                                  tgcGCTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTgg                     <  1:1531568/93‑1 (MQ=255)

                                                                   gcgcTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGt                    <  1:3704988/93‑1 (MQ=255)

                                                                    cgcTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGtt                   <  2:417296/93‑1 (MQ=255)

                                                                     gcTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTc                  <  2:3174376/93‑1 (MQ=255)

                                                                     gcTATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTc                  >  2:4721719/1‑93 (MQ=255)

                                                                       tATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCgg                >  1:2493350/1‑93 (MQ=255)

                                                                       tATTCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCgg                >  2:1596547/1‑93 (MQ=255)

                                                                         ttCAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCa              <  2:3020689/93‑1 (MQ=255)

                                                                           cAGCGTAATTGGCGGCAGCTGCGGATCGTCATCACTCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCaga            <  1:4207304/93‑1 (MQ=255)

                                                                           cAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCaga            <  2:3699490/93‑1 (MQ=255)

                                                                           cAGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCaga            <  2:3215436/93‑1 (MQ=255)

                                                                           cAGCGTAATTGGCGGCAGCTGCGAATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCaga            <  2:3215437/93‑1 (MQ=255)

                                                                            aGCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCagag           >  1:4338876/1‑93 (MQ=255)

                                                                             gCGTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAgaga          <  2:4213952/93‑1 (MQ=255)

                                                                               gTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACt        >  1:1919169/1‑93 (MQ=255)

                                                                               gTAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACt        <  1:4451002/93‑1 (MQ=255)

                                                                                tAATTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACtt       >  2:4318287/1‑93 (MQ=255)

                                                                                 aaTTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACttt      <  1:1028390/93‑1 (MQ=255)

                                                                                 aaTTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACttt      <  2:3927328/93‑1 (MQ=255)

                                                                                 aaTTGGCCGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACttt      <  1:976735/93‑1 (MQ=255)

                                                                                  aTTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGAGttgt     >  1:175729/1‑91 (MQ=255)

                                                                                  aTTGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACtttt     >  2:1912159/1‑93 (MQ=255)

                                                                                   ttGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTa    <  1:1451765/93‑1 (MQ=255)

                                                                                   ttGGCGGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTa    <  2:2917909/93‑1 (MQ=255)

                                                                                    tGGCGGCAGCTGCGGACCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTAt   >  1:1380601/1‑93 (MQ=255)

                                                                                     ggcggcAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATc  <  1:4643075/93‑1 (MQ=255)

ACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCATTGCGCTATTCAGCGTAATTGGCAGCAGCTGCGGATCGTCATCAATCTTCTCTGACAGCGGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATC  >  NC_000913/1168689‑1168866

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