BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	518,863	C→T	15.5%	G96E (GGA→GAA)	ybbO ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	518,863	0	C	T	15.5%	430.7 / 56.3	172	G96E (GGA→GAA)	ybbO	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (63/82); new base T (9/18); total (72/100)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 174 total aligned reads displayed.

GGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAAC  >  NC_000913/518765‑518961

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

 gcgcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCg                                                                                                  <  2:2750242/100‑1 (MQ=255)

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

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

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

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

    cATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCGTACATGCCGGATTCGGCa                                                                                               <  2:1365427/100‑1 (MQ=255)

    cATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCa                                                                                               >  2:499498/1‑100 (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:2142536/100‑4 (MQ=255)

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

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

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

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

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

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

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

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

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

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

            gCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatggc                                                                                       <  1:2779341/100‑8 (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)

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

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

                     cgcCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATa                                                                              >  1:3221909/1‑100 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                   gCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCgg                                                                <  1:2580227/100‑1 (MQ=255)

                                    cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGt                                                               >  1:3754796/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)

                                        aaaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGGAAGGGGGCCATACATGCCGAATCCGTCATTTTTAAAGAGCCCATACAGGCAATGATCGTTCAg                                                            >  2:1125916/1‑99 (MQ=255)

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

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

                                          cgatctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGc                                                         <  2:2549227/96‑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:2487760/1‑100 (MQ=255)

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

                                                gTTCCATCTGCGCACGGTTGATGGTGGAAAGGGGGCGAGACATGCAGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAc                                                   >  1:1122864/1‑100 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                          cgcACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCt                                         <  2:3487397/100‑1 (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)

                                                                  tGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATc                                 <  1:1903219/100‑1 (MQ=255)

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

                                                                   gATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCa                                <  1:2156027/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)

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

                                                                             aGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTg                      <  2:378279/100‑1 (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)

                                                                                    cATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACCCTTTCTGGTGAATc               >  1:3474286/1‑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:2153939/1‑100 (MQ=255)

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

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

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

                                                                                               aaTCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCa    >  1:2997822/1‑100 (MQ=255)

                                                                                               aaTCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCa    <  2:1052010/100‑1 (MQ=255)

                                                                                                aTCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCaa   >  1:427934/1‑100 (MQ=255)

                                                                                                 tCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAac  <  2:2030562/100‑1 (MQ=255)

GGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAAC  >  NC_000913/518765‑518961

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