BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	362,884	A→G	100%	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	96.5%	351.9 / 1.1	115	L99S (TTA→TCA)	lacY	lactose permease
Reads supporting (aligned to +/- strand): ref base A (1/3); new base G (50/58); total (52/63)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 127 total aligned reads displayed.

TTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAG  >  NC_000913/362788‑362980

ttCTCGATAAATGTCTCTACTGCTGTGGCACAGGCGTTAAAACAAAAGCCTAGATAAATTCTTCCTACAATCGATTCTACTATAATGTTGTATTGTGACa                                                                                               <  2:2072338/100‑1 (MQ=255)

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

 tctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAg                                                                                              >  1:842655/1‑100 (MQ=255)

 tctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACAACCAACAATCGATCCTACTAAAATGTTGTATTGTGACa                                                                                               >  1:431264/1‑99 (MQ=255)

   tcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACTAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTg                                                                                            >  2:1497371/1‑100 (MQ=255)

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

    cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTg                                                                                            >  2:838365/1‑99 (MQ=255)

    cAATAAATGCCTCTACTGCTGGAGCCACGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTCCTAAAATGTTGTATTGTGACAGTgg                                                                                           >  2:152289/1‑100 (MQ=255)

     aaTAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGc                                                                                          <  1:2763043/100‑1 (MQ=255)

      aTAAATGCCTCTACTGCTTGCGCACCGGCGTTAAAATAATTTTTTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGcc                                                                                         <  1:849535/100‑1 (MQ=255)

         aaTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCAATACTACTAAAATGTTAAATTGTGACAGTGGCCCGa                                                                                      >  1:129789/1‑100 (MQ=255)

          aTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGaa                                                                                     <  1:2416964/100‑1 (MQ=255)

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

            gCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTATATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAg                                                                                    <  1:93602/99‑1 (MQ=255)

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

             ccTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGAt                                                                                  >  1:2063857/1‑100 (MQ=255)

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

               tctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaa                                                                                <  1:1119601/100‑1 (MQ=255)

               tctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaa                                                                                <  2:1901129/100‑1 (MQ=255)

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

                 tACTGCTAGCGCACCCCCGTTAAAACAAAAGCATAAATAAATACCACCAAAAATTGATACAACTGAAATGTTGTATTGTGACAGTGGCACGagtacaaa                                                                               >  2:1451085/1‑92 (MQ=255)

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

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

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

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

                   ctgctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAAt                                                                            <  2:2777313/100‑1 (MQ=255)

                     gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaa                                                                          >  1:2095924/1‑100 (MQ=255)

                     gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaa                                                                          >  2:2867034/1‑100 (MQ=255)

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

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

                       tgGCGCACCGGCGTTAAAACAAAAGCCGAGATAAATACCACCAAAAATCGAACCTAAAAAAATGATGTATTGTGAAAATGTCCCGAAGAGAAAAATAAAg                                                                        >  1:892583/1‑100 (MQ=255)

                       tgGCGCACCGGCGTTAAAACAAAAGACTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAg                                                                        <  2:794944/100‑1 (MQ=255)

                         gcgcACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGaa                                                                      <  2:3038397/100‑1 (MQ=255)

                            cACCGGCGTTAAAACAAAAGTCAAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACgg                                                                   <  1:1689329/100‑1 (MQ=255)

                            cACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACgg                                                                   <  2:613182/100‑1 (MQ=255)

                                ggCGTTAAAACAAAAGCCTAGATAAATACCTCCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCa                                                               <  1:112907/100‑1 (MQ=255)

                                ggCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGgcgc                                                                <  1:1727055/99‑1 (MQ=255)

                                 gCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaa                                                              <  2:3013259/100‑1 (MQ=255)

                                  cGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCaaa                                                             >  2:2214759/1‑100 (MQ=255)

                                   gTTAAATTATAAGCCTAGATAAATACCACCATCAATCTATCCTACTAAAATGTTTTATTGTGACAGTGGCTTTAAGATAATAATAAAGAACGGCGCaaa                                                             <  2:2722323/99‑1 (MQ=255)

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

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

                                      aaaaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACAt                                                          >  2:2169892/1‑99 (MQ=255)

                                      aaaaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTATGACAGTGGCCCAAAGATAAAAATAAAGAACGGCGCAAACAt                                                          >  1:964282/1‑99 (MQ=255)

                                            aaaGCCTAGATAAATACCACCTACATTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAAc                                                   <  1:2592163/100‑1 (MQ=255)

                                            aaaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAAc                                                   <  2:2174357/100‑1 (MQ=255)

                                            aaaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAAc                                                   <  2:1543805/100‑1 (MQ=255)

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

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

                                                 cTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGcc                                              >  2:750228/1‑100 (MQ=255)

                                                 cTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGcc                                              <  2:2857466/100‑1 (MQ=255)

                                                 cTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGcc                                              <  1:2867066/100‑1 (MQ=255)

                                                   aGATAAATACCACCAACAATCGATCCTACTAAATTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                            <  1:987322/100‑1 (MQ=255)

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

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

                                                     aTAAATACCACCAACAATCGATCCTACTAAAATGTTGGATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGTCGCAAATATCACTAACATGCCGGta                                          <  2:916442/100‑1 (MQ=255)

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

                                                        aaTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGgaata                                       >  2:2830641/1‑100 (MQ=255)

                                                         aTACCACCAACAATCGATCCTACTAAAATGTCGTATTCTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaataa                                      <  2:2043250/100‑1 (MQ=255)

                                                          tACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTaataat                                     >  2:930931/1‑100 (MQ=255)

                                                            ccaccaACGATCGATCCTACTAAAATGTTGTATTGTTACAGTGGCCCGAAGATAAAAAGAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                   <  1:1296765/100‑1 (MQ=255)

                                                            ccaccaACAATCTATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATc                                    <  1:312704/99‑1 (MQ=255)

                                                            ccaccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaataa                                      >  2:1506984/1‑97 (MQ=255)

                                                            ccaccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                                   <  1:565277/100‑1 (MQ=255)

                                                            ccaccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACAGCGCAAACATCACTAACATGCCGGTAATAATc                                    >  2:144817/1‑99 (MQ=255)

                                                             caccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCca                                  <  1:2625072/100‑1 (MQ=255)

                                                             caccaACAATCGATCCTACTAAAATATTGTATTGTGACAGTGGCCCGAAGATAAAAAGAAAGAAAGGAGCAAACATCAATAACATGACAGtaatca                                      >  2:2010794/1‑94 (MQ=255)

                                                              accaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac                                 >  1:2251165/1‑100 (MQ=255)

                                                               ccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGTTAAAAATAAAGAACGGCGCAAAAATCAATAAAAAGCCGGTAATAATTcaca                                >  1:484803/1‑100 (MQ=255)

                                                                caacaaTCGATCCTACTAAAATGTTGTATTGTGACAGT‑GCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                              >  2:1240402/1‑100 (MQ=255)

                                                                   caaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagcag                            >  2:427329/1‑100 (MQ=255)

                                                                   caaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagcag                            >  1:2428776/1‑100 (MQ=255)

                                                                    aaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg                           >  1:298424/1‑100 (MQ=255)

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

                                                                    aaTCGATCATACTAAAATGTTGTATTATGACAGTGGCCCGAAGATAAAAATAAAGAACGTCGCAAACATCACTAACATGCCGGTAATAATCCAcagcag                            <  2:946981/99‑1 (MQ=255)

                                                                     aTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCTCGACGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGt                          <  1:136007/100‑1 (MQ=255)

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

                                                                           ccTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgc                    >  1:2249945/1‑100 (MQ=255)

                                                                           ccTACTAAAATGTTGTATTGTAACAGTGGCACGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCTGGTAATAATCCACAGCAGGTATTTgc                    <  1:1842333/100‑1 (MQ=255)

                                                                            ctactaAAATGTTGTATTGTTACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcg                   <  1:2743723/100‑1 (MQ=255)

                                                                            ctactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcg                   <  2:1741725/100‑1 (MQ=255)

                                                                             tactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc                  >  1:1650120/1‑100 (MQ=255)

                                                                             tactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc                  <  1:511257/100‑1 (MQ=255)

                                                                              actaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCa                 <  1:2758348/100‑1 (MQ=255)

                                                                               ctaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCa                 >  2:1638653/1‑99 (MQ=255)

                                                                               ctaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAg                >  1:1061345/1‑100 (MQ=255)

                                                                                taAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAATAACATACCGGTAATAATCCACAGCAGGTATTTGCGCAg                >  1:66532/1‑99 (MQ=255)

                                                                                 aaaaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc              <  1:559065/100‑1 (MQ=255)

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

                                                                                      gttgtATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAAACCACAGCAGGTATTTGGGAAGACCGa           >  2:1495895/1‑98 (MQ=255)

                                                                                       ttgtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGtt        >  1:708016/1‑100 (MQ=255)

                                                                                       ttgtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACAAACATGCCGGTAATAATACACAGCAGGTATTTGCGCAGCCCGAgat        >  1:2571108/1‑98 (MQ=255)

                                                                                        tgtattgttACATTGGCTTGAATTTAATATTTAAGAATGGCTTAAACATCACTAACATGCCGGTAATAAGCCACAGCAGGTATTTGCGCAGCCCGAGttt       <  1:930937/100‑1 (MQ=255)

                                                                                        tgtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGtt        <  2:107116/99‑1 (MQ=255)

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

                                                                                          tattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGt     >  2:768326/1‑100 (MQ=255)

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

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

                                                                                             tgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAg  >  2:899618/1‑100 (MQ=255)

                                                                                             tgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAg  <  1:469835/100‑1 (MQ=255)

TTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAG  >  NC_000913/362788‑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: `‑`