BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	651,821	A→G	R12R (CGT→CGC)	citC ←	citrate lyase synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	651,821	0	A	G	100.0%	1711.4 / NA	408	R12R (CGT→CGC)	citC	citrate lyase synthetase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (206/202); total (206/202)

Only 100 of 408 total aligned reads displayed.

TCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTG  >  NC_000913/651723‑651918

tCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTGATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCg                                                                                                 >  1:1304087/1‑101 (MQ=255)

tCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCg                                                                                                 <  2:1303204/101‑1 (MQ=255)

  gcgGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGtt                                                                                               >  1:1304092/1‑101 (MQ=255)

   cgGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGttt                                                                                              >  2:1304097/1‑101 (MQ=255)

    gggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGtttt                                                                                             <  2:1303348/101‑1 (MQ=255)

    gggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGtttt                                                                                             <  1:1303758/101‑1 (MQ=255)

    gggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGtttt                                                                                             >  1:1304101/1‑101 (MQ=255)

     ggTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTa                                                                                            <  1:1303407/101‑1 (MQ=255)

       ttACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTAcg                                                                                          >  2:1304105/1‑101 (MQ=255)

        tACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTAcgc                                                                                         >  1:1304108/1‑101 (MQ=255)

           ggTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTa                                                                                            <  2:1303380/95‑1 (MQ=255)

            gTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGt                                                                                     <  2:1303365/101‑1 (MQ=255)

            gTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGt                                                                                     >  1:1304115/1‑101 (MQ=255)

              aataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGa                                                                                   >  2:1304118/1‑101 (MQ=255)

               ataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGaa                                                                                  >  2:1304125/1‑101 (MQ=255)

               ataaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGaa                                                                                  >  1:1304122/1‑101 (MQ=255)

                taaataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGaaa                                                                                 <  1:1303609/101‑1 (MQ=255)

                  aataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAat                                                                               >  2:1304130/1‑101 (MQ=255)

                  aataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAat                                                                               <  1:1303130/101‑1 (MQ=255)

                   ataCTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAata                                                                              <  1:1303448/101‑1 (MQ=255)

                      cTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCa                                                                           >  1:1304134/1‑101 (MQ=255)

                       ttCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCAt                                                                          <  1:1303432/101‑1 (MQ=255)

                       ttCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCAt                                                                          <  1:1303392/101‑1 (MQ=255)

                        tCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCAtt                                                                         <  1:1303363/101‑1 (MQ=255)

                         cGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTg                                                                        >  1:1304140/1‑101 (MQ=255)

                          gACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGc                                                                       <  1:1303476/101‑1 (MQ=255)

                          gACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGc                                                                       >  2:1304145/1‑101 (MQ=255)

                            cTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCg                                                                     >  2:1304150/1‑101 (MQ=255)

                            cTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCg                                                                     >  1:1304149/1‑101 (MQ=255)

                               tggtgTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAAc                                                                  <  1:1303345/101‑1 (MQ=255)

                               tggtgTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAAc                                                                  <  2:1303760/101‑1 (MQ=255)

                                ggtgTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACa                                                                 >  2:1304156/1‑101 (MQ=255)

                                 gtgtCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACat                                                                <  1:1303561/101‑1 (MQ=255)

                                  tgtCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACata                                                               <  2:1303743/101‑1 (MQ=255)

                                  tgtCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACata                                                               <  1:1303550/101‑1 (MQ=255)

                                  tgtCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACata                                                               <  2:1303246/101‑1 (MQ=255)

                                   gtCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACataa                                                              >  1:1304161/1‑101 (MQ=255)

                                   gtCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACataa                                                              <  1:1303695/101‑1 (MQ=255)

                                     cAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACataata                                                            >  1:1304165/1‑101 (MQ=255)

                                      aaCGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCAt                                                                          <  2:1303464/86‑1 (MQ=255)

                                      aaCGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACAtaataa                                                           >  2:1304167/1‑101 (MQ=255)

                                         gCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaata                                                        >  1:1304175/1‑101 (MQ=255)

                                         gCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaata                                                        >  2:1304172/1‑101 (MQ=255)

                                           tCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGt                                                      >  2:1304179/1‑101 (MQ=255)

                                            cAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTa                                                     <  2:1303292/101‑1 (MQ=255)

                                             aaaTCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTAt                                                    >  2:1304184/1‑101 (MQ=255)

                                              aaTCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATc                                                   <  1:1303275/101‑1 (MQ=255)

                                               aTCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGaa                                                                   <  2:1303479/84‑1 (MQ=255)

                                                tCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                 >  2:1304188/1‑101 (MQ=255)

                                                tCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                 <  2:1303660/101‑1 (MQ=255)

                                                tCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                 <  2:1303490/101‑1 (MQ=255)

                                                tCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCt                                                 <  1:1303539/101‑1 (MQ=255)

                                                   ttttCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGaa                                              >  2:1304194/1‑101 (MQ=255)

                                                    tttCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                             >  1:1304195/1‑101 (MQ=255)

                                                    tttCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAg                                             <  2:1303507/101‑1 (MQ=255)

                                                      tCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGt                                           <  1:1303599/101‑1 (MQ=255)

                                                        aTGCAGGAATTGGGGGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGt                                           <  1:1303222/99‑1 (MQ=255)

                                                        aTGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGc                                         >  2:1304203/1‑101 (MQ=255)

                                                        aTGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGc                                         >  1:1304205/1‑101 (MQ=255)

                                                          gCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAt                                       >  2:1304208/1‑101 (MQ=255)

                                                           cAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTCTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtg                                      <  1:1303648/101‑1 (MQ=255)

                                                             ggAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgtt                                    >  2:1304216/1‑101 (MQ=255)

                                                             ggAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgtt                                    >  1:1304215/1‑101 (MQ=255)

                                                               aaTTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttgt                                  >  1:1304220/1‑101 (MQ=255)

                                                                agtGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCggg                                                                                      <  1:1303588/46‑1 (MQ=39)

                                                                aTTGGGCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATgttgtt                                 >  2:1304225/1‑101 (MQ=255)

                                                                    ggCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCg                             >  1:1304230/1‑101 (MQ=255)

                                                                     gCGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGa                            >  2:1304232/1‑101 (MQ=255)

                                                                      cGATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAt                           <  2:1303595/101‑1 (MQ=255)

                                                                       gATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAtt                          <  1:1303851/101‑1 (MQ=255)

                                                                       gATTTCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGAtt                          >  1:1304237/1‑101 (MQ=255)

                                                                         tttCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGt                                           <  2:1303438/82‑1 (MQ=255)

                                                                         tttCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTg                        >  2:1304243/1‑101 (MQ=255)

                                                                         tttCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTg                        <  1:1303678/101‑1 (MQ=255)

                                                                         tttCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTg                        <  1:1303636/101‑1 (MQ=255)

                                                                         tttCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTg                        <  1:1303492/101‑1 (MQ=255)

                                                                          ttCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGc                       >  1:1304246/1‑101 (MQ=255)

                                                                           tCCGCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCa                      >  2:1304248/1‑101 (MQ=255)

                                                                              gCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACg                   >  2:1304253/1‑101 (MQ=255)

                                                                              gCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACg                   >  1:1304256/1‑101 (MQ=255)

                                                                              gCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACg                   <  2:1303717/101‑1 (MQ=255)

                                                                              gCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACg                   <  1:1303802/101‑1 (MQ=255)

                                                                              gCCATTTTTTTATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACg                   <  1:1303529/101‑1 (MQ=255)

                                                                                  tttttttATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtg               >  2:1304261/1‑101 (MQ=255)

                                                                                    tttttATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgtt             <  2:1303528/101‑1 (MQ=255)

                                                                                    tttttATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgtt             <  2:1303848/101‑1 (MQ=255)

                                                                                      tttATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgttgt           <  2:1303559/101‑1 (MQ=255)

                                                                                      tttATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgttgt           >  2:1304268/1‑101 (MQ=255)

                                                                                      tttATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgttgt           >  1:1304267/1‑101 (MQ=255)

                                                                                       ttATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATgttgtt          <  2:1304032/101‑1 (MQ=255)

                                                                                        tATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTc         <  2:1303577/101‑1 (MQ=255)

                                                                                        tATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTc         >  1:1304277/1‑101 (MQ=255)

                                                                                        tATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTc         <  2:1303631/101‑1 (MQ=255)

                                                                                        tATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTc         >  2:1304274/1‑101 (MQ=255)

                                                                                        tATTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTc         <  2:1303783/101‑1 (MQ=255)

                                                                                         aTTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCa        <  2:1303412/101‑1 (MQ=255)

                                                                                         aTTTTCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCa        >  1:1304282/1‑101 (MQ=255)

                                                                                          ttttCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCaa       <  1:1303511/101‑1 (MQ=255)

                                                                                          ttttCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCaa       >  2:1304284/1‑101 (MQ=255)

                                                                                            ttCTGAGCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATg     >  1:1304291/1‑101 (MQ=255)

TCGCGGGTTACGGTAATAAATACTTCGACTGTGGTGTCAACGCTCAAATCATTTTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTG  >  NC_000913/651723‑651918

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