BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,902,494	G→T	17.2%	G200V (GGA→GTA) ‡	ygcE →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,902,494	0	G	T	17.2%	465.4 / 87.5	186	G200V (GGA→GTA)	ygcE	putative kinase
Reads supporting (aligned to +/- strand): ref base G (79/75); new base T (7/25); total (86/100)
Fisher's exact test for biased strand distribution p-value = 3.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-02

Only 100 of 192 total aligned reads displayed.

TTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACA  >  NC_000913/2902395‑2902592

ttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAgg                                                                                                    <  2:1706/100‑1 (MQ=255)

tgaaTGTCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAgg                                                                                                    <  1:3667372/98‑1 (MQ=255)

   atatCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACt                                                                                                 >  2:1130474/1‑100 (MQ=255)

      tCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGcc                                                                                              <  1:2867356/100‑1 (MQ=255)

       cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCg                                                                                             >  2:3139448/1‑100 (MQ=255)

       cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCg                                                                                             <  2:1988458/100‑1 (MQ=255)

       cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCg                                                                                             >  2:1638605/1‑100 (MQ=255)

        ccccGTCATATGCTGTTTGATGTGCAACTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCgg                                                                                            <  2:1880606/100‑1 (MQ=255)

        ccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATGACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCgg                                                                                            >  2:3411974/1‑100 (MQ=255)

         cccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggg                                                                                           <  2:3040567/100‑2 (MQ=255)

         cccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggg                                                                                           <  2:1565026/100‑2 (MQ=255)

         cccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggg                                                                                           <  2:1410322/100‑2 (MQ=255)

         cccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgt                                                                                           <  1:1030027/100‑1 (MQ=255)

         cccGTAATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggg                                                                                           <  2:1200064/100‑2 (MQ=255)

          ccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggga                                                                                          <  1:1701747/100‑3 (MQ=255)

          ccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggga                                                                                          <  2:2846531/100‑3 (MQ=255)

          ccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggga                                                                                          <  2:498444/100‑3 (MQ=255)

          ccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggga                                                                                          <  2:3925929/100‑3 (MQ=255)

          ccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgtt                                                                                          >  1:1849108/1‑100 (MQ=255)

           cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGTCGACACATTTCCCGGCAGGACTGCCg                                                                                             <  1:55611/96‑1 (MQ=255)

           cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaa                                                                                         <  1:1396610/100‑4 (MQ=255)

           cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaa                                                                                         <  2:905474/100‑4 (MQ=255)

           cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaa                                                                                         <  1:3493051/100‑4 (MQ=255)

           cGTCATATGCTGTTTGATGTGCAAATGCCGGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggga                                                                                          <  1:2165015/99‑3 (MQ=255)

            gTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaa                                                                                        <  2:715257/100‑5 (MQ=255)

            gTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgt                                                                                        <  1:4135915/100‑1 (MQ=255)

              cATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTTGCGACACATTTCCCGGCAGGACTGCCGGttgttt                                                                                      <  1:147813/100‑1 (MQ=255)

              cATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatg                                                                                      <  2:1490122/100‑7 (MQ=255)

              cATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGttgttt                                                                                      >  2:2102094/1‑100 (MQ=255)

               atatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGttgtttg                                                                                     <  2:1269478/100‑1 (MQ=255)

                tatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatgtg                                                                                    <  1:1730346/100‑9 (MQ=255)

                  tGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatgtgtc                                                                                  <  1:1267469/100‑11 (MQ=38)

                   gCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACACGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTacc                                                                                 >  2:2472361/1‑100 (MQ=255)

                   gCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTacc                                                                                 >  2:914367/1‑100 (MQ=255)

                    cTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTacca                                                                                >  2:2215150/1‑100 (MQ=255)

                    cTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTacca                                                                                >  1:383745/1‑100 (MQ=255)

                      gTTTGATGTACAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTaccacc                                                                              <  1:1202634/100‑1 (MQ=255)

                        ttGATGTGCAAATGCCTGGCACCGGCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAg                                                                            <  2:3378929/100‑1 (MQ=255)

                          gATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTg                                                                          <  1:2069791/100‑1 (MQ=255)

                           aTGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGa                                                                         <  1:615803/100‑1 (MQ=255)

                            tgtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                        >  2:2997365/1‑100 (MQ=255)

                            tgtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                        >  2:2205825/1‑100 (MQ=255)

                            tgtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                        <  1:2600788/100‑1 (MQ=255)

                            tgtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                        >  2:3467977/1‑100 (MQ=255)

                             gtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATa                                                                       >  2:4314381/1‑100 (MQ=255)

                              tgCATATGCCTGGCACCGTCCTCGTACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGTCAGGACTGCCGGTTGTTTGTACCACCAGTGATa                                                                       <  1:3770147/99‑1 (MQ=255)

                                cAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATaaa                                                                     >  1:2731907/1‑99 (MQ=255)

                                 aaaTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAcc                                                                   <  2:2915827/100‑1 (MQ=255)

                                 aaaTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAcc                                                                   <  2:3846670/100‑1 (MQ=255)

                                  aaggCCTGGCACCGTTTTGGGACATTTCACACCACAAGCCTCACTTGCGACACATTTCTGGGCATGACTGCCGGTTGTTTGTACCACCAGTGATAAAcc                                                                   <  1:2381018/96‑1 (MQ=255)

                                   aTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCgg                                                                 <  2:2084718/100‑1 (MQ=255)

                                   aTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCgg                                                                 <  1:102369/100‑1 (MQ=255)

                                    tGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGt                                                                <  2:957448/100‑1 (MQ=255)

                                    tGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGt                                                                >  1:2226152/1‑100 (MQ=255)

                                     gCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTa                                                               >  2:2411515/1‑100 (MQ=255)

                                     gCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTa                                                               <  1:900618/100‑1 (MQ=255)

                                      ccTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTa                                                               >  1:3017159/1‑99 (MQ=255)

                                         ggCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACGGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAg                                                           >  1:2392481/1‑100 (MQ=255)

                                          gCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGc                                                          <  1:4130947/100‑1 (MQ=255)

                                             ccGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGctct                                                       >  2:1220267/1‑100 (MQ=255)

                                             ccGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGctc                                                        >  1:3417535/1‑99 (MQ=255)

                                               gTCCTCGGACATATCACACCACAAGTCGCTCTGGGGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTgg                                                     <  1:3457/100‑1 (MQ=255)

                                                 ccTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTgggg                                                   <  2:76260/100‑1 (MQ=255)

                                                  cTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggggg                                                  >  2:1431173/1‑100 (MQ=255)

                                                  cTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggggg                                                  >  2:821146/1‑100 (MQ=255)

                                                      gACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCgg                                              <  2:4384725/100‑1 (MQ=255)

                                                        cATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGAt                                            >  1:390846/1‑100 (MQ=255)

                                                          tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATCAACCGGTAGAAGCTCTGGGGGCCGGATTa                                          >  2:2646945/1‑100 (MQ=255)

                                                          tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTa                                          >  2:3934034/1‑100 (MQ=255)

                                                          tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTa                                          >  1:1157233/1‑100 (MQ=255)

                                                           atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTAc                                         >  2:4049082/1‑100 (MQ=255)

                                                           atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTAc                                         <  1:2462443/100‑1 (MQ=255)

                                                                 ccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatg                                   <  1:1442859/100‑1 (MQ=255)

                                                                 ccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatg                                   >  2:1696549/1‑100 (MQ=255)

                                                                  cacaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatga                                  >  1:475395/1‑100 (MQ=255)

                                                                  cacaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatga                                  <  2:4247349/100‑1 (MQ=255)

                                                                   acaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatgat                                 >  2:2332010/1‑100 (MQ=255)

                                                                   acaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatgat                                 <  2:218840/100‑1 (MQ=255)

                                                                     aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGAtgatga                               >  1:1094003/1‑100 (MQ=255)

                                                                     aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGAtgatga                               <  2:4300597/100‑1 (MQ=255)

                                                                      aGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGaa                              <  1:1333960/100‑1 (MQ=255)

                                                                          gCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcgg                          >  2:3567770/1‑100 (MQ=255)

                                                                             cTGGCAACACATTTCTCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggcgg                       <  2:4331261/100‑1 (MQ=255)

                                                                                gCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTA‑CACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAt                   <  1:1424419/100‑1 (MQ=255)

                                                                                gCGACACATTTCACGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAGACCGGTAAAAGCCCTGGGGGCCGGATTACTGGATGATGAAAcggcgg                       <  1:833353/97‑1 (MQ=255)

                                                                                 cGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAt                   >  2:454908/1‑100 (MQ=255)

                                                                                 cGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAt                   <  2:1812661/100‑1 (MQ=255)

                                                                                  gACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAtt                  >  2:2806623/1‑100 (MQ=255)

                                                                                      caTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc              >  2:1036250/1‑100 (MQ=255)

                                                                                      aaacaaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc              >  2:3811223/7‑100 (MQ=255)

                                                                                       aTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTctct             <  1:4223311/100‑1 (MQ=255)

                                                                                         ttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg           <  1:2719987/100‑1 (MQ=255)

                                                                                         caaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg           >  2:1445503/4‑100 (MQ=255)

                                                                                          tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg          >  1:1052347/1‑100 (MQ=255)

                                                                                          aaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg          >  1:1120200/3‑100 (MQ=255)

                                                                                           cccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc         >  2:2103584/1‑100 (MQ=255)

                                                                                             cGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCAc       <  1:1184964/100‑1 (MQ=255)

                                                                                              ggCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACt      >  2:341911/1‑100 (MQ=255)

                                                                                               gcagcaCCGCCGGTTGTTCGTACCACCAGTGATAAACCGGTAGACGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCCGTAATTTCTCTTCGCACtt     <  1:3322289/100‑1 (MQ=255)

                                                                                                 aGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAGCGGCGGTAATTTCTCTTGGCACTTa    <  1:1290840/99‑1 (MQ=255)

TTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACA  >  NC_000913/2902395‑2902592

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