BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,015,850	T→C	Q558R (CAG→CGG)	b0955 ←	putative ATP‑dependent protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,015,850	0	T	C	100.0%	891.5 / NA	252	Q558R (CAG→CGG)	b0955	putative ATP‑dependent protease
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (130/122); total (130/122)

Only 100 of 254 total aligned reads displayed.

GATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCC  >  NC_000913/1015760‑1015939

gATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCgg                                                                                         >  2:216561/1‑93 (MQ=255)

gATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCgg                                                                                         >  2:3207502/1‑93 (MQ=255)

gATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCgg                                                                                         >  1:1949523/1‑93 (MQ=255)

gATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCgg                                                                                         >  2:4383439/1‑93 (MQ=255)

gATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCgg                                                                                         >  2:419232/1‑93 (MQ=255)

 aTCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGa                                                                                        >  1:6356437/1‑93 (MQ=255)

 aTCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGa                                                                                        >  1:5521737/1‑93 (MQ=255)

  tCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGAt                                                                                       >  1:881920/1‑93 (MQ=255)

    aGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATgg                                                                                     >  2:4060110/1‑93 (MQ=255)

     gTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGt                                                                                    >  2:3629976/1‑93 (MQ=255)

     gTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGt                                                                                    >  1:6094853/1‑93 (MQ=255)

     gTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGt                                                                                    >  2:1632810/1‑93 (MQ=255)

      ttCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGtt                                                                                   >  1:1594466/1‑93 (MQ=255)

      ttCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGtt                                                                                   <  1:919729/93‑1 (MQ=255)

       tCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGttt                                                                                  >  2:5215241/1‑93 (MQ=255)

       tCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGttt                                                                                  >  2:2559767/1‑93 (MQ=255)

          gAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCa                                                                               >  2:274779/1‑93 (MQ=255)

          gAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCa                                                                               >  1:489651/1‑93 (MQ=255)

            aTAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATc                                                                             >  1:5399067/1‑93 (MQ=255)

             tAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCa                                                                            >  2:564148/1‑93 (MQ=255)

             tAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCa                                                                            >  2:291242/1‑93 (MQ=255)

             tAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCa                                                                            >  2:1795211/1‑93 (MQ=255)

              aaaCCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAg                                                                           <  1:308833/93‑1 (MQ=255)

              aaaCCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAg                                                                           >  1:1235079/1‑93 (MQ=255)

               aaCCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGc                                                                          <  2:3153602/93‑1 (MQ=255)

               aaCCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTCGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGc                                                                          <  1:4136609/93‑1 (MQ=255)

                 ccAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGcgt                                                                        <  2:1679110/93‑1 (MQ=255)

                 ccAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGcgt                                                                        <  1:4602773/93‑1 (MQ=255)

                 ccAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGcgt                                                                        <  2:5949480/93‑1 (MQ=255)

                    gTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGcgtcgt                                                                     >  2:834449/1‑93 (MQ=255)

                     ttCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGtt                                                                    >  1:748408/1‑93 (MQ=255)

                      tCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGttt                                                                   <  1:4923914/93‑1 (MQ=255)

                      tCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGttt                                                                   <  1:2524744/93‑1 (MQ=255)

                        aGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTgg                                                                 <  2:6313280/93‑1 (MQ=255)

                          ccAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGcc                                                               <  1:5259500/93‑1 (MQ=255)

                          ccAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGcc                                                               <  2:1074857/93‑1 (MQ=255)

                          ccAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGcc                                                               <  2:755629/93‑1 (MQ=255)

                            aaCGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCtt                                                             <  2:5324604/93‑1 (MQ=255)

                            aaCGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCtt                                                             <  1:3597186/93‑1 (MQ=255)

                             aCGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTc                                                            <  2:2063002/93‑1 (MQ=255)

                              cGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCg                                                           <  1:5731323/93‑1 (MQ=255)

                                 aaTGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCa                                                        >  1:1373624/1‑93 (MQ=255)

                                 aaTGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCa                                                        <  2:4528171/93‑1 (MQ=255)

                                  aTGCCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCAt                                                       <  2:6348849/93‑1 (MQ=255)

                                    ggCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATcc                                                     <  2:1932759/93‑1 (MQ=255)

                                     gCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATccc                                                    <  2:6081922/93‑1 (MQ=255)

                                       cATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCcac                                                  >  1:4746901/1‑93 (MQ=255)

                                        aTGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCcaca                                                 >  1:37173/1‑93 (MQ=255)

                                        aTGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCcaca                                                 >  2:176600/1‑93 (MQ=255)

                                         tGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCacac                                                >  2:4728503/1‑93 (MQ=255)

                                         tGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCacac                                                <  1:5863232/93‑1 (MQ=255)

                                          ggAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACAcc                                               <  1:2657946/93‑1 (MQ=255)

                                          ggAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACAcc                                               >  1:5946510/1‑93 (MQ=255)

                                           gAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCa                                              <  1:3128974/93‑1 (MQ=255)

                                            aaaaCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAg                                             >  1:953006/1‑93 (MQ=255)

                                             aaaCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGTATCAGCGTCGTTTGGCCTTCGCCATCCCACGCCAGa                                            <  1:2105044/93‑1 (MQ=255)

                                             aaaCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGa                                            <  1:1470293/93‑1 (MQ=255)

                                             aaaCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGa                                            <  2:2255493/93‑1 (MQ=255)

                                              agcGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGAt                                           <  1:724387/91‑1 (MQ=255)

                                              aaCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGAt                                           >  1:4066097/1‑93 (MQ=255)

                                              aaCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGAt                                           >  1:179417/1‑93 (MQ=255)

                                                 ggTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTta                                        >  2:4340231/1‑93 (MQ=255)

                                                 ggTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTta                                        >  1:2965751/1‑93 (MQ=255)

                                                  gTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTAGCCATCCCACACCAGATTtaa                                       <  1:5309120/93‑1 (MQ=255)

                                                     aCGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaataa                                    <  2:2391612/93‑1 (MQ=255)

                                                     aCGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCAATCCCACACCAGATTtaataa                                    >  1:4383760/1‑93 (MQ=255)

                                                       gTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaataata                                  >  2:4610558/1‑93 (MQ=255)

                                                        tccttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaataataa                                 >  2:2813864/1‑93 (MQ=255)

                                                        tccttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaataataa                                 <  2:3817502/93‑1 (MQ=255)

                                                        tccttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaataataa                                 >  2:5574533/1‑93 (MQ=255)

                                                        tccttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaataataa                                 >  2:5715988/1‑93 (MQ=255)

                                                        tccttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTCTGGCCTTCGCCATCCCACACCAGATTtaataataa                                 <  2:5120903/93‑1 (MQ=255)

                                                         ccttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAAc                                <  2:3482233/93‑1 (MQ=255)

                                                         ccttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAAc                                <  1:2391228/93‑1 (MQ=255)

                                                          cttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACg                               <  2:3078241/93‑1 (MQ=255)

                                                          cttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACg                               >  1:5228525/1‑93 (MQ=255)

                                                           ttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACgg                              >  2:5062731/1‑93 (MQ=255)

                                                           ttcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACgg                              <  2:3741000/93‑1 (MQ=255)

                                                            tcctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATTATAACGGc                             <  1:6299361/93‑1 (MQ=255)

                                                             cctGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCa                            <  2:1201455/93‑1 (MQ=255)

                                                               tGTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGt                          >  1:1018416/1‑93 (MQ=255)

                                                                gTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTg                         >  2:4583617/1‑93 (MQ=255)

                                                                gTTGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTg                         >  1:3522311/1‑93 (MQ=255)

                                                                 ttGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGc                        <  1:2884975/93‑1 (MQ=255)

                                                                  tGCGATGCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCg                       <  2:3378460/93‑1 (MQ=255)

                                                                     gATACTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGAGCgtca                    >  1:3260213/1‑93 (MQ=255)

                                                                        gCTTGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCgtcagtc                 <  2:2618499/93‑1 (MQ=255)

                                                                           tGCGCGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTcaca              >  1:1283053/1‑93 (MQ=255)

                                                                            gcgcGATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACAt             <  2:68508/93‑1 (MQ=255)

                                                                             cgcgATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATc            >  1:1833404/1‑93 (MQ=255)

                                                                               cgATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGt          <  2:4660120/93‑1 (MQ=255)

                                                                               cgATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGt          <  2:1325397/93‑1 (MQ=255)

                                                                                gATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTc         >  2:1317263/1‑93 (MQ=255)

                                                                                gATACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTc         <  1:3808097/93‑1 (MQ=255)

                                                                                 aTACGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCt        <  1:564148/93‑1 (MQ=255)

                                                                                   aCGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTAc      <  2:5424331/93‑1 (MQ=255)

                                                                                   aCGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTAc      <  2:5397938/93‑1 (MQ=255)

                                                                                   aCGTTCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTAc      <  1:6172041/93‑1 (MQ=255)

                                                                                      ttCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGc   <  2:326681/93‑1 (MQ=255)

                                                                                       tCCCGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGcc  >  1:3906890/1‑93 (MQ=255)

GATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCC  >  NC_000913/1015760‑1015939

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