BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,288,263	T→C	E31E (GAA→GAG)	ychJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,288,263	0	T	C	77.7%	530.2 / 176.8	292	E31E (GAA→GAG)	ychJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (37/28); new base C (122/105); total (159/133)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-01

Only 100 of 297 total aligned reads displayed.

CGCCGCTGCTCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCTGAGACAC  >  NC_000913/1288173‑1288355

 gccgcTGCTCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGctct                                                                                           <  1:1565275/93‑1 (MQ=255)

  ccgcTGCTCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTg                                                                                          >  1:4020948/1‑93 (MQ=255)

    gctgctCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGa                                                                                        <  1:2927239/93‑1 (MQ=255)

     ctgctCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGAt                                                                                       <  2:714161/93‑1 (MQ=255)

      tgctCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATc                                                                                      >  1:5586944/1‑93 (MQ=255)

       gctCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCa                                                                                     <  2:7482547/93‑1 (MQ=255)

       gctCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCa                                                                                     >  2:3957783/1‑93 (MQ=255)

        ctCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAg                                                                                    >  1:6372238/1‑93 (MQ=255)

        ctCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAg                                                                                    <  2:1447511/93‑1 (MQ=255)

         tCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAgg                                                                                   >  2:6179398/1‑93 (MQ=255)

         tCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAgg                                                                                   >  2:1957864/1‑93 (MQ=255)

          ccACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGt                                                                                  <  1:3034160/93‑1 (MQ=255)

          ccACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGt                                                                                  >  2:5736260/1‑93 (MQ=255)

           cacaAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTg                                                                                 >  1:2592463/1‑93 (MQ=255)

            acaAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGc                                                                                <  2:5303484/93‑1 (MQ=255)

              aaGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCaa                                                                              <  1:666705/93‑1 (MQ=255)

               agagGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAAc                                                                             <  1:2155923/93‑1 (MQ=255)

               agagGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAAc                                                                             >  1:1678901/1‑93 (MQ=255)

               agagGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAAc                                                                             >  2:6782707/1‑93 (MQ=255)

                gagGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAAcc                                                                            >  1:7284577/1‑93 (MQ=255)

                 agGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCt                                                                           <  2:740282/93‑1 (MQ=255)

                 agGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCt                                                                           <  1:4510969/93‑1 (MQ=255)

                     aTGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTc                                                                       <  1:1011022/93‑1 (MQ=255)

                        ccAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCAcc                                                                    >  1:3206225/1‑93 (MQ=255)

                         cAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCa                                                                   >  2:1643340/1‑93 (MQ=255)

                         cAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCa                                                                   >  2:2703049/1‑93 (MQ=255)

                         cAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCa                                                                   >  1:2500229/1‑93 (MQ=255)

                         cAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCa                                                                   <  1:7205622/93‑1 (MQ=255)

                          aGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAg                                                                  >  1:5883729/1‑93 (MQ=255)

                          aGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAg                                                                  >  1:1078515/1‑93 (MQ=255)

                           ggTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGa                                                                 >  1:4092763/1‑93 (MQ=255)

                            gTCTTTGTTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGac                                                                >  2:4770674/1‑93 (MQ=255)

                             tCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGaca                                                               <  1:5242998/93‑1 (MQ=255)

                              cTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGacac                                                              >  2:5891127/1‑93 (MQ=255)

                               tttATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGAcaca                                                             <  2:4224989/93‑1 (MQ=255)

                               tttATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGAcaca                                                             >  2:5994657/1‑93 (MQ=255)

                               tttATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGAcaca                                                             <  2:793879/93‑1 (MQ=255)

                                   ttaAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAg                                                         >  1:2721167/1‑93 (MQ=255)

                                    taAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAgg                                                        >  2:346614/1‑93 (MQ=255)

                                     aaaTAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAggg                                                       >  2:6287117/1‑93 (MQ=255)

                                       ataatCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTg                                                     <  2:4123380/93‑1 (MQ=255)

                                        taatCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTgg                                                    >  1:1804497/1‑93 (MQ=255)

                                         aatCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGc                                                   <  2:5272680/93‑1 (MQ=255)

                                         aatCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGc                                                   >  1:4687504/1‑93 (MQ=255)

                                           tCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCaa                                                 >  2:1139302/1‑93 (MQ=255)

                                            cTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAAc                                                <  2:2107429/93‑1 (MQ=255)

                                             tGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACa                                               >  2:849150/1‑93 (MQ=255)

                                             tGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACa                                               >  1:5375373/1‑93 (MQ=255)

                                              gCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACAt                                              >  1:6265795/1‑93 (MQ=255)

                                              gCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACAt                                              <  2:1796212/93‑1 (MQ=255)

                                               cGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATa                                             <  2:6475852/93‑1 (MQ=255)

                                                gTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAg                                            <  2:6985569/93‑1 (MQ=255)

                                                 tCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAgg                                           >  1:2305828/1‑93 (MQ=255)

                                                 tCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAgg                                           >  1:1909557/1‑93 (MQ=255)

                                                 tCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAgg                                           <  1:3822262/93‑1 (MQ=255)

                                                    tGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCta                                        >  1:3563594/1‑93 (MQ=255)

                                                    tGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCta                                        >  2:6109631/1‑93 (MQ=255)

                                                    tGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCta                                        <  2:2392583/93‑1 (MQ=255)

                                                     gCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCtat                                       >  1:4177559/1‑93 (MQ=255)

                                                     gCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGGGCCAACAAAGGCaat                                       >  1:4956622/1‑93 (MQ=255)

                                                     gCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACATTTTCACCAGACACATAAGGGTGGCAACATAGGCtat                                       >  1:6427016/1‑93 (MQ=255)

                                                       aTCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATAc                                     <  2:6038691/93‑1 (MQ=255)

                                                        tCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACt                                    >  1:6507386/1‑93 (MQ=255)

                                                        tCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACt                                    >  2:4866650/1‑93 (MQ=255)

                                                        tCACAAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACt                                    <  1:1353235/93‑1 (MQ=255)

                                                        tCACAAAAGCGCAGTAACGCGAACGCATGAGAGACTCTGGATCAGGGGCAACATTATCACCAGACACATAAGGGTGGCAACATCGGCTATACt                                    >  2:2332768/1‑93 (MQ=255)

                                                         cacaAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTc                                   <  1:4257706/93‑1 (MQ=255)

                                                          acaAAAGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCg                                  >  2:1738807/1‑93 (MQ=255)

                                                             aaaGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACa                               >  1:3684331/1‑93 (MQ=255)

                                                             aaaGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACa                               >  1:1209571/1‑93 (MQ=255)

                                                              aaGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAg                              >  1:3321587/1‑93 (MQ=255)

                                                               aGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGc                             >  1:697396/1‑93 (MQ=255)

                                                               aGCGCAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGc                             >  2:2041761/1‑93 (MQ=255)

                                                                gcgcAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCa                            >  1:2083940/1‑93 (MQ=255)

                                                                 cgcAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCAc                           <  2:2964945/93‑1 (MQ=255)

                                                                  gcAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACt                          >  2:4555872/1‑93 (MQ=255)

                                                                  gcAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACt                          >  1:5456843/1‑93 (MQ=255)

                                                                   cAGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTa                         >  2:7109628/1‑93 (MQ=255)

                                                                    aGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTAc                        <  2:3663401/93‑1 (MQ=255)

                                                                    aGTAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTAc                        >  2:657353/1‑93 (MQ=255)

                                                                      tAACGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACca                      >  2:4413817/1‑93 (MQ=255)

                                                                       aaCGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACcac                     <  2:2921776/93‑1 (MQ=255)

                                                                        aCGCGAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACcaca                    <  2:1371709/93‑1 (MQ=255)

                                                                            gAACGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGa                <  2:368138/93‑1 (MQ=255)

                                                                             aaCGCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGAc               >  1:771243/1‑93 (MQ=255)

                                                                                gCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACaaa            >  2:267306/1‑93 (MQ=255)

                                                                                gCATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACaaa            <  2:1763528/93‑1 (MQ=255)

                                                                                 cATGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAg           >  1:6557321/1‑93 (MQ=255)

                                                                                  aTGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGc          >  1:7436207/1‑93 (MQ=255)

                                                                                  aTGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGc          >  1:6859544/1‑93 (MQ=255)

                                                                                  aTGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGc          <  1:5974946/93‑1 (MQ=255)

                                                                                   tGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCt         <  2:3186660/93‑1 (MQ=255)

                                                                                   tGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCt         <  1:6605332/93‑1 (MQ=255)

                                                                                   tGAGATGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCt         >  2:3552490/1‑93 (MQ=255)

                                                                                    gagaTGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCTg        >  2:5090505/1‑93 (MQ=255)

                                                                                     agaTGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCTga       <  1:950163/93‑1 (MQ=255)

                                                                                      gaTGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCTgag      <  2:5441494/93‑1 (MQ=255)

                                                                                       aTGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCTgaga     <  1:5749227/93‑1 (MQ=255)

                                                                                       aTGCTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCTgaga     <  2:3985593/93‑1 (MQ=255)

                                                                                          ttCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCTGAGacac  >  2:5597356/1‑93 (MQ=255)

CGCCGCTGCTCCACAAGAGGGATGCCAGGTCTTTATTAAATAATCTGCGTCTTGCATCACAAAAGCGCAGTAACGCGAACGCATGAGATGTTCTGGATCAGGTGCAACCTTTTCACCAGACACATAAGGGTGGCAACATAGGCTATACTCGACAGCACTACCACAGGGACAAAGCTGAGACAC  >  NC_000913/1288173‑1288355

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