BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,097,268	T→C	H280R (CAC→CGC)	rhaT ←	RhaT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,097,268	0	T	C	99.6%	802.5 / ‑4.5	232	H280R (CAC→CGC)	rhaT	RhaT
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (112/119); total (112/120)
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.03e-01

Only 100 of 236 total aligned reads displayed.

AGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCC  >  NC_000913/4097178‑4097357

aGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCgg                                                                                         >  1:37782/1‑93 (MQ=255)

aGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCgg                                                                                         <  2:611234/93‑1 (MQ=255)

 gCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGc                                                                                        >  2:3728611/1‑93 (MQ=255)

   cccAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGccc                                                                                      <  1:425468/92‑1 (MQ=255)

    ccAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGcccc                                                                                     <  2:5173790/93‑1 (MQ=255)

     cAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCa                                                                                    >  2:4181381/1‑93 (MQ=255)

      aGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAg                                                                                   <  1:2338256/93‑1 (MQ=255)

      aGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAg                                                                                   <  1:361607/93‑1 (MQ=255)

       gCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAgg                                                                                  >  1:1992283/1‑93 (MQ=255)

        cccGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGc                                                                                 <  1:303759/93‑1 (MQ=255)

          cgacgaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCAt                                                                               <  2:5530628/93‑1 (MQ=255)

          cgacgaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCAt                                                                               >  2:2314815/1‑93 (MQ=255)

           gacgaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATa                                                                              >  1:655591/1‑93 (MQ=255)

           gacgaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATa                                                                              <  1:3566179/93‑1 (MQ=255)

            acgaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAg                                                                             >  2:1926606/1‑93 (MQ=255)

            acgaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAg                                                                             >  1:5165433/1‑93 (MQ=255)

            acgaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAg                                                                             <  2:720835/93‑1 (MQ=255)

              gaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGaa                                                                           >  2:1846199/1‑93 (MQ=255)

              gaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGaa                                                                           <  2:5988734/93‑1 (MQ=255)

               aTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGaaa                                                                          >  1:3570166/1‑93 (MQ=255)

                tACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGaaaa                                                                         <  1:1011627/93‑1 (MQ=255)

                tACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGaaaa                                                                         <  2:5378762/93‑1 (MQ=255)

                 aCCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGaaaaa                                                                        >  1:4031936/1‑93 (MQ=255)

                  ccgccgCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAg                                                                       >  2:4696196/1‑93 (MQ=255)

                  ccgccgCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAg                                                                       >  2:5511244/1‑93 (MQ=255)

                   cgccgcACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGa                                                                      >  1:1450334/1‑93 (MQ=255)

                   cgccgcACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGa                                                                      >  2:6359953/1‑93 (MQ=255)

                      cgcACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAAtt                                                                   <  1:3345153/93‑1 (MQ=255)

                      cgcACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAAtt                                                                   <  2:2070539/93‑1 (MQ=255)

                      cgcACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAAtt                                                                   >  2:1488232/1‑93 (MQ=255)

                        cacaATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGc                                                                 >  1:598412/1‑93 (MQ=255)

                         acaATACATAGAAACTCATATGCAGCATCTAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCa                                                                <  2:4852324/93‑1 (MQ=255)

                         acaATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCGCCAGGCATAGAAAAAGAATTGCa                                                                >  1:6378763/1‑93 (MQ=255)

                         acaATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCa                                                                <  2:6334906/93‑1 (MQ=255)

                         acaATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCa                                                                >  1:1559504/1‑93 (MQ=255)

                         acaATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCa                                                                >  1:6015941/1‑93 (MQ=255)

                          caATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAg                                                               >  2:6276619/1‑93 (MQ=255)

                           aaTACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGa                                                              <  2:1534630/93‑1 (MQ=255)

                            aTACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGAt                                                             <  2:3432613/93‑1 (MQ=255)

                            aTACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGAt                                                             >  2:6152923/1‑93 (MQ=255)

                              aCATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATAc                                                           <  1:5041380/93‑1 (MQ=255)

                               cATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATAcc                                                          >  2:3855809/1‑93 (MQ=255)

                                 tAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACcac                                                        >  1:97558/1‑93 (MQ=255)

                                  aGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACcaca                                                       >  2:4264609/1‑93 (MQ=255)

                                   gAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACAt                                                      >  1:5835335/1‑93 (MQ=255)

                                   gAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACAt                                                      >  2:5716480/1‑93 (MQ=255)

                                      aCTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCaa                                                   <  1:2873930/93‑1 (MQ=255)

                                       cTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAAc                                                  <  2:590703/93‑1 (MQ=255)

                                       cTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAAc                                                  <  2:216160/93‑1 (MQ=255)

                                        tCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAAcc                                                 <  2:3350230/93‑1 (MQ=255)

                                        tCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAAcc                                                 <  2:941228/93‑1 (MQ=255)

                                         cATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAAccc                                                <  1:5231720/93‑1 (MQ=255)

                                          atatGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCg                                               >  2:2207824/1‑93 (MQ=255)

                                           tatGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGc                                              >  2:3455539/1‑93 (MQ=255)

                                             tGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGCAAAAGAATTGCAGATACCACATCAACCCGccc                                            >  1:4473119/1‑93 (MQ=255)

                                             tGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGccc                                            >  1:996858/1‑93 (MQ=255)

                                              gcagcaTCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCa                                           >  1:1310110/1‑93 (MQ=255)

                                                agcaTCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAgt                                         >  2:1666994/1‑93 (MQ=255)

                                                 gcaTCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAgtg                                        <  1:3082012/93‑1 (MQ=255)

                                                 gcaTCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAgtg                                        <  2:3815238/93‑1 (MQ=255)

                                                  caTCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAgtgt                                       <  2:2945885/93‑1 (MQ=255)

                                                    tCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCg                                     >  1:4855048/1‑93 (MQ=255)

                                                     ccAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCga                                    <  1:4536819/93‑1 (MQ=255)

                                                     ccAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCga                                    >  2:180694/1‑93 (MQ=255)

                                                       aaCTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCgaga                                  <  1:853951/93‑1 (MQ=255)

                                                        aCTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCTACCCGCCCAGTGTCGagag                                 <  1:6297074/93‑1 (MQ=255)

                                                        aCTGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGagag                                 <  2:130886/93‑1 (MQ=255)

                                                          tGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTa                               >  2:1083753/1‑93 (MQ=255)

                                                          tGATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTa                               <  2:5948792/93‑1 (MQ=255)

                                                           gATGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTaa                              <  2:5070722/93‑1 (MQ=255)

                                                            aTGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAac                             <  1:2471540/93‑1 (MQ=255)

                                                            aTGTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTATCGAGAGTAac                             <  2:2888714/93‑1 (MQ=255)

                                                              gTAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAacac                           >  2:4794653/1‑93 (MQ=255)

                                                               tAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAAcaca                          <  1:5578196/93‑1 (MQ=255)

                                                               tAGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAAcaca                          >  2:6050415/1‑93 (MQ=255)

                                                                aGTCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACAt                         <  1:887978/93‑1 (MQ=255)

                                                                  tCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATtg                       <  1:791400/93‑1 (MQ=255)

                                                                  tCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATtg                       <  1:5658205/93‑1 (MQ=255)

                                                                  tCATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATtg                       >  2:4098762/1‑93 (MQ=255)

                                                                   cATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATtgt                      >  2:4277185/1‑93 (MQ=255)

                                                                   cATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATtgt                      <  1:2773113/93‑1 (MQ=255)

                                                                   cATACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATtgt                      <  2:5703280/93‑1 (MQ=255)

                                                                    aTACTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATtgtg                     >  2:3219658/1‑93 (MQ=255)

                                                                       cTGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAAt                  <  2:5212823/93‑1 (MQ=255)

                                                                        tGCGCCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATg                 >  1:4135620/1‑93 (MQ=255)

                                                                         gcgcCGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGa                <  1:2234324/93‑1 (MQ=255)

                                                                            ccGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCACCCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCa             <  1:4720261/93‑1 (MQ=255)

                                                                            ccGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCa             <  2:2816568/93‑1 (MQ=255)

                                                                            ccGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCa             >  2:1183561/1‑93 (MQ=255)

                                                                             cGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCACCCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAg            <  1:3758597/93‑1 (MQ=255)

                                                                             cGGAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAg            >  1:5423619/1‑93 (MQ=255)

                                                                              ggAATGCGGGCGCGGCGCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCAGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGc           <  2:2496331/93‑1 (MQ=255)

                                                                              ggAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGc           >  1:4771137/1‑93 (MQ=255)

                                                                              ggAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGc           >  2:459291/1‑93 (MQ=255)

                                                                               gAATGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCg          <  1:1700143/93‑1 (MQ=255)

                                                                                 aTGCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGAt        <  2:3276317/93‑1 (MQ=255)

                                                                                   gCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGAttt      <  2:165056/93‑1 (MQ=255)

                                                                                   gCGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGAttt      <  2:3581621/93‑1 (MQ=255)

                                                                                    cGGGCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGAtttt     <  1:1613179/93‑1 (MQ=255)

                                                                                     gggCGCGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTg    >  2:2624463/1‑93 (MQ=255)

AGCACCAGCCCGACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCC  >  NC_000913/4097178‑4097357

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