BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,782,476	A→G	V56A (GTT→GCT)	yibN ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,782,476	0	A	G	89.1%	923.2 / 84.4	350	V56A (GTT→GCT)	yibN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (19/19); new base G (167/145); total (186/164)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

Only 100 of 355 total aligned reads displayed.

ACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGA  >  NC_000913/3782386‑3782568

aCATTGTTGGCGTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAAcc                                                                                            >  2:916008/1‑93 (MQ=255)

  aTTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCcac                                                                                          <  1:1705173/93‑1 (MQ=255)

       tgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGt                                                                                     <  1:5023243/93‑1 (MQ=255)

       tgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGt                                                                                     >  2:1762212/1‑93 (MQ=255)

        ggCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtc                                                                                    <  1:6247919/93‑1 (MQ=255)

         gCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtct                                                                                   <  2:267113/93‑1 (MQ=255)

          cTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtctt                                                                                  <  1:1066205/93‑1 (MQ=255)

          cTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtctt                                                                                  >  1:7048909/1‑93 (MQ=255)

          cTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtctt                                                                                  <  2:4551496/93‑1 (MQ=255)

           tttGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtcttc                                                                                 >  1:1462188/1‑93 (MQ=255)

            ttGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGtcttct                                                                                <  2:4209538/93‑1 (MQ=255)

               aTTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTttg                                                                             >  1:2811918/1‑93 (MQ=255)

               aTTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTttg                                                                             >  2:3452107/1‑93 (MQ=255)

                tttCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTttgt                                                                            <  2:7058502/93‑1 (MQ=255)

                tttCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTttgt                                                                            <  1:5963106/93‑1 (MQ=255)

                tttCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTttgt                                                                            >  1:5088710/1‑93 (MQ=255)

                      tCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCa                                                                      <  2:7401101/93‑1 (MQ=255)

                      tCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGCCTTCTTTGTTGATCa                                                                      >  1:1138628/1‑93 (MQ=255)

                       cGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAg                                                                     <  2:162771/93‑1 (MQ=255)

                       cGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAg                                                                     >  2:2791166/1‑93 (MQ=255)

                         gCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGAc                                                                   >  1:6130266/1‑93 (MQ=255)

                          cAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAATTCATCACGCTTACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGAcg                                                                  <  1:4771724/93‑1 (MQ=255)

                           aaCAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGAcgc                                                                 <  1:7259350/93‑1 (MQ=255)

                           aaCAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGAcgc                                                                 >  1:5148919/1‑93 (MQ=255)

                            aCAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGAcgcg                                                                >  2:1862313/1‑93 (MQ=255)

                            aCAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGAcgcg                                                                >  1:6557064/1‑93 (MQ=255)

                            aCAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGAcgcg                                                                >  1:6972898/1‑93 (MQ=255)

                             cAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGt                                                               <  1:4604311/93‑1 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                              <  1:4314023/93‑1 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                              <  1:3327512/93‑1 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                              >  2:1210178/1‑93 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                              <  2:6824291/93‑1 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTATGTTGATCAGACGCGTa                                                              <  1:2992139/93‑1 (MQ=255)

                              aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                              <  2:5550288/93‑1 (MQ=255)

                                aTTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGc                                                            >  1:4081835/1‑93 (MQ=255)

                                 ttAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCt                                                           >  2:709528/1‑93 (MQ=255)

                                    aTAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTca                                                        >  2:2209395/1‑93 (MQ=255)

                                    aTAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTca                                                        <  1:1388714/93‑1 (MQ=255)

                                    aTAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTca                                                        <  2:2109441/93‑1 (MQ=255)

                                      aGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTcacc                                                      >  2:5361131/1‑93 (MQ=255)

                                       gAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTcacca                                                     <  2:3318364/93‑1 (MQ=255)

                                        aaCCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTcaccac                                                    <  2:4133274/93‑1 (MQ=255)

                                          ccTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGa                                                  <  2:5719588/93‑1 (MQ=255)

                                          ccTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGa                                                  <  2:3583690/93‑1 (MQ=255)

                                           cTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAg                                                 >  1:1780095/1‑93 (MQ=255)

                                            tGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGt                                                >  2:4331413/1‑93 (MQ=255)

                                             gCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTa                                               <  2:4833823/93‑1 (MQ=255)

                                             gCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTa                                               <  1:2423473/93‑1 (MQ=255)

                                              cGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTaa                                              <  2:2926166/93‑1 (MQ=255)

                                               gATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCATCACCACGAGTAAt                                             >  2:198309/1‑93 (MQ=255)

                                               gATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAAt                                             >  2:5222880/1‑93 (MQ=255)

                                               gATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAAt                                             >  2:2586871/1‑93 (MQ=255)

                                                atatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATc                                            <  1:3169708/93‑1 (MQ=255)

                                                atatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATc                                            <  2:5975529/93‑1 (MQ=255)

                                                atatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATc                                            <  2:6760709/93‑1 (MQ=255)

                                                atatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATc                                            >  2:1303001/1‑93 (MQ=255)

                                                atatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATc                                            <  1:6711443/93‑1 (MQ=255)

                                                 tatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCa                                           >  1:3250723/1‑93 (MQ=255)

                                                  atGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCAc                                          >  2:3855915/1‑93 (MQ=255)

                                                  atGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCAc                                          >  2:1484997/1‑93 (MQ=255)

                                                   tGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCAcc                                         >  1:7505460/1‑93 (MQ=255)

                                                    ggCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCt                                        <  1:4020688/93‑1 (MQ=255)

                                                    ggCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCt                                        >  1:967328/1‑93 (MQ=255)

                                                     gCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCtt                                       >  1:7192413/1‑93 (MQ=255)

                                                      ccTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTc                                      >  1:5589095/1‑93 (MQ=255)

                                                       cTTTACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCa                                     >  2:6496958/1‑93 (MQ=255)

                                                        tttACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCAc                                    <  2:4737823/93‑1 (MQ=255)

                                                        tttACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCAc                                    >  2:5046481/1‑93 (MQ=255)

                                                        tttACGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCAc                                    >  1:2661196/1‑93 (MQ=255)

                                                           aCGGAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACttt                                 >  2:6633076/1‑93 (MQ=255)

                                                           aCGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACttt                                 >  1:5793430/1‑93 (MQ=255)

                                                            cGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTc                                >  2:140940/1‑93 (MQ=255)

                                                              gAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGa                              >  2:545676/1‑93 (MQ=255)

                                                              gAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGa                              >  2:7206373/1‑93 (MQ=255)

                                                              gAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGa                              >  1:4829430/1‑93 (MQ=255)

                                                              gAAGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGa                              >  1:4715193/1‑93 (MQ=255)

                                                              gAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGa                              <  2:5293060/93‑1 (MQ=255)

                                                               aaGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAg                             >  2:4723933/1‑93 (MQ=255)

                                                                aGTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAgg                            <  2:2438585/93‑1 (MQ=255)

                                                                 gTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGt                           >  1:6197109/1‑93 (MQ=255)

                                                                 gTCATCACGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGt                           <  1:85505/93‑1 (MQ=255)

                                                                  tcatcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTc                          <  1:4387599/93‑1 (MQ=255)

                                                                    atcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAg                        <  1:1971363/93‑1 (MQ=255)

                                                                    atcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAg                        >  2:7509580/1‑93 (MQ=255)

                                                                    atcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAg                        >  1:4992016/1‑93 (MQ=255)

                                                                    atcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAg                        >  1:425278/1‑93 (MQ=255)

                                                                    atcaCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAg                        >  1:2348734/1‑93 (MQ=255)

                                                                      caCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGc                      <  1:3739323/93‑1 (MQ=255)

                                                                      caCGCTGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGc                      <  2:6321581/93‑1 (MQ=255)

                                                                           tGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTa                 >  1:3492244/1‑93 (MQ=255)

                                                                           tGACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTa                 >  2:442840/1‑93 (MQ=255)

                                                                            gACGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTaa                >  1:7480172/1‑93 (MQ=255)

                                                                             aCGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTaaa               <  2:3175458/93‑1 (MQ=255)

                                                                             aCGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTaaa               >  1:2031761/1‑93 (MQ=255)

                                                                              cGTAAATCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAAc              <  1:1870224/93‑1 (MQ=255)

                                                                                tAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGt            <  1:2209994/93‑1 (MQ=255)

                                                                                 aaaTCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTa           >  2:2376206/1‑93 (MQ=255)

                                                                                  aaTCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAg          >  2:6166004/1‑93 (MQ=255)

                                                                                  aaTCCACAGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAg          >  2:6736631/1‑93 (MQ=255)

                                                                                       acaGCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTcaca     >  1:5491066/1‑93 (MQ=255)

ACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGA  >  NC_000913/3782386‑3782568

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