BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,408,560	G→A	A144T (GCC→ACC)	yjfJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,408,560	0	G	A	99.0%	705.8 / NA	203	A144T (GCC→ACC)	yjfJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (101/100); minor base C (0/2); total (101/102)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

CTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTA  >  NC_000913/4408470‑4408649

  ggAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       >  2:286027/1‑93 (MQ=255)

   gAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACt                                                                                      <  1:290113/93‑1 (MQ=255)

    aaGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTg                                                                                     <  2:2491412/93‑1 (MQ=255)

     aGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGa                                                                                    <  2:2611578/93‑1 (MQ=255)

     aGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGa                                                                                    >  2:3885763/1‑93 (MQ=255)

       tctcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAg                                                                                  <  2:267240/93‑1 (MQ=255)

       tctcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAg                                                                                  >  2:3251640/1‑93 (MQ=255)

       tctcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAg                                                                                  >  1:907018/1‑93 (MQ=255)

        ctcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGc                                                                                 <  1:326640/93‑1 (MQ=255)

           ccGTGATGGGGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCAt                                                                              <  1:1868616/93‑1 (MQ=255)

           ccGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCAt                                                                              >  1:2216642/1‑93 (MQ=255)

            cGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCTATg                                                                             >  1:3113812/1‑93 (MQ=255)

             gTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGc                                                                            <  1:1078853/93‑1 (MQ=255)

              tGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCa                                                                           <  1:848255/93‑1 (MQ=255)

               gATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAg                                                                          <  1:4610185/93‑1 (MQ=255)

                aTGGCGTGGAATAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGc                                                                         >  2:3798477/1‑93 (MQ=255)

                 tGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCg                                                                        >  1:2552742/1‑93 (MQ=255)

                 tGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCg                                                                        >  1:777105/1‑93 (MQ=255)

                  ggCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGt                                                                       <  1:1362233/93‑1 (MQ=255)

                     gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGca                                                                    >  1:2952765/1‑93 (MQ=255)

                       ggAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGcaca                                                                  >  1:4388457/1‑93 (MQ=255)

                         aaaaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaac                                                                >  1:3910216/1‑93 (MQ=255)

                         aaaaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaac                                                                >  1:2829924/1‑93 (MQ=255)

                         aaaaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaac                                                                >  2:3609798/1‑93 (MQ=255)

                          aaaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaaca                                                               >  2:2947423/1‑93 (MQ=255)

                          aaaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaaca                                                               >  1:3354920/1‑93 (MQ=255)

                          aaaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaaca                                                               <  2:4471622/93‑1 (MQ=255)

                            aaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAgg                                                             >  2:2173606/1‑93 (MQ=255)

                             aGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGc                                                            <  1:1585771/93‑1 (MQ=255)

                              gCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCt                                                           <  1:1424219/93‑1 (MQ=255)

                              gCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCt                                                           >  2:3170370/1‑93 (MQ=255)

                               cAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTg                                                          <  1:3527478/93‑1 (MQ=255)

                               cAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTg                                                          <  2:2730114/93‑1 (MQ=255)

                                 gTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTa                                                        >  2:238629/1‑93 (MQ=255)

                                  ttACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTac                                                       >  2:3993467/1‑92 (MQ=255)

                                  ttACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaa                                                       >  1:3027643/1‑93 (MQ=255)

                                   tACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaac                                                      <  2:2213153/93‑1 (MQ=255)

                                    acagcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaaca                                                     >  1:1004381/1‑93 (MQ=255)

                                       gcgacggcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAACCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAAcc                                                  <  2:4188941/93‑1 (MQ=255)

                                       gcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAAcc                                                  >  2:1954546/1‑93 (MQ=255)

                                         gacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctc                                                <  1:3768543/93‑1 (MQ=255)

                                          acagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctct                                               <  1:3715785/93‑1 (MQ=255)

                                          acagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctct                                               >  1:458422/1‑93 (MQ=255)

                                          acagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctct                                               >  2:637762/1‑93 (MQ=255)

                                           cagcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTa                                              <  1:581424/93‑1 (MQ=255)

                                           cagcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTa                                              <  1:2022998/93‑1 (MQ=255)

                                            agcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAc                                             <  1:2179748/93‑1 (MQ=255)

                                             gcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAcc                                            >  2:486903/1‑93 (MQ=255)

                                               gcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGt                                          <  2:2053974/93‑1 (MQ=255)

                                                cagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGtt                                         <  2:4494593/93‑1 (MQ=255)

                                                 agcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTg                                        <  2:158497/93‑1 (MQ=255)

                                                 agcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTg                                        <  2:609133/93‑1 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       >  2:3515577/1‑93 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       >  1:2433532/1‑93 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       <  2:3018709/93‑1 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       >  2:565257/1‑93 (MQ=255)

                                                   cgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGgc                                      <  1:4090897/93‑1 (MQ=255)

                                                   cgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGgc                                      >  1:1655727/1‑93 (MQ=255)

                                                    gTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGgcg                                     >  1:98440/1‑93 (MQ=255)

                                                    gTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGgcg                                     >  2:6864/1‑93 (MQ=255)

                                                      aTCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCa                                   >  2:2931514/1‑93 (MQ=255)

                                                      aTCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGCTGGCGCa                                   >  2:4349402/1‑93 (MQ=255)

                                                         gCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCt                                <  2:4270354/93‑1 (MQ=255)

                                                          cTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCtt                               >  1:4627340/1‑93 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              <  2:4648290/93‑1 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              >  1:2366032/1‑93 (MQ=255)

                                                            cAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcc                             >  2:707106/1‑93 (MQ=255)

                                                            cAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcc                             >  1:4121461/1‑93 (MQ=255)

                                                             aGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCa                            <  2:1248098/93‑1 (MQ=255)

                                                             aGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCa                            <  2:112983/93‑1 (MQ=255)

                                                              gTTTGAGCAGCAAATGGAAGTCGTTAAACCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTACCAACCTCTACCGTTGGCGCATCTTCCAg                           <  2:2142460/93‑1 (MQ=255)

                                                              gTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAg                           >  2:1773953/1‑93 (MQ=255)

                                                               tttGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGc                          <  2:4276396/93‑1 (MQ=255)

                                                               tttGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGc                          <  2:278228/93‑1 (MQ=255)

                                                               tttGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGc                          >  1:2265112/1‑93 (MQ=255)

                                                                  gAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGtt                       <  1:2132680/93‑1 (MQ=255)

                                                                   agcagcAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGttt                      <  2:1832221/93‑1 (MQ=255)

                                                                   agcagcAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGttt                      <  2:4385357/93‑1 (MQ=255)

                                                                    gcagcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTc                     <  2:1510668/93‑1 (MQ=255)

                                                                    gcagcaAATGGAAGTCGTCAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCAAGCGTTTc                     <  2:4573242/93‑1 (MQ=255)

                                                                     cagcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCCTCTTCCAGCGTTTCg                    >  2:1294504/1‑93 (MQ=255)

                                                                       gcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGAc                  <  2:1173028/93‑1 (MQ=255)

                                                                        caAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACa                 >  1:145407/1‑93 (MQ=255)

                                                                        caAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACa                 <  1:1125774/93‑1 (MQ=255)

                                                                        caAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACa                 >  2:4048336/1‑93 (MQ=255)

                                                                         aaaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAg                >  1:533295/1‑93 (MQ=255)

                                                                         aaaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAg                >  2:1326466/1‑93 (MQ=255)

                                                                         aaaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAg                <  2:2830139/93‑1 (MQ=255)

                                                                          aaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGc               <  1:1959371/93‑1 (MQ=255)

                                                                           aTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCg              <  1:2088121/93‑1 (MQ=255)

                                                                            tGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCgg             >  1:25707/1‑93 (MQ=255)

                                                                             ggAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGc            >  2:4155316/1‑93 (MQ=255)

                                                                             ggAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGc            >  1:688584/1‑93 (MQ=255)

                                                                             ggAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGc            >  2:139485/1‑93 (MQ=255)

                                                                               aaGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAg          <  1:2733360/93‑1 (MQ=255)

                                                                                aGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGa         >  1:2721377/1‑93 (MQ=255)

                                                                                 gTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGaa        >  1:4148018/1‑93 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    >  1:1754075/1‑93 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    <  1:1249126/93‑1 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCCTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    >  1:3560326/1‑93 (MQ=255)

CTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTA  >  NC_000913/4408470‑4408649

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