BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,902,495	A→C	17.1%	G200G (GGA→GGC) ‡	ygcE →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,902,495	0	A	C	17.1%	430.7 / 92.6	188	G200G (GGA→GGC)	ygcE	putative kinase
Reads supporting (aligned to +/- strand): ref base A (81/75); new base C (7/25); total (88/100)
Fisher's exact test for biased strand distribution p-value = 1.87e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-02

Only 100 of 190 total aligned reads displayed.

AATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATC  >  NC_000913/2902397‑2902594

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

  tatCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTg                                                                                                  >  1:1409283/1‑100 (MQ=255)

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

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

      ccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCgg                                                                                              >  2:4376991/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)

        ccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGTCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgtt                                                                                            <  1:2345369/100‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

           tCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaat                                                                                         <  2:34711/100‑6 (MQ=255)

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

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

            cATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatg                                                                                        <  1:1350467/100‑7 (MQ=255)

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

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

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

             atatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCAGttgtttg                                                                                       <  2:1958865/100‑1 (MQ=255)

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

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

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

                 gCTGTTGGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatg                                                                                        <  1:3457901/95‑7 (MQ=255)

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

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

                      ttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAg                                                                              >  2:4322941/1‑100 (MQ=255)

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

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

                          tgtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                          >  2:2997365/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                                                                     <  1:2383674/100‑1 (MQ=255)

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

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

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

                                  tGCTTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCgg                                                                   <  2:1819173/99‑1 (MQ=255)

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

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

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

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

                                     cTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGa                                                               >  2:2866404/1‑100 (MQ=255)

                                     cTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGa                                                               <  2:2325220/100‑1 (MQ=255)

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

                                        gCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGc                                                            >  1:3451947/1‑100 (MQ=255)

                                           ccGTCCTCGGAGATATCACACCACAAGCCGCACTGGCGACACATTGCCCGGACGGACTCCGGGTGGTTTGTACCACAAGTGATAACCCGGTAGAAg                                                             >  1:3607907/1‑96 (MQ=255)

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

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

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

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

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

                                                  cGGACATATCACACCACAAGCCGCACTGGCGACACAGTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggggg                                                    <  2:4309114/98‑1 (MQ=255)

                                                   ggACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCg                                                 >  1:2421911/1‑100 (MQ=255)

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

                                                      cATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGAt                                              <  2:2190008/100‑1 (MQ=255)

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

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

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

                                                        tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTa                                            <  2:3417553/100‑1 (MQ=255)

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

                                                             caccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGga                                       >  2:2061052/1‑100 (MQ=255)

                                                              accacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgat                                      <  2:2818586/100‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                           cTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggcgg                         >  2:2576673/1‑100 (MQ=255)

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

                                                                               cGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAt                     >  2:454908/1‑100 (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:4360443/1‑100 (MQ=255)

                                                                                       ttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg             >  2:1728/1‑100 (MQ=255)

                                                                                       ttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg             >  1:497936/1‑100 (MQ=255)

                                                                                       caaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGACACGGCGATAATTTCTCtt              >  1:4162429/4‑99 (MQ=255)

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

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

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

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

                                                                                         accGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc           >  2:3145507/2‑100 (MQ=255)

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

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

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

                                                                                                ggACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACa    >  2:2104439/1‑100 (MQ=255)

                                                                                                ggACTGCCGGGTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACa    >  2:2396251/1‑100 (MQ=255)

AATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATC  >  NC_000913/2902397‑2902594

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