BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	518,860	A→G	15.7%	F97S (TTC→TCC)	ybbO ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	518,860	0	A	G	15.7%	397.4 / 60.1	170	F97S (TTC→TCC)	ybbO	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (65/78); new base G (9/18); total (74/96)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 175 total aligned reads displayed.

CAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCA  >  NC_000913/518763‑518959

cAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAAt                                                                                                   >  2:2523435/1‑100 (MQ=255)

  ggCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATcc                                                                                                 >  1:3243696/1‑100 (MQ=255)

  ggCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATcc                                                                                                 <  2:2910382/100‑1 (MQ=255)

    cgcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCgg                                                                                               <  2:437977/100‑1 (MQ=255)

    cgcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCgg                                                                                               <  2:1581612/100‑1 (MQ=255)

     gcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGc                                                                                              <  1:1194754/100‑1 (MQ=255)

      cATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCa                                                                                             >  2:499498/1‑100 (MQ=255)

       aTGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTGCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAAtactg                                                                                               >  1:2555664/1‑93 (MQ=255)

       aTGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCAt                                                                                            <  1:952773/100‑1 (MQ=255)

        tGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtg                                                                                           <  1:3418757/100‑2 (MQ=255)

         ggTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTg                                                                                          >  1:410093/1‑100 (MQ=255)

          gTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgta                                                                                         <  1:100700/100‑4 (MQ=255)

          gTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgta                                                                                         <  1:2142536/100‑4 (MQ=255)

          gTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                         >  1:2462897/1‑100 (MQ=255)

           tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtat                                                                                        <  2:2775938/100‑5 (MQ=255)

           tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtat                                                                                        <  2:3719011/100‑5 (MQ=255)

           tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtat                                                                                        <  2:1969705/100‑5 (MQ=255)

           tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt                                                                                        <  2:1687709/100‑1 (MQ=255)

            gAGCTGGTGTGGGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATACCGGATTCGGCAtgtatg                                                                                       <  1:2396341/100‑6 (MQ=255)

            gAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatg                                                                                       <  2:3464674/100‑6 (MQ=255)

            gAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatg                                                                                       <  2:1959380/100‑6 (MQ=255)

            gAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTa                                                                                       >  1:1352233/1‑100 (MQ=255)

            gAGCTGGTGTCCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGACAGGGGGCCATACATGCCGGATTCGGCAtgtatg                                                                                       <  1:2498946/100‑6 (MQ=255)

             aGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaa                                                                                      <  1:2689238/100‑1 (MQ=255)

              gCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatggc                                                                                     <  1:2955621/100‑8 (MQ=255)

              gCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaaa                                                                                     <  1:642094/100‑1 (MQ=255)

              gCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaaa                                                                                     <  1:44340/100‑1 (MQ=255)

                tGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatggccc                                                                                   <  2:1569172/100‑10 (MQ=255)

                  gtgtGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATc                                                                                 <  2:612961/100‑1 (MQ=255)

                   tgtgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATcc                                                                                <  2:3398089/100‑1 (MQ=255)

                     tgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCGCa                                                                              >  2:1874427/1‑100 (MQ=255)

                        gcCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATAc                                                                           >  1:2456014/1‑100 (MQ=255)

                        gcCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATAc                                                                           <  1:3546131/100‑1 (MQ=255)

                        gcCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATAc                                                                           <  2:2722697/100‑1 (MQ=255)

                         ccGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACa                                                                          <  2:2892609/100‑1 (MQ=255)

                         ccGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACa                                                                          >  1:3629659/1‑100 (MQ=255)

                         ccGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACa                                                                          <  1:1957480/100‑1 (MQ=255)

                            aaaaaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGAc                                                                       <  1:88172/100‑1 (MQ=255)

                               aaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAAt                                                                    <  2:3567417/100‑1 (MQ=255)

                                 gTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTa                                                                  >  1:1009643/1‑100 (MQ=255)

                                   tGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCCGCATTGTTAAAGATCCCATACAGACAATTATc                                                                <  1:2325683/100‑1 (MQ=255)

                                    ggCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTAttg                                                               <  2:1105229/100‑3 (MQ=255)

                                    ggCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTCAAAGCTCCCATACAGCAAATTATCg                                                               >  1:1969523/1‑100 (MQ=255)

                                      cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGt                                                             >  2:3254983/1‑100 (MQ=255)

                                      cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGt                                                             >  1:1808614/1‑100 (MQ=255)

                                        gAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCa                                                           <  1:144195/100‑1 (MQ=255)

                                         aaaaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAg                                                          <  1:2900157/100‑1 (MQ=255)

                                          aaaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAgg                                                         <  2:2160356/100‑1 (MQ=255)

                                           aaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAggg                                                        <  2:3012855/100‑1 (MQ=255)

                                            cgatctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGc                                                       <  2:2549227/96‑1 (MQ=255)

                                            aCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGc                                                       <  1:220029/100‑1 (MQ=255)

                                              tgctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGa                                                     >  2:3039284/1‑100 (MQ=255)

                                               gctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGAt                                                    >  1:3215225/1‑100 (MQ=255)

                                                ctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATc                                                   >  1:2669352/1‑100 (MQ=255)

                                                 tgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCa                                                  <  2:1808623/100‑1 (MQ=255)

                                                  gTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAc                                                 >  2:3044548/1‑100 (MQ=255)

                                                  gTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAc                                                 >  1:2261807/1‑100 (MQ=255)

                                                    tCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAcc                                                <  2:1140349/99‑1 (MQ=255)

                                                    tCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                               >  1:1116471/1‑100 (MQ=255)

                                                    tCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                               >  2:424665/1‑100 (MQ=255)

                                                      cATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcg                                             >  2:3161481/1‑100 (MQ=255)

                                                        tCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcgtc                                           <  2:211664/100‑1 (MQ=255)

                                                          tGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCgg                                         <  1:3783849/100‑1 (MQ=255)

                                                           gcgcACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGc                                        <  2:1640052/100‑1 (MQ=255)

                                                            cgcACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCt                                       <  2:3487397/100‑1 (MQ=255)

                                                            cgcACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCt                                       >  2:1380032/1‑100 (MQ=255)

                                                             gcACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTg                                      >  1:3523324/1‑100 (MQ=255)

                                                             gcACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTg                                      >  2:77973/1‑100 (MQ=255)

                                                             gcACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTg                                      >  2:2715813/1‑100 (MQ=255)

                                                                   cTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGCTCCCATACAGACGATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGAt                                <  2:1145073/100‑1 (MQ=255)

                                                                    tGATTGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCAGACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATc                               <  1:207181/100‑1 (MQ=255)

                                                                     gATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAGGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCa                              <  1:1142363/100‑1 (MQ=255)

                                                                     gATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCa                              <  2:2462485/100‑1 (MQ=255)

                                                                       tggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGTGATCACCTCGTCGGCTGCGCGATCAac                            <  1:3026237/100‑1 (MQ=255)

                                                                       tggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAac                            <  2:2943253/100‑1 (MQ=255)

                                                                         gtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAacac                          >  2:522017/1‑100 (MQ=255)

                                                                          tggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACCATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACt                         >  2:2815065/1‑100 (MQ=255)

                                                                            gAAAGGGGGCCATAGATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACttt                       >  1:505739/1‑100 (MQ=255)

                                                                            gAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACttt                       >  1:1421503/1‑100 (MQ=255)

                                                                               aGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTg                    <  2:378279/100‑1 (MQ=255)

                                                                               aGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTg                    >  1:1338978/1‑100 (MQ=255)

                                                                                gggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTgg                   <  2:1862529/100‑1 (MQ=255)

                                                                                gggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTgg                   >  1:3286950/1‑100 (MQ=255)

                                                                                 ggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                  >  1:1536138/1‑100 (MQ=255)

                                                                                   ctttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGa                >  2:869261/9‑100 (MQ=255)

                                                                                   ctttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGa                >  2:180583/9‑100 (MQ=255)

                                                                                    tttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGaa               >  2:3449985/8‑100 (MQ=255)

                                                                                     ttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAAt              >  2:3304938/7‑100 (MQ=255)

                                                                                     ttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAAt              >  1:2776000/7‑100 (MQ=255)

                                                                                     ttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAAt              >  2:1037161/7‑100 (MQ=255)

                                                                                       aacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATcc            >  1:3756350/5‑100 (MQ=255)

                                                                                          cATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAga         >  1:3737123/1‑100 (MQ=255)

                                                                                            tGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc       <  2:310735/100‑1 (MQ=255)

                                                                                            tGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc       <  1:1053683/100‑1 (MQ=255)

                                                                                             gCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcg      >  1:2951125/1‑100 (MQ=255)

                                                                                             gCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcg      >  2:1600620/1‑100 (MQ=255)

                                                                                             gCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcg      >  1:933695/1‑100 (MQ=255)

                                                                                               cGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGTGATGACCCCGTCGGCTGCGAGATCAACACTTTCTGGTGAATCCAgatcgat    <  1:3376608/100‑1 (MQ=255)

                                                                                                gAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcgatc   <  2:2349173/100‑1 (MQ=255)

                                                                                                gAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcgatc   <  2:1343634/100‑1 (MQ=255)

CAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCA  >  NC_000913/518763‑518959

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