BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,246,351	T→C	S68G (AGT→GGT)	lysP ←	lysine‑specific permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,246,351	0	T	C	99.6%	856.6 / ‑4.8	265	S68G (AGT→GGT)	lysP	lysine‑specific permease
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (139/125); total (139/126)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

Only 100 of 267 total aligned reads displayed.

GCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGG  >  NC_000913/2246263‑2246440

  gAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcat                                                                                     <  1:1139612/93‑1 (MQ=255)

  gAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcat                                                                                     >  2:5691333/1‑93 (MQ=255)

   aaGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcatc                                                                                    >  1:928124/1‑93 (MQ=255)

    aGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcatca                                                                                   <  1:145539/93‑1 (MQ=255)

    aGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGtcatca                                                                                   >  2:1991776/1‑93 (MQ=255)

      ccAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAgg                                                                                 >  2:5029804/1‑93 (MQ=255)

       cAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGa                                                                                >  1:1215467/1‑93 (MQ=255)

       cAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGa                                                                                >  2:6637121/1‑93 (MQ=255)

       cAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGa                                                                                >  2:4581483/1‑93 (MQ=255)

        aaaGCCTTTTTCAACATAGTTCTGACCGTAAGTGGCAAACGGACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGaa                                                                               >  1:6512295/1‑93 (MQ=255)

        aaaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGaa                                                                               <  2:3114783/93‑1 (MQ=255)

        aaaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGaa                                                                               >  1:2366017/1‑93 (MQ=255)

         aaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAg                                                                              >  2:1031041/1‑93 (MQ=255)

         aaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAg                                                                              >  2:380609/1‑93 (MQ=255)

          aGCCTTCTTCAACATAGTTCTGACCGTAAGTGACAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGt                                                                             >  1:772523/1‑93 (MQ=255)

            ccTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTaa                                                                           >  1:1020248/1‑93 (MQ=255)

             cttcttCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTaaa                                                                          >  2:1812223/1‑93 (MQ=255)

             cttcttCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTaaa                                                                          >  2:2332157/1‑93 (MQ=255)

             cttcttCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTaaa                                                                          >  1:3368705/1‑93 (MQ=255)

              ttcttcAACATAGTTCTGTCCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAc                                                                         <  2:5526945/93‑1 (MQ=255)

              ttcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAc                                                                         <  1:986279/93‑1 (MQ=255)

              ttcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAc                                                                         >  2:860795/1‑93 (MQ=255)

               tcttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAcc                                                                        <  2:6137823/93‑1 (MQ=255)

                cttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAAcaa                                                                       >  2:3058449/1‑91 (MQ=255)

                 ttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCAt                                                                      >  1:3914690/1‑93 (MQ=255)

                  tcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATc                                                                     <  2:3316188/93‑1 (MQ=255)

                  tcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATc                                                                     <  1:1798743/93‑1 (MQ=255)

                   cAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCa                                                                    <  2:2422990/93‑1 (MQ=255)

                   cAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCa                                                                    <  2:5283955/93‑1 (MQ=255)

                    aaCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAg                                                                   >  1:5019587/1‑93 (MQ=255)

                    aaCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAg                                                                   >  1:3711090/1‑93 (MQ=255)

                     aCATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAgg                                                                  <  2:6390139/93‑1 (MQ=255)

                        tAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCa                                                               >  1:1700503/1‑93 (MQ=255)

                        tAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCa                                                               >  2:4680144/1‑93 (MQ=255)

                          gTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAAt                                                             <  1:5032994/93‑1 (MQ=255)

                           ttCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATc                                                            >  1:1339582/1‑93 (MQ=255)

                           ttCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATc                                                            >  2:1385543/1‑93 (MQ=255)

                             cTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAg                                                          <  1:2692334/93‑1 (MQ=255)

                               gACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCa                                                        <  1:6408444/93‑1 (MQ=255)

                               gACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCa                                                        >  2:1598893/1‑93 (MQ=255)

                                aCCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCat                                                       <  2:6061890/93‑1 (MQ=255)

                                aCCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCat                                                       <  1:1288072/93‑1 (MQ=255)

                                  cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCatat                                                     >  2:1701533/1‑93 (MQ=255)

                                   gTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCAtata                                                    >  1:3003843/1‑93 (MQ=255)

                                    tAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATAc                                                   <  1:2957995/93‑1 (MQ=255)

                                    tAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATAc                                                   <  1:1305138/93‑1 (MQ=255)

                                      aGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACga                                                 <  1:5744430/93‑1 (MQ=255)

                                      aGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACga                                                 <  1:6578334/93‑1 (MQ=255)

                                        tGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACgaga                                               <  2:3030675/93‑1 (MQ=255)

                                         ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGagag                                              <  2:128175/93‑1 (MQ=255)

                                         ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGagag                                              <  2:75581/93‑1 (MQ=255)

                                         ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGagag                                              >  2:2690979/1‑93 (MQ=255)

                                          gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGc                                             <  2:5438108/93‑1 (MQ=255)

                                          gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGc                                             <  1:3624431/93‑1 (MQ=255)

                                          gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGc                                             <  1:4665887/93‑1 (MQ=255)

                                              aCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATg                                         >  1:2910775/1‑93 (MQ=255)

                                                gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGcc                                       <  2:3870250/93‑1 (MQ=255)

                                                gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGcc                                       <  1:4059686/93‑1 (MQ=255)

                                                 aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGccc                                      <  1:2201363/93‑1 (MQ=255)

                                                 aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGccc                                      >  1:4749809/1‑93 (MQ=255)

                                                    cGGAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgc                                   >  1:2528598/1‑93 (MQ=255)

                                                     ggAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgcc                                  >  2:3517346/1‑93 (MQ=255)

                                                     ggAAACCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgcc                                  >  1:1942/1‑93 (MQ=255)

                                                        aaCCGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCgg                               >  2:4166919/1‑93 (MQ=255)

                                                          ccGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAc                             >  2:2152150/1‑93 (MQ=255)

                                                          ccGGCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAc                             <  2:4133575/93‑1 (MQ=255)

                                                            ggCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcct                           >  2:3212833/1‑93 (MQ=255)

                                                             gCATATATGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctg                          >  2:6455/1‑93 (MQ=255)

                                                               atatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcc                        >  1:2868775/1‑93 (MQ=255)

                                                               atatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATAAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcc                        <  2:4322105/93‑1 (MQ=255)

                                                                tatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcct                       <  2:6152738/93‑1 (MQ=255)

                                                                tatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcct                       <  2:3688112/93‑1 (MQ=255)

                                                                tatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcct                       >  1:6082499/1‑93 (MQ=255)

                                                                tatatGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAcctgcct                       >  1:5400008/1‑93 (MQ=255)

                                                                   atGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTgag                    <  2:5749514/93‑1 (MQ=255)

                                                                    tGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTgaga                   <  2:2086351/93‑1 (MQ=255)

                                                                    tGCAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTgaga                   >  2:4467762/1‑93 (MQ=255)

                                                                      cAGCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCACCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGaaa                 <  2:1768889/93‑1 (MQ=255)

                                                                        gCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGGCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATc               >  1:5236528/1‑93 (MQ=255)

                                                                        gCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATc               >  1:558991/1‑93 (MQ=255)

                                                                        gCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATc               >  1:4704193/1‑93 (MQ=255)

                                                                        gCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATc               >  1:4089536/1‑93 (MQ=255)

                                                                        gCCAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATc               <  2:2242506/93‑1 (MQ=255)

                                                                          cAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGt             <  1:5543482/93‑1 (MQ=255)

                                                                          cAGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGt             <  2:5777145/93‑1 (MQ=255)

                                                                          cAGTCCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGt             <  1:5291363/93‑1 (MQ=255)

                                                                           aGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGtt            <  1:3484482/93‑1 (MQ=255)

                                                                           aGTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGtt            <  1:3220013/93‑1 (MQ=255)

                                                                            gTTCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTg           >  1:207635/1‑93 (MQ=255)

                                                                              tCACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgcc         >  1:6099417/1‑92 (MQ=255)

                                                                               cACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgcgc        >  1:3021030/1‑93 (MQ=255)

                                                                               cACCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgcgc        <  2:3427119/93‑1 (MQ=255)

                                                                                aCCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGcc       >  1:4382187/1‑93 (MQ=255)

                                                                                aCCGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGcc       >  2:652330/1‑93 (MQ=255)

                                                                                 ccGAGACCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCa      >  2:5903256/1‑93 (MQ=255)

                                                                                   gagaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCaga    >  2:6626450/1‑93 (MQ=255)

                                                                                    agaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCagag   >  1:5102760/1‑93 (MQ=255)

                                                                                    agaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCagag   >  1:3303047/1‑93 (MQ=255)

                                                                                     gaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAgg  >  2:4941185/1‑93 (MQ=255)

                                                                                     gaCCGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAgg  >  2:2816069/1‑93 (MQ=255)

GCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGG  >  NC_000913/2246263‑2246440

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