BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,137,219	C→T	G24R (GGG→AGG)	talC ←	putative transaldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,219	0	C	T	97.5%	559.0 / ‑0.3	161	G24R (GGG→AGG)	talC	putative transaldolase
Reads supporting (aligned to +/- strand): ref base C (2/2); new base T (81/76); total (83/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

Only 100 of 164 total aligned reads displayed.

CCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTT  >  NC_000913/4137129‑4137307

ccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTgg                                                                                        >  2:1233532/1‑93 (MQ=255)

 cTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGc                                                                                       >  2:529808/1‑93 (MQ=255)

 cTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGc                                                                                       >  2:1641545/1‑93 (MQ=255)

 cTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGc                                                                                       >  1:3783907/1‑93 (MQ=255)

   catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCaa                                                                                     >  1:2729687/1‑93 (MQ=255)

   catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCaa                                                                                     >  2:1474235/1‑93 (MQ=255)

   catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCaa                                                                                     <  2:436630/93‑1 (MQ=255)

    atcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAAt                                                                                    >  2:1119793/1‑93 (MQ=255)

    atcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAAt                                                                                    <  1:3123978/93‑1 (MQ=255)

    atcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAAt                                                                                    <  2:1233288/93‑1 (MQ=255)

    atcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAAt                                                                                    >  1:3819643/1‑93 (MQ=255)

       aCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATggg                                                                                 <  1:2169731/93‑1 (MQ=255)

        ccAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATgggg                                                                                <  1:4361656/93‑1 (MQ=255)

         cAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGa                                                                               >  1:4138014/1‑93 (MQ=255)

         cAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGa                                                                               >  2:2422914/1‑93 (MQ=255)

          aaTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGaa                                                                              >  2:3424084/1‑93 (MQ=255)

          aaTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGaa                                                                              >  1:3362178/1‑93 (MQ=255)

          aaTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGaa                                                                              <  2:708393/93‑1 (MQ=255)

           aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAat                                                                             <  2:209311/93‑1 (MQ=255)

            tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                            <  1:11285/93‑1 (MQ=255)

            tCGCTGTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                            <  2:2176266/93‑1 (MQ=255)

             cGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                           <  2:899915/93‑1 (MQ=255)

             cGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                           >  1:2726336/1‑93 (MQ=255)

             cGCTTTTTGCAGAAGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                           <  1:1583106/93‑1 (MQ=255)

              gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                          >  2:2536178/1‑93 (MQ=255)

                tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                                        >  2:1280992/1‑93 (MQ=255)

                tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                                        >  2:3811493/1‑93 (MQ=255)

                 ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGt                                                                       <  2:1444689/93‑1 (MQ=255)

                     gCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGccg                                                                   >  1:1408929/1‑92 (MQ=255)

                       aGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGa                                                                 >  1:789989/1‑93 (MQ=255)

                       aGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGa                                                                 <  1:916039/93‑1 (MQ=255)

                       aGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGa                                                                 <  1:2475258/93‑1 (MQ=255)

                        gACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGAc                                                                <  1:1763628/93‑1 (MQ=255)

                        gACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGAc                                                                <  1:1520681/93‑1 (MQ=255)

                        gACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGAc                                                                <  1:697536/93‑1 (MQ=255)

                         aCGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACg                                                               <  2:1130995/93‑1 (MQ=255)

                          cgcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGt                                                              >  2:4630092/1‑93 (MQ=255)

                          cgcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGt                                                              <  1:228448/93‑1 (MQ=255)

                           gcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGtt                                                             <  1:1987028/93‑1 (MQ=255)

                            cgtcagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTc                                                            <  2:1265575/90‑1 (MQ=255)

                            cgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTc                                                            <  2:24154/93‑1 (MQ=255)

                             ggCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCg                                                           <  2:212258/93‑1 (MQ=255)

                               cagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGAc                                                         <  2:2104825/93‑1 (MQ=255)

                                agcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACt                                                        >  1:1894916/1‑93 (MQ=255)

                                 gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACtt                                                       <  1:1022526/93‑1 (MQ=255)

                                  caCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTc                                                      <  2:4095415/93‑1 (MQ=255)

                                   aCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCt                                                     <  2:18283/93‑1 (MQ=255)

                                     ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                                   <  2:919727/93‑1 (MQ=255)

                                     ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                                   >  1:575280/1‑93 (MQ=255)

                                        ccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGAc                                                <  2:1309245/93‑1 (MQ=255)

                                         cATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACg                                               >  1:1040031/1‑93 (MQ=255)

                                          atatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTTTGCGACGt                                              <  2:3839589/93‑1 (MQ=255)

                                           tatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                             >  2:2712614/1‑93 (MQ=255)

                                           tatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                             <  2:4661547/93‑1 (MQ=255)

                                           tatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                             >  1:2016213/1‑93 (MQ=255)

                                            atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGGCGTTa                                            >  2:1709696/1‑93 (MQ=255)

                                            atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTa                                            >  2:1979422/1‑93 (MQ=255)

                                            atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTa                                            <  2:2932738/93‑1 (MQ=255)

                                            atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTa                                            <  1:2200255/93‑1 (MQ=255)

                                             tGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAg                                           <  2:3365971/93‑1 (MQ=255)

                                             tGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAg                                           >  1:3324840/1‑93 (MQ=255)

                                              ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                          <  1:1984158/93‑1 (MQ=255)

                                              ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                          <  2:2743588/93‑1 (MQ=255)

                                                  tCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgtg                                      <  2:3115967/93‑1 (MQ=255)

                                                  tCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgtg                                      >  1:4607022/1‑93 (MQ=255)

                                                     ttGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                                   <  2:4542038/93‑1 (MQ=255)

                                                     ttGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                                   >  2:321783/1‑93 (MQ=255)

                                                     ttGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTcc                                   >  2:3132904/1‑93 (MQ=255)

                                                      tGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCa                                  <  2:3484929/93‑1 (MQ=255)

                                                        cTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGa                                >  1:3868998/1‑93 (MQ=255)

                                                         tGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGAt                               >  1:59875/1‑93 (MQ=255)

                                                           gCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATAc                             >  2:4406430/1‑93 (MQ=255)

                                                            cAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACa                            <  1:3539446/93‑1 (MQ=255)

                                                            cAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACa                            <  1:2553566/93‑1 (MQ=255)

                                                             aGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAg                           <  2:3710117/93‑1 (MQ=255)

                                                              gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                          <  2:3946946/93‑1 (MQ=255)

                                                              gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                          >  1:1800684/1‑93 (MQ=255)

                                                              gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                          >  1:4675514/1‑93 (MQ=255)

                                                                gATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTc                        <  1:502806/93‑1 (MQ=255)

                                                                 aTAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTcc                       <  1:1845668/93‑1 (MQ=255)

                                                                  tAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGACTATGCGTGGCAGACGGTCGACTTCTGCGACGTGAGCGGTGTCCAGGTACAGTTGCa                      >  2:2178878/1‑93 (MQ=255)

                                                                  tAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGATCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTTCa                      >  2:2812083/1‑93 (MQ=255)

                                                                   aaTGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCat                     <  1:31751/93‑1 (MQ=255)

                                                                   aaTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCat                     >  2:2137492/1‑93 (MQ=255)

                                                                    aTGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatg                    >  1:1002453/1‑93 (MQ=255)

                                                                    aTGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatg                    <  1:2583035/93‑1 (MQ=255)

                                                                        tCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGt                <  2:1703220/93‑1 (MQ=255)

                                                                        tCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGt                >  2:1187752/1‑93 (MQ=255)

                                                                         cGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGtt               >  2:855424/1‑93 (MQ=255)

                                                                            gttggttgTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTcc            <  2:2231449/93‑1 (MQ=255)

                                                                             ttAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCt           <  2:4021495/93‑1 (MQ=255)

                                                                              tggTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtt          >  1:1926890/3‑93 (MQ=255)

                                                                                gTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtttt        >  1:862267/1‑93 (MQ=255)

                                                                                 ttGTCACCCTGGCAATGGGGAATATGCGTGCCAGCCGTTCGACTTATGCTACGTTCGCGGGGTCCAGATACAGTTCCATGATGTTTCCTTTTa       >  1:4528266/1‑93 (MQ=255)

                                                                                   gTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGCGTCCAGATACAGTTCCATGATGTTTCCTTTTAca     >  1:1356995/1‑93 (MQ=255)

                                                                                    tCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcag    >  1:2824913/1‑93 (MQ=255)

                                                                                    tCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcag    >  2:3008112/1‑93 (MQ=255)

                                                                                     cACCCTGGCAATGGGGAATAGGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGCTACAGTTCCATGATGTTTCCTTTTAcagt   >  2:3548125/1‑93 (MQ=255)

                                                                                      cccctGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagtt  <  1:392786/92‑1 (MQ=255)

                                                                                      aCCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTATAcagtt  >  1:682180/1‑93 (MQ=255)

CCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTT  >  NC_000913/4137129‑4137307

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