BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	362,884	A→G	16.2%	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	16.2%	233.9 / 39.1	107	L99S (TTA→TCA)	lacY	lactose permease
Reads supporting (aligned to +/- strand): ref base A (51/38); new base G (9/8); total (60/47)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-01

Only 100 of 121 total aligned reads displayed.

CTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAG  >  NC_000913/362786‑362980

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

 tttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAAc                                                                                                <  2:2229048/100‑1 (MQ=255)

  ttCTTAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGttgtattgta                                                                                                  <  2:1004195/97‑1 (MQ=255)

  ttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACa                                                                                               >  1:991233/1‑100 (MQ=255)

   tctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAg                                                                                              <  1:601265/100‑1 (MQ=255)

   tctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAg                                                                                              >  1:1481481/1‑100 (MQ=255)

     tcAATAAATGCCTCTACTGCTGGAGCACCGGAGTTAAAACAAAATCCTAGATAAATACCACCAACAAACAATCCTACTAAAATGTAGTATTGTAACa                                                                                               >  2:2279503/1‑97 (MQ=255)

      cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGttg                                                                                           >  1:1675549/1‑98 (MQ=255)

        aTAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGcc                                                                                         <  1:1106150/100‑1 (MQ=255)

        aTAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGAATTGTAACAGTGGc                                                                                          >  1:2332235/1‑99 (MQ=255)

         tAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGccc                                                                                        >  1:2002436/1‑100 (MQ=255)

          aaaTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCg                                                                                       >  1:1789362/1‑100 (MQ=255)

           aaTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGa                                                                                      >  1:249910/1‑100 (MQ=255)

            aTGCCTCTACTGCTGGCACACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGaa                                                                                     <  2:1147292/100‑1 (MQ=255)

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

              gCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAg                                                                                    <  2:1882522/99‑1 (MQ=255)

                ctctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATa                                                                                 <  2:1668016/100‑1 (MQ=255)

                 tctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaa                                                                                >  1:2027405/1‑100 (MQ=255)

                 tcgacggttgGCGCACGGGCGTTAAAACAAAAGCCTATATAAATAGCGCCAACAAGTGATCCGACTAAAATGTTGTAGTGTAACGGTGGCCCGAAGATaa                                                                                <  2:2175959/92‑1 (MQ=255)

                  ctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaa                                                                               >  1:83942/1‑100 (MQ=255)

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

                    aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaaaa                                                                             >  1:1597307/1‑100 (MQ=255)

                      ttctggattaCCGGCGTTAAAACAAAAGCCGATATATGTTCCTCCAACAGTCGATCCTTCTAAATTGTTGTATTGTAATAGTGGTCCAAAGATAAAAAt                                                                            <  2:950886/90‑1 (MQ=255)

                          ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGCTCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGa                                                                       >  1:214182/1‑100 (MQ=255)

                          ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTAGCCCGAAGAAAAAAATAAAGa                                                                       >  1:679780/1‑100 (MQ=255)

                          ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGa                                                                       >  2:182304/1‑100 (MQ=255)

                           gcgcACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGaa                                                                      >  1:769347/1‑100 (MQ=255)

                              cACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACgg                                                                   >  2:1062214/1‑100 (MQ=255)

                                ccGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACgg                                                                   >  2:1496246/1‑98 (MQ=255)

                                   gCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaa                                                              >  1:487487/1‑100 (MQ=255)

                                      ttAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACa                                                           <  1:1547899/100‑1 (MQ=255)

                                        aaaaCAAAAGCCTAGATAAATACCACAAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGACCCGAAGATAAAAATAAAGAACGGCACAAACATc                                                         >  2:1046081/1‑100 (MQ=255)

                                         aaaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCa                                                        >  1:2393121/1‑100 (MQ=255)

                                         aaaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGAGGCCCGAAGATAAAAATAAAGAACAGCGCAAACATCa                                                        >  2:955434/1‑100 (MQ=255)

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

                                          aaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATcgc                                                       >  2:1451040/1‑98 (MQ=255)

                                          aaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAc                                                       <  1:1244246/100‑1 (MQ=255)

                                          aaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAc                                                       >  2:606716/1‑100 (MQ=255)

                                             aaaaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTaa                                                    >  1:1195731/1‑100 (MQ=255)

                                              aaaGCCTAGATAAATACCACCAACCGTCGATCCTACTAAAATGTTGTATTGTGACTGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAAc                                                   <  1:536818/100‑1 (MQ=255)

                                               aaGCCTAGATAAATACCCCCCACAATCGACCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATCAAGAACGGCGCAAACATCACTAACa                                                  <  2:1191369/100‑1 (MQ=255)

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

                                                aGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACAt                                                 <  2:1597287/100‑1 (MQ=255)

                                                aGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACAt                                                 >  1:953770/1‑100 (MQ=255)

                                                    tAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCg                                             >  2:586912/1‑100 (MQ=255)

                                                    tAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCg                                             >  1:1599666/1‑100 (MQ=255)

                                                    tAGATAAATACCACCAACAATCGAGCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCg                                             <  1:1213021/100‑1 (MQ=255)

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

                                                     aGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                            >  2:311882/1‑100 (MQ=255)

                                                     aGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                            <  2:249907/100‑1 (MQ=255)

                                                     aGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                            <  1:1451056/100‑1 (MQ=255)

                                                      gATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGt                                           <  2:1223894/100‑1 (MQ=255)

                                                        tAAATACCACCAACAATCGATCCTACTAAAATGATGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAGCGGCGCAAACATCACTAACATGCCGGtaa                                         >  2:1440072/1‑100 (MQ=255)

                                                         aaaTACCACCAACAATCGATCCTACTAAATTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaat                                        <  2:985819/100‑1 (MQ=255)

                                                            tACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTaataat                                     <  1:412343/100‑1 (MQ=255)

                                                            tACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAAAACGGCGCAAAAATCACTAACATGCCGGTaataat                                     >  2:1098175/1‑100 (MQ=255)

                                                             accaccAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATc                                    <  1:998682/100‑1 (MQ=255)

                                                              ccaccaACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                   <  2:796888/100‑1 (MQ=255)

                                                              ccaccaACAATATAACCTACTAAAATATTGTATTCTAACAGTGGCCCGACGATAACAATAAAGAACGGCGCAAACATCACTAACATGCCGGTaataataa                                   <  2:1304279/100‑3 (MQ=255)

                                                                accaACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac                                 >  1:373882/1‑100 (MQ=255)

                                                                accaACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac                                 >  1:380297/1‑100 (MQ=255)

                                                                accaACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCca                                  >  1:764126/1‑99 (MQ=255)

                                                                 ccgacaaTCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAGAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                   <  2:1216076/94‑1 (MQ=255)

                                                                  caacaaTCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcaca                                >  2:2092204/1‑99 (MQ=255)

                                                                   aacaaTCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                              <  1:594897/100‑1 (MQ=255)

                                                                    acaaTCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagca                             >  2:2269789/1‑100 (MQ=255)

                                                                      aaTCGATCCTACTAAAATGTTTTATTGGAACAGTGTCCCTAAGATAAAAATAAAGAACGTCGCAAACATCACTAACATGCAGGTAATAATCCACAGCAgg                           <  2:2007499/100‑1 (MQ=255)

                                                                      aaTCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg                           <  2:1761975/100‑1 (MQ=255)

                                                                        tCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTa                         <  1:532146/100‑1 (MQ=255)

                                                                         tgATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAt                        >  1:235682/2‑100 (MQ=255)

                                                                         cGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAt                        >  1:151120/1‑100 (MQ=255)

                                                                         cGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAt                        >  2:1155378/1‑100 (MQ=255)

                                                                         cGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAAAAATCCACAGCAGGTAt                        >  2:2362879/1‑100 (MQ=255)

                                                                          gATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTa                         >  1:1094108/1‑98 (MQ=255)

                                                                          gATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAtt                       >  2:1933366/1‑100 (MQ=255)

                                                                          gATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATACAcagcatcatgc                       >  2:406561/1‑94 (MQ=255)

                                                                            tCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTg                     <  1:411434/100‑1 (MQ=255)

                                                                            tCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTg                     >  1:1689871/1‑100 (MQ=255)

                                                                             ccTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAttt                      >  1:49012/1‑98 (MQ=255)

                                                                             ccTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgc                    >  2:1373250/1‑100 (MQ=255)

                                                                              ctactaAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcg                   <  1:757698/100‑1 (MQ=255)

                                                                                  tcaaaaGTTGTCTTGTAACAATGGCCACAAAAAAAAAATAAAAAACTCCTTAAACATCACTAACATGCCGGTAAACATCCACAGCACGCAtt                       >  2:465927/3‑92 (MQ=255)

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

                                                                                  taAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGc               <  1:885684/100‑1 (MQ=255)

                                                                                   aaaaTGTTGTATTGTTACAGTGGCCAAAAGATAAAAATAAAGAACGGCGCAAACAGCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGc               <  2:1966294/99‑1 (MQ=255)

                                                                                   aaaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc              >  2:566388/1‑100 (MQ=255)

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

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

                                                                                   aaaaTGTTGTATTGTAACAGTGGCAAGAAGATAAAAATAAAGAAAGGCGCAAACATCACTAACATGCCGGTAATAATACACAGCAAGTATTTGCGCAGc               >  2:2141283/1‑99 (MQ=255)

                                                                                   aaaaTGTTGTATTGTAACAGTGACCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc              >  2:2063891/1‑100 (MQ=255)

                                                                                     aaTGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc              <  2:2027383/98‑1 (MQ=255)

                                                                                        gttgtATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGt         <  2:1599646/100‑1 (MQ=255)

                                                                                         ttgtattgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGtt        <  1:613387/100‑1 (MQ=255)

                                                                                           gtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTg      <  2:429967/100‑1 (MQ=255)

                                                                                           gtattgtaACAGTGACCCTAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCATGTATTTGCGCAGCCCGAGttgg      <  2:621180/100‑3 (MQ=255)

                                                                                             attgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTc    >  2:1048668/1‑100 (MQ=255)

                                                                                              ttgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGTCCGAGTTTGTCa   <  2:2332207/100‑1 (MQ=255)

                                                                                              ttgtaACAATGGCCCGAATATAAAAATAAAGAACGGCTCAAACATCACTAACATGCCGTTAATAATCCACATCATTTATTTGCGCAGCCCGAGTTTg      <  2:1820784/97‑1 (MQ=255)

                                                                                              tgttaACAAAGGACCCAAAATAAAAATAAAGAACGGCGCCATCATCACACACATCACGGTAATAATCGACAGCAAGTATTTGGGCAGaccaa           >  2:1350913/4‑90 (MQ=255)

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

CTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAG  >  NC_000913/362786‑362980

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