BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,597,539	T→C	intergenic (‑250/+21)	b3460 ← / ← rpoH	LivJ/sigma(32) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,597,539	0	T	C	97.4%	720.0 / 1.9	231	intergenic (‑250/+21)	b3460/rpoH	LivJ/sigma(32) factor
Reads supporting (aligned to +/- strand): ref base T (5/1); new base C (107/118); total (112/119)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628

aCAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtc                                                                                         <  1:2283911/93‑1 (MQ=255)

  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       <  2:1560212/93‑1 (MQ=255)

  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       >  2:3003642/1‑93 (MQ=255)

  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       >  2:4389180/1‑93 (MQ=255)

   aTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGctctct                                                                                      <  1:790690/93‑1 (MQ=255)

     cTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGc                                                                                    <  2:4682478/93‑1 (MQ=255)

          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCCTCCAGGGCTCTCTGCTTAAt                                                                               >  2:650873/1‑93 (MQ=255)

          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAt                                                                               <  2:2246175/93‑1 (MQ=255)

          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAt                                                                               <  1:3760295/93‑1 (MQ=255)

           gaggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATa                                                                              <  1:926371/91‑1 (MQ=255)

               aTTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCgg                                                                          <  1:3435614/93‑1 (MQ=255)

                ttGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGa                                                                         <  2:3936586/93‑1 (MQ=255)

                  ggCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGaaa                                                                       >  1:3872593/1‑93 (MQ=255)

                  ggCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGaaa                                                                       <  1:221214/93‑1 (MQ=255)

                   gCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAt                                                                      >  1:3362497/1‑93 (MQ=255)

                   gCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAt                                                                      <  1:571354/93‑1 (MQ=255)

                    cGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAtt                                                                     >  2:4290984/1‑93 (MQ=255)

                      ttttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTAc                                                                   <  1:3526943/93‑1 (MQ=255)

                      ttttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTAc                                                                   <  1:1769346/93‑1 (MQ=255)

                       tttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACg                                                                  >  2:100141/1‑93 (MQ=255)

                        ttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGc                                                                 <  2:262263/93‑1 (MQ=255)

                          gCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCtt                                                               <  1:2842483/93‑1 (MQ=255)

                          gCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCtt                                                               <  2:3316039/93‑1 (MQ=255)

                           ccGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTc                                                              <  1:3037951/93‑1 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            <  2:4455237/93‑1 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            <  2:1877756/93‑1 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            >  2:4680161/1‑93 (MQ=255)

                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           <  1:3937238/93‑1 (MQ=255)

                               aGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATg                                                          <  1:3118966/93‑1 (MQ=255)

                                ggggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATgg                                                         <  1:3204258/93‑1 (MQ=255)

                                 gggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgc                                                        >  1:1587523/1‑93 (MQ=255)

                                   gAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcag                                                      <  1:2490009/93‑1 (MQ=255)

                                   gAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcag                                                      <  2:586371/93‑1 (MQ=255)

                                    aaaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagc                                                     >  2:3465088/1‑93 (MQ=255)

                                    aaaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagc                                                     >  1:248320/1‑93 (MQ=255)

                                    aaaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagc                                                     <  1:3358860/93‑1 (MQ=255)

                                     aaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagca                                                    <  2:4510961/93‑1 (MQ=255)

                                          attattACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCaa                                               >  1:3703986/1‑93 (MQ=255)

                                             attaCAGAGGCCCAAAAAACAAAACCACCGGACTCTCATCCAGGGCTCTGTGGTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAttt                                            <  1:533323/93‑1 (MQ=255)

                                             attaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAttt                                            >  2:2520152/1‑93 (MQ=255)

                                              ttaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttt                                           <  2:3974730/93‑1 (MQ=255)

                                               taCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAttttt                                          <  2:2107868/93‑1 (MQ=255)

                                                 cAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTc                                        <  1:1686455/93‑1 (MQ=255)

                                                  agagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCa                                       >  2:2905595/1‑93 (MQ=255)

                                                   gagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCAt                                      >  1:4038552/1‑93 (MQ=255)

                                                    agGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATc                                     <  1:1383355/93‑1 (MQ=255)

                                                    agGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATc                                     <  2:1418844/93‑1 (MQ=255)

                                                     ggCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcg                                    >  2:1117129/1‑93 (MQ=255)

                                                      gCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcgc                                   >  1:6128/1‑93 (MQ=255)

                                                      gCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcgc                                   >  2:589800/1‑93 (MQ=255)

                                                        ccAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGt                                 >  2:238298/1‑93 (MQ=255)

                                                        ccAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGt                                 <  2:3012064/93‑1 (MQ=255)

                                                          aaaaaaCAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTc                               <  2:3260517/93‑1 (MQ=255)

                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              >  1:4209532/1‑93 (MQ=255)

                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              <  1:1140147/93‑1 (MQ=255)

                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              <  2:3497209/93‑1 (MQ=255)

                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              <  1:2055187/93‑1 (MQ=255)

                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              <  1:3993396/93‑1 (MQ=255)

                                                            aaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtt                             >  2:293491/1‑93 (MQ=255)

                                                            aaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtt                             >  2:2702175/1‑93 (MQ=255)

                                                             aaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCttt                            <  2:192220/93‑1 (MQ=255)

                                                             aaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCttt                            >  2:3881346/1‑93 (MQ=255)

                                                             aaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCttt                            <  2:3131228/93‑1 (MQ=255)

                                                               acaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTc                          <  2:3773311/93‑1 (MQ=255)

                                                                caaaaacCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTcc                         >  2:1592178/1‑93 (MQ=255)

                                                                caaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTcc                         >  2:912938/1‑93 (MQ=255)

                                                                caaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTcc                         <  1:3161808/93‑1 (MQ=255)

                                                                  aaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAg                       <  2:1205894/93‑1 (MQ=255)

                                                                  aaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAg                       >  1:4395185/1‑93 (MQ=255)

                                                                    aacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCTAGCt                     >  1:3603341/1‑93 (MQ=255)

                                                                    aacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCt                     <  1:2687636/93‑1 (MQ=255)

                                                                    aacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCt                     <  1:89096/93‑1 (MQ=255)

                                                                     acCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTg                    >  1:197664/1‑93 (MQ=255)

                                                                     acCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTg                    <  1:2400729/93‑1 (MQ=255)

                                                                       cccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGc                  <  2:2726742/93‑1 (MQ=255)

                                                                       cccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGc                  >  1:4720386/1‑93 (MQ=255)

                                                                        ccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCCATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGca                 >  1:4351147/1‑92 (MQ=255)

                                                                        ccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCg                 >  1:3825905/1‑93 (MQ=255)

                                                                         cGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGt                >  1:4235321/1‑93 (MQ=255)

                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               >  1:2257788/1‑93 (MQ=255)

                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               <  2:3682478/93‑1 (MQ=255)

                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               <  1:738669/93‑1 (MQ=255)

                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              <  2:1409140/93‑1 (MQ=255)

                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              >  2:4099779/1‑93 (MQ=255)

                                                                             ctctCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTacac            <  2:524013/93‑1 (MQ=255)

                                                                              tctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACg           >  2:3074486/1‑93 (MQ=255)

                                                                              tctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACg           >  2:740448/1‑93 (MQ=255)

                                                                               ctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGc          <  1:4601043/93‑1 (MQ=255)

                                                                                tcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCt         >  2:2426379/1‑93 (MQ=255)

                                                                                  aTCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCa       <  1:3270103/93‑1 (MQ=255)

                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      >  2:3607761/1‑93 (MQ=255)

                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      >  1:2602398/1‑93 (MQ=255)

                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACACTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      >  2:2332980/1‑93 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  2:1164826/1‑93 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    <  1:4112730/93‑1 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    <  2:1704786/93‑1 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  1:130520/1‑93 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  1:1033115/1‑93 (MQ=255)

                                                                                       gggCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGa  >  1:1915725/1‑93 (MQ=255)

                                                                                       gggCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGa  >  1:185722/1‑93 (MQ=255)

ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628

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