BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,413,167	A→C	E47A (GAG→GCG)	yjfC →	ATP‑Grasp family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,413,167	0	A	C	97.6%	732.4 / ‑5.0	212	E47A (GAG→GCG)	yjfC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base A (1/4); new base C (119/87); total (120/92)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

ATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCAACTGGCGA  >  NC_000913/4413070‑4413262

 tCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCa                                                                                              <  1:4399320/100‑1 (MQ=255)

 tCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCa                                                                                              <  1:669574/100‑1 (MQ=255)

  cgccgcTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAg                                                                                             <  2:1509312/100‑1 (MQ=255)

    ccgcTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGAt                                                                                           >  1:4313502/1‑100 (MQ=255)

    ccgcTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGAt                                                                                           >  1:495802/1‑100 (MQ=255)

    ccgcTGACAACGGTTTCGACTTACATATCAACGACAATGAAAACTAGTGGGATGAGAGTCTGGCTTAACGATGTACTCTGCCCCTGATTGAAGCGCAGAt                                                                                           <  2:3977160/100‑1 (MQ=255)

     cgcTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATc                                                                                          <  1:764501/100‑1 (MQ=255)

     cgcTGACAACGGTGTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATc                                                                                          <  2:4307878/100‑1 (MQ=255)

      gcTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCg                                                                                         >  2:3116794/1‑100 (MQ=255)

       cTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGa                                                                                        >  1:2320726/1‑100 (MQ=255)

        tGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATGGGGATGGGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGaa                                                                                       <  1:2655749/100‑1 (MQ=255)

         gACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGaaa                                                                                      >  2:2390341/1‑100 (MQ=255)

           cAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGaaaaa                                                                                    >  2:3048275/1‑100 (MQ=255)

            aaCGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCGGATCGAAAAAc                                                                                   >  2:3612003/1‑100 (MQ=255)

              cGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCg                                                                                 <  2:741475/100‑1 (MQ=255)

              cGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCg                                                                                 >  1:35539/1‑100 (MQ=255)

              cGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCg                                                                                 >  1:594377/1‑100 (MQ=255)

                gTTTCGACTTTCATAAAATCGACAATGAAATCTATTGGGATGAGAGCCGGGCTTACAGCTTTACTCTGCGGCAGATTGAAGCGCAGATTGAAAAACCGAc                                                                               >  2:4995041/1‑100 (MQ=255)

                     gACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCgg                                                                          <  1:1258057/100‑1 (MQ=255)

                      aCTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGa                                                                         <  2:1404007/100‑1 (MQ=255)

                       cTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGaa                                                                        >  2:1847956/1‑100 (MQ=255)

                       cTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGaa                                                                        >  2:857302/1‑100 (MQ=255)

                         ttCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACt                                                                      >  1:3897140/1‑100 (MQ=255)

                            atatCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCa                                                                   <  2:959505/100‑1 (MQ=255)

                            atatCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCa                                                                   <  2:4727499/100‑1 (MQ=255)

                             tatCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCGGATCGAAAAACCGACTGCGGAACTGCAt                                                                  >  2:4767093/1‑100 (MQ=255)

                             tatCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCAt                                                                  >  1:429529/1‑100 (MQ=255)

                               tcatcGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCa                                                                <  1:3956382/100‑1 (MQ=255)

                               tcatcGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCa                                                                >  1:5204135/1‑100 (MQ=255)

                                catcGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAg                                                               >  1:2446286/1‑100 (MQ=255)

                                catcGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAg                                                               <  1:1742468/100‑1 (MQ=255)

                                   cGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGCTTTACTCTGCGCCAGATAGAAGCGCAGATCGAAAAAGCGACTGCGGAACTCCATCAGAtg                                                            >  1:4856552/1‑100 (MQ=255)

                                    gACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGAtgt                                                           >  1:1143746/1‑100 (MQ=255)

                                      cAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGc                                                         >  2:4263828/1‑100 (MQ=255)

                                       aaTGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGcc                                                        <  1:3602815/100‑1 (MQ=255)

                                       aaTGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGcc                                                        <  2:5393280/100‑1 (MQ=255)

                                        aTGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCt                                                       >  1:3703594/1‑100 (MQ=255)

                                         tGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGAAGATAGAAAAACCGACTGCGGGACTGCGTCAGATGTGccga                                                      >  2:3176101/1‑98 (MQ=255)

                                           aaaTCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGa                                                    <  1:4484856/100‑1 (MQ=255)

                                           aaaTCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGa                                                    >  2:618030/1‑100 (MQ=255)

                                            aaTCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAg                                                   >  1:2591469/1‑100 (MQ=255)

                                            aaTCTATTGGGATGAGAGTCGGGCTGACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAg                                                   <  2:3732266/100‑1 (MQ=255)

                                             aTCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAgg                                                  <  2:4468985/100‑1 (MQ=255)

                                             aTCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAgg                                                  <  1:5527841/100‑1 (MQ=255)

                                             aTCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAgg                                                  >  1:1984908/1‑100 (MQ=255)

                                                tATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAggtgg                                               <  2:4808551/100‑1 (MQ=255)

                                                 aTTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAggtggt                                              <  2:379155/100‑1 (MQ=255)

                                                 aTTGAGATGAGAGTCGAGCTTACGGTTGTACTCTGTGTCTGATTTAAGCGTAGATCGAAAAACCGCCTGCGGAACTGCATCAGATGTTCCTTGAggtggt                                              <  2:1150152/100‑1 (MQ=255)

                                                   tGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCGGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGtggtgg                                            <  1:3645604/100‑1 (MQ=255)

                                                   tGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGtggtgg                                            >  1:3030333/1‑100 (MQ=255)

                                                     ggATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGAt                                          >  2:3000411/1‑100 (MQ=255)

                                                     ggATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGAt                                          <  1:5341177/100‑1 (MQ=255)

                                                     ggATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGAt                                          >  2:442722/1‑100 (MQ=255)

                                                      gATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATc                                         >  1:3999777/1‑100 (MQ=255)

                                                           gagTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATcgcgcg                                    <  2:1749242/100‑1 (MQ=255)

                                                           gagTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATcgcgcg                                    >  1:2112626/1‑100 (MQ=255)

                                                             gTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTCGAGGTGGTGGATcgcgcgct                                  >  2:1944251/1‑98 (MQ=255)

                                                             gTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGAATGCGGAACTGCATCAGCTGTGCCTTGAGGTGGTGGATCGCGCGGt                                  >  1:4189777/1‑100 (MQ=255)

                                                             gTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAAAAAAAAAAAAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGt                                  <  2:5121911/100‑1 (MQ=255)

                                                                gggCTTACCGTTTAACTCTGCGCCAGATTGAACCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGCGGTGGATCGCGCGGTaca                               >  2:2009294/1‑98 (MQ=255)

                                                                 ggCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTaaaa                              <  2:3936386/100‑1 (MQ=255)

                                                                  gCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAg                             >  1:4529715/1‑100 (MQ=255)

                                                                  gCTTACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAg                             >  1:1680225/1‑100 (MQ=255)

                                                                    ttACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGAt                           <  1:2499814/100‑1 (MQ=255)

                                                                     tACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATg                          >  1:1092216/1‑100 (MQ=255)

                                                                     tACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATg                          >  1:3148735/1‑100 (MQ=255)

                                                                     tACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATg                          >  1:3389591/1‑100 (MQ=255)

                                                                     tACCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATg                          >  1:5042109/1‑100 (MQ=255)

                                                                      aCCGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGa                         <  1:4201677/100‑1 (MQ=255)

                                                                       ccGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa                        >  1:3157737/1‑100 (MQ=255)

                                                                       ccGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa                        >  1:4780277/1‑100 (MQ=255)

                                                                       ccGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGaa                        >  2:4203856/1‑100 (MQ=255)

                                                                        cGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAag                       >  2:1234685/1‑100 (MQ=255)

                                                                        cGTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAag                       >  2:4504296/1‑100 (MQ=255)

                                                                         gTTTTACTCTGCGCCAGGTTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGTTGTTGGACCGCGCGGTACAAGATGAaga                      >  2:1260381/1‑100 (MQ=255)

                                                                         gTTTTACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAaga                      >  1:1872280/1‑100 (MQ=255)

                                                                            ttACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGAt                   <  2:4191306/100‑1 (MQ=255)

                                                                             tACTCTGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATc                  <  1:2023048/100‑1 (MQ=255)

                                                                               ctctGGGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCt                >  2:4839748/1‑100 (MQ=255)

                                                                               ctctGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCt                <  1:5068977/100‑1 (MQ=255)

                                                                               ctctGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCt                >  2:4940819/1‑100 (MQ=255)

                                                                               ctctGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAAAAATCCt                >  2:224103/1‑100 (MQ=255)

                                                                                 ctGCGCCTGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGa              >  2:4081456/1‑100 (MQ=255)

                                                                                 ctGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTGGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGa              >  2:828557/1‑100 (MQ=255)

                                                                                  tGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGAc             <  2:3510598/100‑1 (MQ=255)

                                                                                  tGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGAc             <  2:2421940/100‑1 (MQ=255)

                                                                                  tGCGCCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCCCGGTAAAAGATGAAGAGATCCTGAc             <  1:2352452/100‑1 (MQ=255)

                                                                                   gcgcCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACg            <  2:2795683/100‑1 (MQ=255)

                                                                                   gcgcCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACg            <  1:1427709/100‑1 (MQ=255)

                                                                                   gcgcCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACg            <  1:97924/100‑1 (MQ=255)

                                                                                   gcgcCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACg            <  2:2795862/100‑1 (MQ=255)

                                                                                    cgcCAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGc           >  1:180932/1‑100 (MQ=255)

                                                                                      ccAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCaa         >  1:2136063/1‑100 (MQ=255)

                                                                                       cAGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCAAc        <  1:1032568/100‑1 (MQ=255)

                                                                                        aGATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCAACt       <  2:2086893/100‑1 (MQ=255)

                                                                                         gATTGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCAACTg      <  2:47180/100‑1 (MQ=255)

                                                                                           ttGAAGCACAGAAAGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCAACTGGc    <  2:367402/100‑1 (MQ=255)

                                                                                            tGAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCAACTGGCg   <  1:3130215/100‑1 (MQ=255)

                                                                                             gAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCAACTGGCGa  >  1:4691530/1‑100 (MQ=255)

                                                                                             gAAGCGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCAACTGGCGa  >  2:5467831/1‑100 (MQ=255)

ATCGCCGCTGACAACGGTTTCGACTTTCATATCATCGACAATGAAATCTATTGGGATGAGAGTCGGGCTTACCGTTTTACTCTGCGCCAGATTGAAGAGCAGATCGAAAAACCGACTGCGGAACTGCATCAGATGTGCCTTGAGGTGGTGGATCGCGCGGTAAAAGATGAAGAGATCCTGACGCAACTGGCGA  >  NC_000913/4413070‑4413262

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