BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,975,207	T→C	intergenic (‑44/+83)	motA ← / ← flhC	MotA/FlhC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,975,207	0	T	C	99.2%	387.7 / ‑4.5	127	intergenic (‑44/+83)	motA/flhC	MotA/FlhC
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (63/62); total (64/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 138 total aligned reads displayed.

CCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACA  >  NC_000913/1975117‑1975296

cctaaaaCGGTACCGAGAACAACCATGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCa                                                                                         <  2:1755220/90‑1 (MQ=255)

ccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCa                                                                                         <  2:1833637/93‑1 (MQ=255)

 cgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAg                                                                                        <  1:556259/93‑1 (MQ=255)

   aaaaCTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGc                                                                                      >  2:925165/1‑93 (MQ=255)

   aaaaCTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGc                                                                                      >  1:182272/1‑93 (MQ=255)

     aaCTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGt                                                                                    >  2:1192122/1‑93 (MQ=255)

     aaCTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGt                                                                                    <  1:4378172/93‑1 (MQ=255)

      aCTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGTATGTTCAGTCGCCAGGCGtt                                                                                   <  2:1085670/93‑1 (MQ=255)

      aCTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGtt                                                                                   >  1:238522/1‑93 (MQ=255)

       cTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTa                                                                                  >  2:2558658/1‑93 (MQ=255)

        tGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTaa                                                                                 <  2:3883701/93‑1 (MQ=255)

        tGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTaa                                                                                 <  1:761748/93‑1 (MQ=255)

        tGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTaa                                                                                 <  1:3382690/93‑1 (MQ=255)

        tGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTaa                                                                                 >  2:3101773/1‑93 (MQ=255)

         gTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAAc                                                                                <  1:4280323/93‑1 (MQ=255)

          tACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAAcg                                                                               >  1:2403135/1‑93 (MQ=255)

           aCCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAAcgc                                                                              >  2:710860/1‑93 (MQ=255)

               agaACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGa                                                                          <  2:1262071/93‑1 (MQ=255)

                gaACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGAt                                                                         >  2:36654/1‑93 (MQ=255)

                gaACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGAt                                                                         >  1:356350/1‑93 (MQ=255)

                gaACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGAt                                                                         >  1:3106723/1‑93 (MQ=255)

                gaACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGAt                                                                         <  1:2216413/93‑1 (MQ=255)

                 aacaccCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGAtt                                                                        <  2:3738770/93‑1 (MQ=255)

                 aacaacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGAtt                                                                        <  1:2216698/93‑1 (MQ=255)

                    aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTggg                                                                     <  2:3951607/93‑1 (MQ=255)

                    aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTggg                                                                     <  2:1600277/93‑1 (MQ=255)

                    aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTggg                                                                     >  2:2069061/1‑93 (MQ=255)

                    aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTggg                                                                     >  2:1847689/1‑93 (MQ=255)

                     acCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTgggg                                                                    >  1:11941/1‑93 (MQ=255)

                     acCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTgggg                                                                    >  2:3436121/1‑93 (MQ=255)

                         ggTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCaaa                                                                <  1:4479131/93‑1 (MQ=255)

                             aCCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCaaaaaaa                                                            >  2:4416497/1‑93 (MQ=255)

                             aCCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCaaaaaaa                                                            >  2:4416498/1‑93 (MQ=255)

                               cTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAag                                                          >  1:9936/1‑93 (MQ=255)

                                taataaGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAagc                                                         >  1:3856939/1‑93 (MQ=255)

                                taataaGATAAGCACGACATCATCCTTCCACTGCTGCCCATAACAGGATGTTCAGTCGCCAGCCGTTAACGCGCGATTGGGGCAAAAAAAagc                                                         <  1:1823251/93‑1 (MQ=255)

                                 aataaGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAagca                                                        <  2:3359089/93‑1 (MQ=255)

                                     agataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCgg                                                    >  2:4592741/1‑93 (MQ=255)

                                      gataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGt                                                   >  2:2338820/1‑93 (MQ=255)

                                       ataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTa                                                  <  1:3993522/93‑1 (MQ=255)

                                        taagCACGACATCATCCGTCCACTGGTGACCATGACAGGATGGTCAGTCGGCAGGCGTTAACGCGCGATTTGGGCAGAAAGTAGCAGCGGTag                                                 >  2:1296057/1‑92 (MQ=255)

                                          agCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAcgt                                               >  2:2258353/1‑93 (MQ=255)

                                          agCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAcgt                                               <  1:2582102/93‑1 (MQ=255)

                                             aCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGTGGCAAACAAAAGCAGCGGTAcgtcgt                                            >  2:2859600/1‑93 (MQ=255)

                                             aCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAcgtcgt                                            <  2:1129441/93‑1 (MQ=255)

                                             aCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAcgtcgt                                            >  2:1760959/1‑93 (MQ=255)

                                              cGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGtt                                           <  2:389002/93‑1 (MQ=255)

                                              cGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGtt                                           <  1:2552301/93‑1 (MQ=255)

                                                aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTAc                                         <  1:3879287/93‑1 (MQ=255)

                                                aCATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTAc                                         <  1:4251086/93‑1 (MQ=255)

                                                   tcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgc                                      >  1:3015463/1‑93 (MQ=255)

                                                   tcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgc                                      <  1:1367231/93‑1 (MQ=255)

                                                     atcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgctg                                    >  2:1428234/1‑93 (MQ=255)

                                                     atcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGAGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgctc                                    >  1:4381589/1‑92 (MQ=255)

                                                      tgcTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgctgc                                   >  2:2958272/3‑93 (MQ=255)

                                                      tcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgctgc                                   <  1:1459684/93‑1 (MQ=255)

                                                        cTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGctgctg                                 >  2:2012792/1‑93 (MQ=255)

                                                        cTTCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGctgctg                                 >  1:4620770/1‑93 (MQ=255)

                                                         ttCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTgg                                >  2:932787/1‑93 (MQ=255)

                                                          tCCACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGa                               >  2:803425/1‑93 (MQ=255)

                                                           ccACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGaa                              >  2:4471417/1‑93 (MQ=255)

                                                            cACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAAt                             >  1:2903799/1‑93 (MQ=255)

                                                            cACTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAAt                             >  1:431314/1‑93 (MQ=255)

                                                             aCTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATg                            <  1:1719249/93‑1 (MQ=255)

                                                             aCTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATg                            <  2:3488452/93‑1 (MQ=255)

                                                              cTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGCATGt                           >  1:3284546/1‑93 (MQ=255)

                                                              cTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGt                           >  2:1126357/1‑93 (MQ=255)

                                                              cTGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGt                           <  1:1050965/93‑1 (MQ=255)

                                                               tGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGtt                          >  1:569824/1‑93 (MQ=255)

                                                               tGTTGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGtt                          >  1:2113634/1‑93 (MQ=255)

                                                                 ttGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTgc                        >  1:1026977/1‑93 (MQ=255)

                                                                 ttGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTgc                        <  1:4640639/93‑1 (MQ=255)

                                                                  tGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACAGCTGCTGGAATGTTgcg                       >  1:832306/1‑93 (MQ=255)

                                                                  tGACCATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGC‑AAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTgcgc                      <  2:875375/93‑1 (MQ=255)

                                                                     ccATGACAGGATGTTCAGTGGCCAGGCGTTGACGCGCGATTGGGCCAAAAAAAAGCAGCGGTACGTCGTTCCCTCTGCTGGAATGTTGCGCCt                    <  2:279261/93‑1 (MQ=255)

                                                                      cATGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTc                   >  1:3510618/1‑93 (MQ=255)

                                                                       aTGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCa                  <  1:2404806/93‑1 (MQ=255)

                                                                       aTGACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCa                  >  1:3839264/1‑93 (MQ=255)

                                                                         gACAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAcc                <  2:4424525/93‑1 (MQ=255)

                                                                          aCAGGATGTTCAGTCGCCAGGCGTTTACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg               >  2:4618854/1‑93 (MQ=255)

                                                                          aCAGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg               <  2:3718003/93‑1 (MQ=255)

                                                                            aGGATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTa             <  2:4102873/93‑1 (MQ=255)

                                                                             ggATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTAt            <  2:2226409/93‑1 (MQ=255)

                                                                              gATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATc           >  1:2127259/1‑93 (MQ=255)

                                                                              gATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATc           >  1:1488425/1‑93 (MQ=255)

                                                                              gATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATc           >  1:3824757/1‑93 (MQ=255)

                                                                              gATGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATc           >  1:1794073/1‑93 (MQ=255)

                                                                               aTGTTCATTCGCCAGGCGTTACCGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa          <  2:1742127/93‑1 (MQ=255)

                                                                               aTGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa          <  1:327898/93‑1 (MQ=255)

                                                                                tGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGGGCCTCACCGTATCAg         >  2:418642/1‑93 (MQ=255)

                                                                                tGTTCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg         <  2:477040/93‑1 (MQ=255)

                                                                                  ttCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt       >  2:2682673/1‑93 (MQ=255)

                                                                                   tCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTa      <  1:2121759/93‑1 (MQ=255)

                                                                                   tCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTa      <  2:663604/93‑1 (MQ=255)

                                                                                   tCAGTCGCCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTa      <  1:1318852/93‑1 (MQ=255)

                                                                                    cAGTCGCCAGGCGTTAACGCGCGATTGGGGC‑AAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaaa    <  1:2968127/93‑1 (MQ=255)

                                                                                      gtcgccAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAc   >  2:3486155/1‑93 (MQ=255)

                                                                                      gtcgccAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAc   <  1:3930014/93‑1 (MQ=255)

                                                                                      gtcgccAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAc   <  1:3870272/93‑1 (MQ=255)

                                                                                       tcgccAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACa  <  2:1873920/93‑1 (MQ=255)

CCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACA  >  NC_000913/1975117‑1975296

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