BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,208,331	A→G	S404S (TCA→TCG)	b2122 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,208,331	0	A	G	99.1%	327.8 / ‑4.4	105	S404S (TCA→TCG)	b2122	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (30/74); total (30/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 128 total aligned reads displayed.

GCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTG  >  NC_000913/2208241‑2208420

gCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTcgca                                                                                         <  2:2545722/93‑1 (MQ=255)

 ccGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCaa                                                                                        <  1:1686054/93‑1 (MQ=255)

  cgcgAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAag                                                                                       <  2:4733474/93‑1 (MQ=255)

   gcgAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAaga                                                                                      >  2:1662575/1‑92 (MQ=255)

    cgAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAagca                                                                                     <  2:323795/93‑1 (MQ=255)

    cgAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAagca                                                                                     <  1:2775173/93‑1 (MQ=255)

     gAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGGCAATTGTTTTGCcagcag                                                                                    >  2:227762/1‑93 (MQ=255)

     gAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCCAAATTTTTCGCaaaccg                                                                                    >  2:2123450/1‑89 (MQ=255)

     gAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGCGCAAATTGTTTCGCAaacag                                                                                    >  1:812963/1‑93 (MQ=255)

      aaGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAagcagc                                                                                   <  2:2710134/93‑1 (MQ=255)

        gTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCCAATTGGTTTGCAAGGAGCaa                                                                                 >  2:2203228/1‑93 (MQ=255)

        gTGAAAAGGAGAACTCGCCGGCAGATGCGCGTAGGCTGTTTAGTTTTCAGGCTGGCGCTCTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCaa                                                                                 <  2:496057/93‑1 (MQ=255)

        atGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCaa                                                                                 <  1:1380932/92‑1 (MQ=255)

           aaCAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGAGCACTGGCGGGCGGGCAAATTGCTTGGCAAGGAGCAAAAc                                                                              >  1:1860541/1‑93 (MQ=255)

            aCAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACg                                                                             >  2:383052/1‑93 (MQ=255)

              aGGAAACCGCACCTGCAGCGACGCGTCGGCGGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTGCCGCAAGCAGCAAAACgaa                                                                           <  2:3871744/93‑3 (MQ=255)

               ggAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGGCAATTGTTTCGCAAGCAGCAAAACGCAg                                                                          >  1:4329369/1‑93 (MQ=255)

               ggAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAg                                                                          <  2:3882838/93‑1 (MQ=255)

                gAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGt                                                                         <  1:3839263/93‑1 (MQ=255)

                gAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGt                                                                         <  1:3143576/93‑1 (MQ=255)

                gAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTCTCACGCTGGCGCACTGGCGGGCGGGCGAACTGTTTCTGCAGGCGCAAAACGCacg                                                                         >  1:2092184/1‑91 (MQ=255)

                 aaaaCTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTg                                                                        <  2:803726/93‑1 (MQ=255)

                 aaaaCTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTg                                                                        <  2:803725/93‑1 (MQ=255)

                 aaaaCTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTg                                                                        <  1:3276956/93‑1 (MQ=255)

                   aaCTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGcc                                                                      <  2:4685889/93‑1 (MQ=255)

                   aaCTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGcc                                                                      <  1:320995/93‑1 (MQ=255)

                   aaCTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAACAAAACGCAGTGcc                                                                      >  2:2298319/1‑93 (MQ=255)

                      tCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTTGCCAGCCACAAAACGCAGTGCCGAt                                                                   >  2:2015518/1‑93 (MQ=255)

                      tCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGAt                                                                   <  2:3560684/93‑1 (MQ=255)

                      tCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGAt                                                                   <  1:3570560/93‑1 (MQ=255)

                       cACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCCAGCAGCAAAAAGCAGTGCCGATg                                                                  >  1:1268676/1‑93 (MQ=255)

                       cACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATg                                                                  <  1:156968/93‑1 (MQ=255)

                         ccGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGc                                                                >  2:1375445/1‑93 (MQ=255)

                         ccGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGc                                                                <  1:1948796/93‑1 (MQ=255)

                           ggCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTTGCAAGCAGCAAAACGCAGTGCCGATGGCgg                                                              >  1:2566742/1‑92 (MQ=255)

                           ggCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGa                                                              <  2:545257/93‑1 (MQ=255)

                            gCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAg                                                             <  1:597015/93‑1 (MQ=255)

                            gCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGGCGTGCCGATGGCGAg                                                             >  2:3227303/1‑93 (MQ=255)

                             cAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAgc                                                            <  2:3267971/93‑1 (MQ=255)

                              aGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAgct                                                           <  1:3249557/93‑1 (MQ=255)

                               gCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAgctg                                                          >  2:2562985/1‑93 (MQ=255)

                               gCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAgctg                                                          >  2:1622297/1‑93 (MQ=255)

                                    gcgTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGtt                                                     <  1:2038134/93‑1 (MQ=255)

                                    gcgTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGtt                                                     >  2:1907107/1‑93 (MQ=255)

                                      gTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAg                                                   >  1:3540352/1‑93 (MQ=255)

                                      gTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAACAGCAAAACGCAGTGCCGATCGCGAGCTGCTGTTAg                                                   >  2:760700/1‑93 (MQ=255)

                                       tCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGc                                                  <  1:3025961/93‑1 (MQ=255)

                                       tCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGc                                                  <  2:1188281/93‑1 (MQ=255)

                                        tggCTTTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCt                                                 <  2:503811/92‑1 (MQ=255)

                                        cGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCt                                                 >  1:1879235/1‑93 (MQ=255)

                                         ggCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTa                                                >  1:2974894/1‑93 (MQ=255)

                                           cTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAGATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAcc                                              >  2:1572475/1‑93 (MQ=255)

                                           cTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAcc                                              <  1:1890300/93‑1 (MQ=255)

                                           cTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAcc                                              <  1:2966041/93‑1 (MQ=255)

                                           cTGTTTAGTTTTCAGGCCGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAcc                                              <  2:3124796/93‑1 (MQ=255)

                                            tGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTTTTAGCTActc                                             >  1:2296241/1‑91 (MQ=255)

                                            tGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTAccc                                             <  1:1264757/93‑1 (MQ=255)

                                              tttAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGc                                           >  1:4677880/1‑93 (MQ=255)

                                                tAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCaa                                         <  1:1044258/93‑1 (MQ=255)

                                                    gttCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGc                                     <  2:2168559/92‑1 (MQ=255)

                                                     ttCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCt                                    <  1:4422786/93‑1 (MQ=255)

                                                     ttCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCt                                    <  1:4383110/93‑1 (MQ=255)

                                                        aGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTaa                                 >  1:1048430/1‑93 (MQ=255)

                                                        aGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTaa                                 <  1:938778/93‑1 (MQ=255)

                                                          gCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAagc                               <  2:1678582/93‑1 (MQ=255)

                                                          gCTGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAagc                               <  1:789848/93‑1 (MQ=255)

                                                            tGGCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAagcag                             <  2:3675756/93‑1 (MQ=255)

                                                             ggCGCACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAagcagc                            <  1:3256042/93‑1 (MQ=255)

                                                               cgcACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGt                          <  1:4051306/93‑1 (MQ=255)

                                                                 cACTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTtg                        >  1:1073037/1‑93 (MQ=255)

                                                                  aCTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTtgt                       >  1:1769961/1‑93 (MQ=255)

                                                                   cTGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTtgtg                      >  1:2744514/1‑93 (MQ=255)

                                                                    tGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGc                     >  1:1481373/1‑93 (MQ=255)

                                                                    tGGCGGGCGGGCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGc                     <  1:3754357/93‑1 (MQ=255)

                                                                       cgggcgggCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCGGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGc                  <  2:2827124/93‑1 (MQ=255)

                                                                       cgggcgggCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGc                  >  2:2336054/1‑93 (MQ=255)

                                                                       cgggcgggCAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGc                  <  2:2369379/93‑1 (MQ=255)

                                                                        gggcgggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAATCGCTTAAGCAGCGTTGTGCCGCt                 >  2:3241608/1‑93 (MQ=255)

                                                                        gggcgggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCt                 <  2:2675504/93‑1 (MQ=255)

                                                                         ggcgggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTg                >  2:222421/1‑93 (MQ=255)

                                                                          gcgggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGt               <  1:158647/93‑1 (MQ=255)

                                                                          gcgggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGt               <  2:3100328/93‑1 (MQ=255)

                                                                          gcgggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGt               <  2:4033421/93‑1 (MQ=255)

                                                                           cgggcAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTc              <  1:3816851/93‑1 (MQ=255)

                                                                            gggcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCg             >  2:1337505/1‑93 (MQ=255)

                                                                              gcAAATTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGcc           <  1:1111832/93‑1 (MQ=255)

                                                                                aaaTTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTg         >  1:1783191/1‑93 (MQ=255)

                                                                                aaaTTGTTTCGAAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTg         <  1:4303827/93‑1 (MQ=255)

                                                                                 aaTTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGa        >  2:2719767/1‑93 (MQ=255)

                                                                                 aaTTGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGa        >  1:1605920/1‑93 (MQ=255)

                                                                                  aTTGTTTCGCAAGCAGCAAAACGCAGTGCGGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGAt       <  2:1971876/93‑1 (MQ=255)

                                                                                   ttGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATg      >  2:1185448/1‑93 (MQ=255)

                                                                                   ttGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATg      >  2:3054362/1‑93 (MQ=255)

                                                                                   ttGTTTCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATg      <  1:2639220/93‑1 (MQ=255)

                                                                                     gTTTCGCAAGCAGCACACCGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGcc    <  2:4204615/93‑1 (MQ=255)

                                                                                      tttCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCTTTGTGCCGCTGTCGCCTGATGCCt   <  2:4686810/93‑1 (MQ=255)

                                                                                      tttCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCt   <  2:4257414/93‑1 (MQ=255)

                                                                                       ttCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTg  <  1:498685/93‑1 (MQ=255)

                                                                                       ttCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTg  <  1:3573/93‑1 (MQ=255)

                                                                                       ttCGCAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTg  <  2:1921225/93‑1 (MQ=255)

GCCGCGAAGTGAACAGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTG  >  NC_000913/2208241‑2208420

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