BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,069,008	C→T	G110D (GGC→GAC)	yihS ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,069,008	0	C	T	99.6%	781.3 / ‑5.8	228	G110D (GGC→GAC)	yihS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (116/110); total (117/111)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 249 total aligned reads displayed.

TTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCT  >  NC_000913/4068916‑4069097

tttCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGATGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGc                                                                                           >  1:4614269/1‑93 (MQ=255)

  tCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGtcc                                                                                         >  1:3509758/1‑93 (MQ=255)

  tCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGtcc                                                                                         <  1:188582/93‑1 (MQ=255)

   cAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCt                                                                                        <  1:397302/93‑1 (MQ=255)

     aTGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGt                                                                                      >  1:3278188/1‑93 (MQ=255)

      tGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCCAAGAAGTGTTGATAGTCCTGtt                                                                                     >  2:124647/1‑93 (MQ=255)

       ggTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGttt                                                                                    >  1:2025301/1‑93 (MQ=255)

         tgtAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTgg                                                                                  >  2:4270385/1‑93 (MQ=255)

           tAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAg                                                                                <  1:2272448/93‑1 (MQ=255)

           tAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAg                                                                                <  2:3732945/93‑1 (MQ=255)

            aaTCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAgg                                                                               >  1:588236/1‑93 (MQ=255)

            aaTCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAgg                                                                               >  2:2815420/1‑93 (MQ=255)

             aTCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGc                                                                              >  1:2082828/1‑93 (MQ=255)

             aTCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGc                                                                              <  1:319099/93‑1 (MQ=255)

             aTCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGc                                                                              >  2:3523283/1‑93 (MQ=255)

             aTCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGc                                                                              <  2:112059/93‑1 (MQ=255)

             aTCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGc                                                                              <  2:533554/93‑1 (MQ=255)

              tCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCa                                                                             <  2:4227442/93‑1 (MQ=255)

              tCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCa                                                                             >  1:4606420/1‑93 (MQ=255)

               cGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGGGTTGATAGTCCTGTTTGGAGGCAt                                                                            >  2:2931880/1‑93 (MQ=255)

                 agcagcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGGGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATcc                                                                          >  1:2008756/1‑93 (MQ=255)

                 agcagcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATcc                                                                          <  1:1390703/93‑1 (MQ=255)

                    agcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATccacc                                                                       <  1:2619429/93‑1 (MQ=255)

                     gcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATccaccc                                                                      >  1:1044277/1‑92 (MQ=255)

                     gcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATccacca                                                                      <  1:912059/93‑1 (MQ=255)

                     gcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATccacca                                                                      >  1:3566127/1‑93 (MQ=255)

                     gcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATccacca                                                                      >  1:1319211/1‑93 (MQ=255)

                     gcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATccacca                                                                      >  2:3277906/1‑93 (MQ=255)

                       ttGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACg                                                                    >  1:2533096/1‑93 (MQ=255)

                       ttGCGCGCTTCCGGGTGACCCGTTGTGAAGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGGTTAGAGGAATCAACCACg                                                                    >  2:2014629/1‑93 (MQ=255)

                        tGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGc                                                                   <  2:2521837/93‑1 (MQ=255)

                           ccgcTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCt                                                                <  2:4342900/92‑1 (MQ=255)

                               ttCCGGGTGACCCGTTGTGACGGCGCTGGCGGCATCCAGCAGAGCAAAGAAGCGTTGATAGTCCTGTTTGGATGCATCCACCACGCCCTCGta                                                            >  2:2304365/1‑92 (MQ=255)

                               ttCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTc                                                            <  2:1641949/93‑1 (MQ=255)

                                   gggTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTc                                                        <  2:3133416/93‑1 (MQ=255)

                                   gggTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTc                                                        <  2:2885974/93‑1 (MQ=255)

                                    ggTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCa                                                       <  1:372088/93‑1 (MQ=255)

                                    ggTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCa                                                       >  2:4403574/1‑93 (MQ=255)

                                     gTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCAc                                                      >  2:1451942/1‑93 (MQ=255)

                                     gTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCAc                                                      <  1:521133/93‑1 (MQ=255)

                                      tGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACg                                                     >  1:1000168/1‑93 (MQ=255)

                                       gACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGc                                                    <  1:4292936/93‑1 (MQ=255)

                                        aCCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCa                                                   <  1:282167/93‑1 (MQ=255)

                                        aCCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGGCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCa                                                   >  1:2213211/1‑93 (MQ=255)

                                          ccGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAgg                                                 <  2:3962932/93‑1 (MQ=255)

                                            gTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCa                                               >  2:3031047/1‑93 (MQ=255)

                                            gTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCa                                               <  1:930280/93‑1 (MQ=255)

                                             ttGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCAt                                              >  1:4152521/1‑93 (MQ=255)

                                             ttGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCAt                                              >  2:700310/1‑93 (MQ=255)

                                             gtgtgACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCAt                                              <  1:1677088/92‑1 (MQ=255)

                                                tgACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATAcc                                           <  1:1445809/93‑1 (MQ=255)

                                                 gACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCa                                          <  2:3477825/93‑1 (MQ=255)

                                                 gACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCa                                          >  2:4131470/1‑93 (MQ=255)

                                                 gACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCa                                          <  2:1872224/93‑1 (MQ=255)

                                                  aCGGCGCTGGCGGCCCCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAg                                         <  2:2277603/93‑1 (MQ=255)

                                                       gcTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAgccgcc                                    <  1:1142898/93‑1 (MQ=255)

                                                        cTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCa                                   <  1:4474771/93‑1 (MQ=255)

                                                        cTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCa                                   >  1:692336/1‑93 (MQ=255)

                                                          ggcggcACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCata                                 <  2:1579698/93‑1 (MQ=255)

                                                          ggcggcACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCata                                 <  2:347046/93‑1 (MQ=255)

                                                          ggcggcACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCata                                 <  2:843165/93‑1 (MQ=255)

                                                          ggcggcACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCata                                 >  1:313636/1‑93 (MQ=255)

                                                          ggcggcACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCata                                 >  1:3768472/1‑93 (MQ=255)

                                                           gcggcACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCatat                                <  1:3081075/93‑1 (MQ=255)

                                                           gcggcACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCatat                                <  1:845943/93‑1 (MQ=255)

                                                           gcggcACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCatat                                <  2:2194328/93‑1 (MQ=255)

                                                            cggcCCCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATAtt                               <  2:4547547/93‑1 (MQ=255)

                                                              gcACCCAGCAGAGCAAATAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATAtttt                             <  2:272411/93‑1 (MQ=255)

                                                              gcACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATAtttt                             <  2:2708254/93‑1 (MQ=255)

                                                               cACCCAGCAGAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATAttttt                            >  2:759455/1‑93 (MQ=255)

                                                                 cccAGCAGAGCAAAGAGGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTa                          <  1:2373127/93‑1 (MQ=255)

                                                                   cagcagAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATc                        <  1:3623958/93‑1 (MQ=255)

                                                                   cagcagAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATc                        <  1:3922748/93‑1 (MQ=255)

                                                                    agcagAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCa                       <  1:1158776/93‑1 (MQ=255)

                                                                     gcagAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCAc                      >  2:2047784/1‑93 (MQ=255)

                                                                     gcagAGCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCAc                      >  1:791366/1‑93 (MQ=255)

                                                                      cagAGCAAAGAAGTGTTGATAGTCCTGTTGGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACg                     <  2:3240499/93‑1 (MQ=255)

                                                                       agagCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGc                    >  1:2898765/1‑93 (MQ=255)

                                                                       agagCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGc                    <  1:2742096/93‑1 (MQ=255)

                                                                        gagCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCa                   <  2:2379244/93‑1 (MQ=255)

                                                                          gCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGt                 >  2:898470/1‑93 (MQ=255)

                                                                          gCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGt                 <  1:4214050/93‑1 (MQ=255)

                                                                          gCAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTAGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGt                 <  2:3820893/93‑1 (MQ=255)

                                                                           cAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTg                <  1:2132636/93‑1 (MQ=255)

                                                                           cAAAGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTAGCACGCAGTg                >  2:3632514/1‑93 (MQ=255)

                                                                            aaaGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTgc               <  1:4376450/93‑1 (MQ=255)

                                                                            aaaGAAGTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTgc               <  2:2337285/93‑1 (MQ=255)

                                                                             aagaagTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTgcg              <  2:626682/93‑1 (MQ=255)

                                                                             aagaagTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTgcg              >  1:352951/1‑93 (MQ=255)

                                                                               gaagTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGcc            >  2:3756724/1‑93 (MQ=255)

                                                                               gaagTGTTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGGCATTCACGCAGGCATAGCAGCCGCCATATTTTTTATCACGCAGTGCGcc            >  2:3660662/1‑93 (MQ=255)

                                                                                  gtgtTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGtt         >  2:3564533/1‑93 (MQ=255)

                                                                                  gtgtTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGtt         >  1:1877624/1‑93 (MQ=255)

                                                                                    gtTGATAGTCCTGTTTGGAGGCATCCACCTCGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCa       <  2:383220/93‑1 (MQ=255)

                                                                                    gtTGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCa       >  2:2944447/1‑93 (MQ=255)

                                                                                      tGATAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATg     >  1:565852/1‑93 (MQ=255)

                                                                                        aTAGTCCTGTTTGGAGGCATCCACCACGCCTTCGTCATTCACGCAGCCATCCCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGc   <  2:4594580/93‑1 (MQ=255)

                                                                                        aTAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATCCCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGc   <  2:1869169/93‑1 (MQ=255)

                                                                                        aTAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGc   <  1:1504890/93‑1 (MQ=255)

                                                                                        aTAGTCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGc   <  1:4534628/93‑1 (MQ=255)

TTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCT  >  NC_000913/4068916‑4069097

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