BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	707,634	A→G	A26A (GCA→GCG)	ybfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	707,634	0	A	G	98.2%	530.9 / ‑1.0	171	A26A (GCA→GCG)	ybfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (86/81); total (87/84)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

Only 100 of 179 total aligned reads displayed.

AGAGGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCAATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGA  >  NC_000913/707544‑707723

agagGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGAt                                                                                         <  1:2947998/93‑1 (MQ=255)

agagGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGAt                                                                                         <  1:390150/93‑1 (MQ=255)

agagGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGAt                                                                                         >  2:4106888/1‑93 (MQ=255)

  agGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGAGATGc                                                                                       >  1:4367658/1‑93 (MQ=255)

   ggATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGcc                                                                                      >  1:1198776/1‑93 (MQ=255)

   ggATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGAGATGcc                                                                                      >  2:1514647/1‑93 (MQ=255)

      ttAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCggg                                                                                   >  1:3781451/1‑92 (MQ=255)

      ttAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGa                                                                                   <  2:4587110/93‑1 (MQ=255)

      ttAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGa                                                                                   >  2:2749181/1‑93 (MQ=255)

       tAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAg                                                                                  <  1:2921278/93‑1 (MQ=255)

        aaCCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAgg                                                                                 <  1:3388838/93‑1 (MQ=255)

         aCCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGc                                                                                >  1:690445/1‑93 (MQ=255)

          cccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGGGGCt                                                                               >  2:508056/1‑93 (MQ=255)

          cccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCt                                                                               <  2:3654253/93‑1 (MQ=255)

           ccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTc                                                                              <  2:3674106/93‑1 (MQ=255)

           ccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTc                                                                              >  2:2426186/1‑93 (MQ=255)

           ccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTc                                                                              >  2:89573/1‑93 (MQ=255)

            cATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCAATGCCGGAGGCTcg                                                                             >  2:3328831/1‑93 (MQ=255)

               gCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCgcgc                                                                          >  2:3929092/1‑93 (MQ=255)

                cGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGcc                                                                         <  2:2542556/93‑1 (MQ=255)

                 gTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCg                                                                        >  1:251749/1‑93 (MQ=255)

                  tACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGa                                                                       >  2:1825880/1‑93 (MQ=255)

                  tACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCAATGCCGGAGGCTCGCGCCGa                                                                       <  2:3805987/93‑1 (MQ=255)

                  tACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATCGAGACGCCGGCGGCTCGCGCCGa                                                                       >  2:4074969/1‑93 (MQ=255)

                  tACGTTTAGTGGCAAACGCTGTACGCTGGCGCTGGCTATCGCCGGTTTTACTGCAATGACGGGCTCTATGGCGATGCTGGTGGCTCGCGTCGa                                                                       >  2:4338718/1‑93 (MQ=255)

                    cGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAg                                                                     <  1:1046096/93‑1 (MQ=255)

                    cGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAg                                                                     >  1:2790147/1‑93 (MQ=255)

                    cGTTTAGTGGCAAACGTAGAACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGCGGCTCGCGCCGTAg                                                                     >  2:3930679/1‑93 (MQ=255)

                     gTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAgg                                                                    <  2:1914217/93‑1 (MQ=255)

                      tttAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGGGATCCCGGAGGCTCGGGCCGAAGGa                                                                   >  1:1633800/1‑93 (MQ=255)

                      tttAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGa                                                                   <  2:1429195/93‑1 (MQ=255)

                      tttAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGa                                                                   >  2:4694120/1‑93 (MQ=255)

                      ttgaGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGa                                                                   <  2:326203/90‑1 (MQ=255)

                        tAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGAtt                                                                 <  1:3198369/93‑1 (MQ=255)

                         aGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTc                                                                <  1:2664801/93‑1 (MQ=255)

                         aGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTc                                                                <  2:1391059/93‑1 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGAttga                                                               >  2:3743594/1‑91 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               <  1:4389516/93‑1 (MQ=255)

                            ggCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGGTCTCGGCGAAGGATTCATc                                                             <  2:790512/93‑1 (MQ=255)

                             gCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcg                                                            >  2:628988/1‑93 (MQ=255)

                               aaaCGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcgac                                                          >  2:1073683/1‑93 (MQ=255)

                                aaCGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcgacg                                                         >  1:597416/1‑93 (MQ=255)

                                 aCGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcgacga                                                        >  2:1630439/1‑93 (MQ=255)

                                    tagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTc                                                     >  1:1797449/1‑93 (MQ=255)

                                     agtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGAGGATTCa                                                    >  2:4186864/1‑93 (MQ=255)

                                      gtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCaa                                                   >  1:753183/1‑93 (MQ=255)

                                       taCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAAc                                                  <  1:4725335/93‑1 (MQ=255)

                                        aCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAAcc                                                 <  2:3742814/93‑1 (MQ=255)

                                        aCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAAcc                                                 <  2:2386903/93‑1 (MQ=255)

                                         cgctggcgctggCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCt                                                <  2:1303742/93‑1 (MQ=255)

                                          gctggctctggcTATCGCCGGTGTTACAGCAATTTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGTTTCAACCtt                                               <  1:841173/93‑1 (MQ=255)

                                          gctggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCtt                                               >  2:1476129/1‑93 (MQ=255)

                                           ctggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTa                                              <  2:319428/93‑1 (MQ=255)

                                           ctggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGAGATGCCGGAGGCTCGCGCCGAAGGATTCATTCGCGATTCAACCTTa                                              >  1:1913614/1‑93 (MQ=255)

                                            tggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTaa                                             >  2:2915900/1‑93 (MQ=255)

                                            tggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTaa                                             <  1:1448336/93‑1 (MQ=255)

                                            tggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTaa                                             >  1:1243101/1‑93 (MQ=255)

                                               cgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcg                                          <  2:1364083/93‑1 (MQ=255)

                                                gctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcgg                                         >  1:3462265/1‑93 (MQ=255)

                                                   ggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGTAGGATTCATCGACGATTCAACCTTAACcggcgg                                      >  1:2072556/1‑93 (MQ=255)

                                                   ggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcggcgg                                      <  2:2820470/93‑1 (MQ=255)

                                                   ggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcggcgg                                      <  1:1039916/93‑1 (MQ=255)

                                                   ggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcggcgg                                      >  1:3345322/1‑93 (MQ=255)

                                                    gcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGt                                     <  1:2346026/93‑1 (MQ=255)

                                                     cTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTa                                    <  1:2152724/93‑1 (MQ=255)

                                                      tATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTAt                                   >  2:3035185/1‑93 (MQ=255)

                                                       aTCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATc                                  <  1:1483975/93‑1 (MQ=255)

                                                         cGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTa                                <  2:756058/93‑1 (MQ=255)

                                                         cGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTa                                <  1:3544325/93‑1 (MQ=255)

                                                          gCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGCTTCAACCTTAACCGGCGGTATCTAt                               >  1:1601737/1‑93 (MQ=255)

                                                           ccGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTAtt                              <  1:2876709/93‑1 (MQ=255)

                                                           ccGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTAtt                              <  1:2364375/93‑1 (MQ=255)

                                                           ccGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTAtt                              >  1:2657885/1‑93 (MQ=255)

                                                            cGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTa                             <  2:3579784/93‑1 (MQ=255)

                                                            cGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTa                             <  2:837009/93‑1 (MQ=255)

                                                             ggTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTAc                            >  1:3103245/1‑93 (MQ=255)

                                                             ggTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTAc                            >  2:2205395/1‑93 (MQ=255)

                                                             ggTGTTACAGCAATGTCGGCCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTAc                            <  1:4616085/93‑1 (MQ=255)

                                                               tgtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTg                          <  1:2200833/93‑1 (MQ=255)

                                                                gtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTgg                         >  1:381050/1‑93 (MQ=255)

                                                                 ttACAGCAATGTCGGGCTTTATTGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGCTTCAACCTTAACCGGCGGTATCTATTACTGGc                        >  2:2832476/1‑93 (MQ=255)

                                                                   aCAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGGTTCAACCTTAACCGGCGGTATCTATTACTGGCAg                      >  1:167517/1‑93 (MQ=255)

                                                                      gCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGt                   >  1:2231657/1‑93 (MQ=255)

                                                                       cAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTg                  <  1:3665626/93‑1 (MQ=255)

                                                                        aaTGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGa                 >  1:1967092/1‑93 (MQ=255)

                                                                         aTGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGaa                <  1:4158864/93‑1 (MQ=255)

                                                                         aTGTAGGACTATATGGCGATACCGGAGGCTCGCGACGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGaa                <  2:4309906/93‑1 (MQ=255)

                                                                           gTCGGGCTTTATGGCGATGCCGGAGGCTTGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAAcg              <  1:3152033/93‑1 (MQ=255)

                                                                            tCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAAcgc             <  1:4288674/93‑1 (MQ=255)

                                                                            tCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAAcgc             >  1:228707/1‑93 (MQ=255)

                                                                            tCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAAcgc             <  2:1737143/93‑1 (MQ=255)

                                                                             cGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAAcgcg            >  2:2087912/1‑93 (MQ=255)

                                                                                gCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGAcc         <  2:1296535/93‑1 (MQ=255)

                                                                                  tttATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACAGt       >  2:2144474/1‑93 (MQ=255)

                                                                                    tATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTTTCTATTACTGGCAGCGTGAACGCGACCGtca     >  1:2172267/1‑91 (MQ=255)

                                                                                    tATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTaa     <  2:2530743/93‑1 (MQ=255)

                                                                                      tGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAg   >  1:3923108/1‑93 (MQ=255)

                                                                                      tGGCAACCCCGAAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAg   <  2:3118578/93‑1 (MQ=255)

                                                                                       ggCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGa  <  1:1791304/93‑1 (MQ=255)

                                                                                       ggCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGa  >  1:906192/1‑93 (MQ=255)

AGAGGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCAATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGA  >  NC_000913/707544‑707723

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