BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	362,884	A→G	71.8%	L99S (TTA→TCA)	lacY ←	lactose permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	362,884	0	A	G	71.8%	232.0 / 93.8	146	L99S (TTA→TCA)	lacY	lactose permease
Reads supporting (aligned to +/- strand): ref base A (21/20); new base G (59/45); total (81/65)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

CTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAGAA  >  NC_000913/362786‑362982

cTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTaa                                                                                                   >  1:591784/1‑100 (MQ=255)

 tttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACAACCAACAATCGATCCTACTAAAATGTTGTATTGTGAc                                                                                                  <  1:1803613/100‑1 (MQ=255)

  ttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACa                                                                                                 >  1:443986/1‑100 (MQ=255)

  ttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACa                                                                                                 >  1:2270630/1‑100 (MQ=255)

   tctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAg                                                                                                <  1:963497/100‑1 (MQ=255)

    ctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGt                                                                                               >  2:1632957/1‑100 (MQ=255)

     tcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTg                                                                                              >  2:2455034/1‑100 (MQ=255)

      cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTgg                                                                                             >  1:1112532/1‑100 (MQ=255)

      cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTgg                                                                                             >  2:2812026/1‑100 (MQ=255)

      cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTgg                                                                                             <  1:2954757/100‑1 (MQ=255)

      cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTgg                                                                                             >  1:1834805/1‑100 (MQ=255)

      cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTgg                                                                                             >  2:1145798/1‑100 (MQ=255)

        aTAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGcc                                                                                           <  2:952655/100‑1 (MQ=255)

        aTAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGcc                                                                                           >  1:117041/1‑100 (MQ=255)

           aaTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCg                                                                                         >  2:602321/1‑99 (MQ=255)

           aaTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGa                                                                                        >  2:1443232/1‑100 (MQ=255)

             tGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAg                                                                                      <  1:1373639/100‑1 (MQ=255)

             tGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAg                                                                                      <  2:88514/100‑1 (MQ=255)

              gCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGa                                                                                     >  2:3137409/1‑100 (MQ=255)

                ctctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATa                                                                                   >  2:314900/1‑100 (MQ=255)

                 tctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGATAGTGGCCCGTAGATac                                                                                  <  1:1487378/100‑2 (MQ=255)

                  ctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaa                                                                                 <  2:33092/100‑1 (MQ=255)

                  ctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaa                                                                                 <  2:34494/100‑1 (MQ=255)

                   tACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaa                                                                                >  2:1243155/1‑100 (MQ=255)

                    aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaaa                                                                               <  1:2115201/100‑1 (MQ=255)

                    aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaaa                                                                               <  2:664474/100‑1 (MQ=255)

                     ctgctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAAt                                                                              >  2:784120/1‑100 (MQ=255)

                        ctgGCGCACCGTTGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaaa                                                                           <  1:653846/100‑1 (MQ=255)

                         tgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAg                                                                          <  1:1243151/100‑1 (MQ=255)

                         tgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAg                                                                          <  2:1775865/100‑1 (MQ=255)

                         tgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAg                                                                          <  1:2508499/100‑1 (MQ=255)

                          ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGa                                                                         >  2:1496198/1‑100 (MQ=255)

                             gcACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGGTCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACg                                                                      >  1:1877103/1‑100 (MQ=255)

                             gcACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACg                                                                      <  2:1719219/100‑1 (MQ=255)

                              cACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGaa                                                                        >  1:2209119/1‑97 (MQ=255)

                                 cGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGgcgc                                                                  <  1:1632950/100‑1 (MQ=255)

                                 cGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGgcgc                                                                  <  1:1705555/100‑1 (MQ=255)

                                    cGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaaa                                                               >  1:301606/1‑100 (MQ=255)

                                     gTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAAc                                                              >  2:792930/1‑100 (MQ=255)

                                      ttAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACa                                                             <  2:2047821/100‑1 (MQ=255)

                                      ttAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACa                                                             <  2:2270646/100‑1 (MQ=255)

                                      ttAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTCACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACa                                                             >  1:2604074/1‑100 (MQ=255)

                                      ttAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACa                                                             >  2:1384343/1‑100 (MQ=255)

                                       tAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACAt                                                            >  2:2091530/1‑100 (MQ=255)

                                        aaaaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATc                                                           >  2:1455479/1‑100 (MQ=255)

                                          aaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAc                                                         >  2:877592/1‑100 (MQ=255)

                                          aaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAc                                                         <  2:1851525/100‑1 (MQ=255)

                                           aCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGTCGCAAACATCACt                                                        <  2:371374/100‑1 (MQ=255)

                                            cAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTa                                                       <  1:1611524/100‑1 (MQ=255)

                                             aaaaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTaa                                                      >  2:493593/1‑100 (MQ=255)

                                               aaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACa                                                    >  1:423098/1‑100 (MQ=255)

                                                aGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACAt                                                   >  1:1203415/1‑100 (MQ=255)

                                                 gCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATg                                                  >  1:2767181/1‑100 (MQ=255)

                                                  ccTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGc                                                 <  1:2575926/100‑1 (MQ=255)

                                                  ccTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGc                                                 >  1:1894941/1‑100 (MQ=255)

                                                   cTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGcc                                                <  1:1551912/100‑1 (MQ=255)

                                                     aGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                              <  2:697710/100‑1 (MQ=255)

                                                      gATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGt                                             <  2:304692/100‑1 (MQ=255)

                                                      gATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGt                                             >  2:2575782/1‑100 (MQ=255)

                                                      gATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGt                                             >  1:2851165/1‑100 (MQ=255)

                                                      gATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGt                                             >  2:484445/1‑100 (MQ=255)

                                                      gATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGt                                             >  1:1060544/1‑100 (MQ=255)

                                                        tAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaa                                           >  1:836424/1‑100 (MQ=255)

                                                         aaaTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaat                                          <  2:2851188/100‑1 (MQ=255)

                                                            tACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTaataat                                       <  1:696720/100‑1 (MQ=255)

                                                             accaccAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTaataat                                       <  1:555985/99‑1 (MQ=255)

                                                              ccaccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                     >  1:897485/1‑100 (MQ=255)

                                                               caccaACAATCGATCCTACTAACATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGATCGGCGCAAACATCACTAACATGCCGGTAATAATcc                                     >  1:1967345/1‑99 (MQ=255)

                                                                 ccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcaca                                  >  2:1915400/1‑100 (MQ=255)

                                                                 ccaACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcaca                                  <  2:754277/100‑1 (MQ=255)

                                                                 ccaACAATCGATCCTACTAAAAAGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGATAATAATCaaca                                  >  1:1337342/1‑100 (MQ=255)

                                                                  caacaaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCca                                    >  1:273195/1‑97 (MQ=255)

                                                                   aacaaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                                <  1:133824/100‑1 (MQ=255)

                                                                     caaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACGTCACTAACATGCCGGTATTAATCCAcagctg                              <  2:2280634/100‑3 (MQ=255)

                                                                      aaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg                             <  1:996717/100‑1 (MQ=255)

                                                                       aTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGt                            >  1:3133146/1‑100 (MQ=255)

                                                                       aTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGt                            >  1:460589/1‑100 (MQ=255)

                                                                          gATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAtt                         >  2:765723/1‑100 (MQ=255)

                                                                          gATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAtt                         >  2:2797213/1‑100 (MQ=255)

                                                                               tactaAAATGTTGTATTGTGTCAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc                    >  1:72435/1‑100 (MQ=255)

                                                                                actaAAATGTTGTATTGTAACAGTCGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc                    <  1:467445/99‑1 (MQ=255)

                                                                                 ctaAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAg                  >  2:2520832/1‑100 (MQ=255)

                                                                                  taAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGc                 >  1:3087730/1‑100 (MQ=255)

                                                                                  naaaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGc                 >  1:1783102/2‑100 (MQ=255)

                                                                                   naaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGAAGc                 >  1:2261862/2‑99 (MQ=255)

                                                                                   aaaaTGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc                >  2:3979/1‑100 (MQ=255)

                                                                                    aaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGccc               >  2:3180290/1‑100 (MQ=255)

                                                                                      aTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCAGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGa             >  1:2911422/1‑100 (MQ=255)

                                                                                        gttgtATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGt           >  1:59593/1‑100 (MQ=255)

                                                                                         ttgtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGtt          <  2:2748980/100‑1 (MQ=255)

                                                                                           gtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACGAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTg        <  1:430498/100‑1 (MQ=255)

                                                                                           gtattgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTg        >  2:1315222/1‑100 (MQ=255)

                                                                                             attgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTc      <  1:858893/100‑1 (MQ=255)

                                                                                              ttgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa     <  1:737178/100‑1 (MQ=255)

                                                                                              ttgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa     >  2:2249404/1‑100 (MQ=255)

                                                                                              ttgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa     >  2:3257419/1‑100 (MQ=255)

                                                                                              ttgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa     >  2:1010487/1‑100 (MQ=255)

                                                                                               tgtgACAGTGGACCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCCGGTATTTGCGCAGCCCGAGTTTGTCa     <  1:939870/99‑1 (MQ=255)

                                                                                               tgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa     >  1:2137088/1‑99 (MQ=255)

                                                                                               tgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAg    <  2:295605/100‑1 (MQ=255)

CTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAGAA  >  NC_000913/362786‑362982

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: `‑`