BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	518,863	C→T	17.1%	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	17.1%	319.3 / 47.3	139	G96E (GGA→GAA)	ybbO	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (55/59); new base T (13/11); total (69/70)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 140 total aligned reads displayed.

AGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATC  >  NC_000913/518764‑518958

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

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

   cgcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCgg                                                                                              >  2:1687441/1‑100 (MQ=255)

   cgcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCgg                                                                                              >  1:1938515/1‑100 (MQ=255)

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

       tGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCAt                                                                                           >  2:1781774/1‑99 (MQ=255)

        ggTAAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtg                                                                                          <  1:2661881/99‑2 (MQ=255)

         gTGAGCTGGTGTGCGCCGAAAAAGTTGGGGGAAAACTGCTGTTCCATCTGCGCACGGCTGATTGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                        <  1:1549217/100‑1 (MQ=255)

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

          tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtat                                                                                       <  1:413319/100‑5 (MQ=255)

          tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt                                                                                       >  1:2218693/1‑100 (MQ=255)

          tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt                                                                                       >  2:2450713/1‑100 (MQ=255)

          tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                        >  1:2218705/1‑99 (MQ=255)

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

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

           gAGCTGGTGTGCCCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatg                                                                                      <  1:544644/100‑6 (MQ=255)

            tgCTGGTGTGCGCCGAATAAGTTGTCGGAAGACGGCGGTTCCGTCTGCGCACGGCTGAGGGTGGAAGGGGGGCCATACATGCCGGATTCGGTAtgtatgg                                                                                     <  1:1435035/99‑7 (MQ=255)

            aGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaa                                                                                     >  1:367863/1‑100 (MQ=255)

             gCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaa                                                                                     >  1:1291611/1‑99 (MQ=255)

              cTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAg                                                                                   <  2:215791/100‑1 (MQ=255)

               tGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGa                                                                                  >  2:1626861/1‑100 (MQ=255)

                ggTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGAt                                                                                 >  2:2891432/1‑100 (MQ=255)

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

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

                  tgtgCGCCGAAAAAGTTGGCAGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATcc                                                                               <  1:898227/100‑1 (MQ=255)

                    tgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCa                                                                             >  2:1379105/1‑100 (MQ=255)

                    tgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCa                                                                             >  2:2913797/1‑100 (MQ=255)

                      cgcCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGTCGAATCCGGCATTGTTAAAGATCCCATa                                                                           >  1:841606/1‑100 (MQ=255)

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

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

                       gcCGAAAAAGTTGTCGGAAAACTGCTGTTCCATCTGCCCACGGCTGATTGTGGAAAGCGCGCCATACATGCCGAATCCGGCATTGTTaacggtca                                                                               <  1:2505869/95‑7 (MQ=255)

                       gacGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATa                                                                           <  1:1358358/97‑1 (MQ=255)

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

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

                         aaaaaaaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCAt                                                                            <  1:2792498/94‑1 (MQ=255)

                          gAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGa                                                                       >  2:512046/1‑100 (MQ=255)

                          gAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCAGAATCCGACATTGTGAAAGATCCCATACAGa                                                                       >  1:1701765/1‑100 (MQ=255)

                               aGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAAtt                                                                  >  1:2366715/1‑100 (MQ=255)

                                 ttGGCGGAAAACTGCTGTTCCATCTGCGCACGGCGGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTAt                                                                <  2:1290816/100‑1 (MQ=255)

                                  tGGCGGAACACTGCTGTTCCATCTGCGCACGGCTGAGGGTGGAAAGAGGGCCATACATGCCGAATCCGGCATTGTTAAAGCTCCCATACAGACAATTATc                                                               <  1:2058280/100‑1 (MQ=255)

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

                                    gCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCgg                                                             >  2:1887923/1‑100 (MQ=255)

                                     cGGAAAGCTCCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCg                                                              <  1:2564341/98‑1 (MQ=255)

                                     cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCgg                                                             >  1:107341/1‑99 (MQ=255)

                                      ggAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTc                                                           <  1:1278655/100‑1 (MQ=255)

                                      ggAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTc                                                           <  2:1031180/100‑1 (MQ=255)

                                      ggAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTc                                                           >  2:1084742/1‑100 (MQ=255)

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

                                           aCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTCTCGGTCAGGGc                                                      >  1:2745335/1‑100 (MQ=255)

                                               ctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATc                                                  <  1:566118/100‑1 (MQ=255)

                                               ctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATc                                                  <  2:1521880/100‑1 (MQ=255)

                                               ctgTTCCAGCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCACACAGACAATTATCGGTCAGGGCGAtn                                                  <  1:1072818/100‑2 (MQ=255)

                                                tgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATc                                                  <  2:1193807/99‑1 (MQ=255)

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

                                                  ttCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAcc                                               <  2:597595/100‑1 (MQ=255)

                                                   tCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATAGATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                              >  1:1800009/1‑100 (MQ=255)

                                                     cATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcg                                            <  2:2218718/100‑1 (MQ=255)

                                                     cATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcg                                            <  2:2218706/100‑1 (MQ=255)

                                                        cTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcgtcg                                         <  1:1645089/100‑1 (MQ=255)

                                                        cTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcgtcg                                         >  1:1279125/1‑100 (MQ=255)

                                                        cTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcgtcg                                         <  1:673185/100‑1 (MQ=255)

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

                                                         tGCGCACGGCTGATGGTGGAAAGGGGGCCATAAATGCCGAATCCGGAATTGTTAAAGATCACATAAAGACAATTATCGGTAAGGGCGATAACCACGTCgg                                        >  1:116582/1‑100 (MQ=255)

                                                          gcgcACGGCTGATGGGGGAAAGGGGGCCATACATGTCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGc                                       <  2:841609/100‑1 (MQ=255)

                                                             cACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgc                                    >  2:2270131/1‑100 (MQ=255)

                                                             cACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgc                                    >  2:642652/1‑100 (MQ=255)

                                                              aCGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgcg                                   >  1:356582/1‑100 (MQ=255)

                                                                  cTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgcgc                                  <  1:1901154/97‑1 (MQ=255)

                                                                   tGATGGTGGAAAGGGAGGCATACATGCAGAATCCGGCATTGTTAAAGATCCCATACCGACAATTATCGGTCAGGGTGATCACCTCGTCGTCTGCCCGATc                              <  1:2430737/100‑1 (MQ=255)

                                                                       ggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAaca                          <  2:2632497/100‑1 (MQ=255)

                                                                       ggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAaca                          <  2:1567450/100‑1 (MQ=255)

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

                                                                         tggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACt                        >  2:2332338/1‑100 (MQ=255)

                                                                             aaGGGGGCCATACATGCCGATTCCGGCATTGTTAAAGATCCCATACAGTCAATTATCGGTCATTGCGATGATACCGCCGGCTCCGCGATCAACACTTTCt                    <  1:755416/100‑1 (MQ=255)

                                                                              aGGGGGCCATCCATGCCGAATCCGGCATTGTTAAAGAACCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTg                   <  2:2734019/100‑1 (MQ=255)

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

                                                                               gggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCAATCAACACTTTCTgg                  <  1:1379100/100‑1 (MQ=255)

                                                                                ggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                 <  1:1893033/100‑1 (MQ=255)

                                                                                ggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                 <  1:2915491/100‑1 (MQ=255)

                                                                                atctttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                 >  2:2441323/11‑100 (MQ=38)

                                                                                  ggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGa               >  2:1056901/1‑100 (MQ=255)

                                                                                   tttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGaa              >  1:2872824/8‑100 (MQ=255)

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

                                                                                     taacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAAt             >  1:2205026/6‑99 (MQ=255)

                                                                                      aacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                 >  1:192586/5‑94 (MQ=255)

                                                                                      aacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATcc           >  2:85204/5‑100 (MQ=255)

                                                                                       tACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCa          >  2:521522/1‑100 (MQ=255)

                                                                                         cATGCCGAATTCTGCATTGTTAAAGATCCCATCCAGACAATGATCTTTCAGGTTGATGATCCGGCGGGCTCCGCGATCTACACTTTCTGGTGAATcc           <  1:1829962/97‑1 (MQ=255)

                                                                                         cATGCCGAATGCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTAAGCTGAGCGACCAACACTTTCTGGTGAATCCaaa        >  1:2095152/1‑98 (MQ=255)

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

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

                                                                                          aTGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgat       >  1:2116175/1‑100 (MQ=255)

                                                                                           tGCCGAATCCGGTATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc      <  1:970494/100‑1 (MQ=255)

                                                                                           tGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATGACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc      <  1:110940/100‑1 (MQ=255)

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

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

                                                                                            gCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcg     <  2:320506/100‑1 (MQ=255)

                                                                                             ccGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcga    >  2:1991658/1‑100 (MQ=255)

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

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

AGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATC  >  NC_000913/518764‑518958

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