BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	715,849	A→G	63.6%	L146L (TTA→TTG)	ybfP →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	715,849	0	A	G	63.6%	127.4 / 129.3	132	L146L (TTA→TTG)	ybfP	lipoprotein
Reads supporting (aligned to +/- strand): ref base A (14/34); new base G (46/38); total (60/72)
Fisher's exact test for biased strand distribution p-value = 6.22e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-02

Only 100 of 140 total aligned reads displayed.

GCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATG  >  NC_000913/715752‑715944

gCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcg                                                                                               >  1:216078/1‑100 (MQ=255)

   aGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTacgacg                                                                                            >  1:2759340/1‑100 (MQ=255)

     ggCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCt                                                                                          >  1:4107038/1‑100 (MQ=255)

      gCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTg                                                                                         <  1:521051/100‑1 (MQ=255)

       cAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGc                                                                                        <  2:368621/100‑1 (MQ=255)

        atctccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGcc                                                                                       <  1:1362130/96‑1 (MQ=255)

        aaGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGACGACGCATGTCGGGATTATTGAGTTGCGACGCTGcc                                                                                       >  2:3280880/1‑100 (MQ=255)

        aaGCCCAGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGTATTTTTTGGCTCGATGATGCATGTCGGGATTATTGGGTTGCGACGCTGcc                                                                                       <  1:320879/100‑1 (MQ=255)

          gCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAg                                                                                     <  2:3163419/100‑1 (MQ=255)

          gCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAg                                                                                     <  2:30453/100‑1 (MQ=255)

          gCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAg                                                                                     <  2:2813928/100‑1 (MQ=255)

           cccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGGTTATTGAGTTGCGACGCTGCCAGt                                                                                    >  1:1227430/1‑100 (MQ=255)

           cccGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGt                                                                                    >  2:81939/1‑100 (MQ=255)

            ccGGCACGCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGAGGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTc                                                                                   <  2:4190730/100‑1 (MQ=255)

             cGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCa                                                                                  <  1:1580225/100‑1 (MQ=255)

              ggCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGAGTATTGAGTTGCGACGCTGCCAGTCAg                                                                                 >  1:672337/1‑100 (MQ=255)

                 aCTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGaaa                                                                              <  2:2835202/100‑1 (MQ=255)

                 aCTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGaaa                                                                              <  2:2390983/100‑1 (MQ=255)

                 aCTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCAGGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGaaa                                                                              <  2:1458655/100‑1 (MQ=255)

                  cTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGaaaa                                                                             >  2:4232893/1‑100 (MQ=255)

                  cTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGaaaa                                                                             >  2:349/1‑100 (MQ=255)

                  cTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTGTGGCTCTATGATGCATGTCGGGATTATTGATTTGCGAAGCTGCCAGTAAGaaaa                                                                             >  1:884835/1‑100 (MQ=255)

                   tCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAAt                                                                            <  2:2480366/100‑1 (MQ=255)

                      aCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGc                                                                         >  2:2158633/1‑100 (MQ=255)

                        cGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGcc                                                                        >  2:2986645/1‑99 (MQ=255)

                        cGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTCACGACGCTGCCAGTCAGAAAATTGCCa                                                                       <  1:2613730/100‑1 (MQ=255)

                         gATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCa                                                                       <  1:1101100/99‑1 (MQ=255)

                         gATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAg                                                                      <  1:1901984/100‑1 (MQ=255)

                          aTAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGTCAGTCAGAAAATTGCCAGc                                                                     <  1:3760462/100‑1 (MQ=255)

                            aGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCaa                                                                   <  2:1712990/100‑1 (MQ=255)

                            aGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCaa                                                                   <  2:492825/100‑1 (MQ=255)

                             gTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAg                                                                  <  2:765346/100‑1 (MQ=255)

                             gTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAg                                                                  >  2:739406/1‑100 (MQ=255)

                                  tGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTg                                                             >  2:2721071/1‑100 (MQ=255)

                                   gAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGa                                                            <  2:3866521/100‑1 (MQ=255)

                                    aGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaa                                                           <  1:1731660/100‑1 (MQ=255)

                                     ggATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa                                                          <  2:1389523/100‑1 (MQ=255)

                                      gATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAg                                                         <  2:3604305/100‑1 (MQ=255)

                                       aTTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGc                                                        <  1:974344/100‑1 (MQ=255)

                                       aTTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGc                                                        <  1:97156/100‑1 (MQ=255)

                                       aTTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGc                                                        <  1:2382297/100‑1 (MQ=255)

                                          aTGTCAGCGAAAGGGGATTTTTTGGCTCGATGATGCATATCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                     >  1:2723642/1‑100 (MQ=255)

                                          aTGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                     >  1:2212776/1‑100 (MQ=255)

                                          aTGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                     <  1:1920127/100‑1 (MQ=255)

                                            gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGCTTATTGAGTTACGACGCTGCCCGTCCGACCATCGCCCGCacgcatcg                                                             >  1:3548782/1‑83 (MQ=255)

                                            gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaa                                                   >  2:2557734/1‑100 (MQ=255)

                                              cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                 >  2:1486679/1‑100 (MQ=255)

                                              cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                 >  2:3832290/1‑100 (MQ=255)

                                              cAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt                                                 <  2:3953619/100‑1 (MQ=255)

                                               aGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATa                                                >  2:1358255/1‑100 (MQ=255)

                                                 cGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca                                              <  1:3412134/100‑1 (MQ=255)

                                                 cGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca                                              <  1:674087/100‑1 (MQ=255)

                                                  gAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac                                             >  2:2925112/1‑100 (MQ=255)

                                                   aaaGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca                                            <  1:1032049/100‑1 (MQ=255)

                                                   aaaGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca                                            >  2:1119091/1‑100 (MQ=255)

                                                    aaGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAt                                           <  1:903630/100‑1 (MQ=255)

                                                      gCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTa                                         >  1:4038087/1‑100 (MQ=255)

                                                       cGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaa                                        <  1:1249443/100‑1 (MQ=255)

                                                       cGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaa                                        <  2:1150272/100‑1 (MQ=255)

                                                        ggATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaaa                                       >  2:319695/1‑100 (MQ=255)

                                                        ggATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaaa                                       >  1:2130060/1‑100 (MQ=255)

                                                          aTTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAAt                                     >  2:3892296/1‑100 (MQ=255)

                                                            tttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                   <  2:1667103/100‑1 (MQ=255)

                                                            tttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                   <  1:2426981/100‑1 (MQ=255)

                                                             ttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                   >  1:2863841/1‑99 (MQ=255)

                                                              tttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATGAATACACATTAAAATAGTg                                 <  2:1904967/100‑1 (MQ=255)

                                                                tGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAg                               <  1:4165967/100‑1 (MQ=255)

                                                                 ggCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGt                              >  2:574519/1‑100 (MQ=255)

                                                                  gCTCGATGATGCATGTCGGGATTATTGAGTTGGGACGATGTCAGTGAGAGAATTGAAGGCAAGAATTGAAAGCGATAAATGGACATTAAAATAGTggg                               >  1:630477/1‑96 (MQ=255)

                                                                  gCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTc                             >  1:4195425/1‑100 (MQ=255)

                                                                  gCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTc                             >  2:2553324/1‑100 (MQ=255)

                                                                  gCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTc                             >  1:4356879/1‑100 (MQ=255)

                                                                    tCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAgatcggaag                           >  1:3920013/1‑92 (MQ=255)

                                                                    tCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCg                           >  1:1044281/1‑100 (MQ=255)

                                                                    tCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCg                           >  2:78877/1‑100 (MQ=255)

                                                                     cGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCgg                          <  1:651738/100‑1 (MQ=255)

                                                                      gatgatGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCggg                         >  1:3544546/1‑100 (MQ=255)

                                                                      gatgatGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCggg                         <  1:1648490/100‑1 (MQ=255)

                                                                      gatgatGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCgg                          >  1:3165214/1‑99 (MQ=255)

                                                                        tgatgCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGaa                       <  1:1254314/100‑1 (MQ=255)

                                                                         gatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAAt                      <  1:1359352/100‑1 (MQ=255)

                                                                         gatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAAt                      <  2:4131904/100‑1 (MQ=255)

                                                                          atgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATg                     <  1:2623341/100‑1 (MQ=255)

                                                                          atgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATg                     <  2:3322584/100‑1 (MQ=255)

                                                                          atgCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATg                     <  1:793209/100‑1 (MQ=255)

                                                                               tGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAcc                >  2:664242/1‑100 (MQ=255)

                                                                                gTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAccc               >  2:2577029/1‑100 (MQ=255)

                                                                                gTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAccc               <  2:1919592/100‑1 (MQ=255)

                                                                                gTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAccc               >  1:2995191/1‑100 (MQ=255)

                                                                                  cGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCgg             >  1:1246429/1‑100 (MQ=255)

                                                                                    ggATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAAACCgg             >  1:405062/1‑98 (MQ=255)

                                                                                     gATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGAcc          <  2:216078/100‑1 (MQ=255)

                                                                                        tattGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGa       >  1:3931612/1‑100 (MQ=255)

                                                                                         attGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGaa      >  2:2383447/1‑100 (MQ=255)

                                                                                         attGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGaa      <  2:1764732/100‑1 (MQ=255)

                                                                                           tGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACa    <  2:4038117/100‑1 (MQ=255)

                                                                                           tGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACa    <  2:2397790/100‑1 (MQ=255)

                                                                                           tGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACa    <  1:1041375/100‑1 (MQ=255)

                                                                                            gAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACAt   <  1:442094/100‑1 (MQ=255)

                                                                                              gTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAGAAATACACATTAAAATAGTGAGTCCGGGAGTGAAACCCGGACCTGAACATg  >  1:1629001/1‑99 (MQ=255)

GCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACATG  >  NC_000913/715752‑715944

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