BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,785,435	T→C	T560A (ACT→GCT)	ppsA ←	phosphoenolpyruvate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,785,435	0	T	C	100.0%	899.6 / NA	219	T560A (ACT→GCT)	ppsA	phosphoenolpyruvate synthetase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (98/121); total (98/121)

Only 100 of 232 total aligned reads displayed.

ACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGC  >  NC_000913/1785337‑1785532

cccGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTtcagcca                                                                                                 <  2:2328422/100‑1 (MQ=255)

accGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTtcagcca                                                                                                 >  1:2328962/1‑101 (MQ=255)

accGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTtcagcca                                                                                                 <  1:2328612/101‑1 (MQ=255)

   gACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAgac                                                                                              <  1:2328880/101‑1 (MQ=255)

    aCCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAgacg                                                                                             >  1:2328965/1‑101 (MQ=255)

     ccAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAgacga                                                                                            <  1:2328463/101‑1 (MQ=255)

     ccAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAgacga                                                                                            >  2:2328967/1‑101 (MQ=255)

      cAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAGCGCGGCACCCAGCGTCGCGATCCCTTCAGCCAgacgac                                                                                           <  1:2328598/101‑1 (MQ=255)

      cAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAgacgac                                                                                           <  2:2328586/101‑1 (MQ=255)

       aGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGAcc                                                                                          >  1:2328969/1‑101 (MQ=255)

         gTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCaa                                                                                        >  1:2328974/1‑101 (MQ=255)

          ttGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAAc                                                                                       >  1:2328976/1‑101 (MQ=255)

           tGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACg                                                                                      >  2:2328977/1‑101 (MQ=255)

              cATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTaa                                                                                   >  2:2328980/1‑101 (MQ=255)

              cATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTaa                                                                                   <  1:2328457/101‑1 (MQ=255)

               atatTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTaaa                                                                                  >  1:2328984/1‑101 (MQ=255)

                tatTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTaaaa                                                                                 >  2:2328985/1‑101 (MQ=255)

                  ttcgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAAtt                                                                               >  1:2328987/1‑101 (MQ=255)

                  ttcgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAAtt                                                                               <  1:2328556/101‑1 (MQ=255)

                   tcgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTc                                                                              <  2:2328680/101‑1 (MQ=255)

                     gttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCAc                                                                            <  2:2328646/101‑1 (MQ=255)

                      ttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCAcg                                                                           >  1:2328991/1‑101 (MQ=255)

                       tcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCAcgc                                                                          <  1:2328790/101‑1 (MQ=255)

                        cgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCAcgcg                                                                         >  2:2328995/1‑101 (MQ=255)

                         gATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCgg                                                                        >  1:2328997/1‑101 (MQ=255)

                          aTTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGga                                                                       <  1:2328670/101‑1 (MQ=255)

                           tttAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGgag                                                                      <  2:2328744/101‑1 (MQ=255)

                           tttAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGgag                                                                      <  1:2328630/101‑1 (MQ=255)

                           tttAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGgag                                                                      <  1:2328658/101‑1 (MQ=255)

                            ttAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGgaga                                                                     >  2:2328998/1‑101 (MQ=255)

                              aaaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCaa                                                                                        <  1:2328578/80‑1 (MQ=255)

                               aaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAgacgac                                                                                           >  2:2329002/1‑76 (MQ=255)

                               aaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATc                                                                  >  1:2329001/1‑101 (MQ=255)

                               aaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATc                                                                  <  2:2328569/101‑1 (MQ=255)

                                 aTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCaa                                                                >  1:2329005/1‑101 (MQ=255)

                                   cAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCaaaa                                                              <  2:2328894/101‑1 (MQ=255)

                                     gagagaCGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAAcc                                                            <  1:2328622/101‑1 (MQ=255)

                                     gagagaCGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAAcc                                                            >  2:2329007/1‑101 (MQ=255)

                                      agagaCGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCt                                                           >  2:2329010/1‑101 (MQ=255)

                                       gagaCGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCtt                                                          <  1:2328749/101‑1 (MQ=255)

                                         gaCGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTa                                                                                    >  1:2329015/1‑73 (MQ=255)

                                         gaCGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTtc                                                        >  1:2329013/1‑101 (MQ=255)

                                         gaCGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTtc                                                        <  2:2328697/101‑1 (MQ=255)

                                           cGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTtcat                                                      <  2:2328683/101‑1 (MQ=255)

                                           cGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTtcat                                                      <  1:2328762/101‑1 (MQ=255)

                                            ggACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTtcatc                                                     <  2:2328517/101‑1 (MQ=255)

                                            ggACAATGACGCGCGTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGga                                                                       <  2:2328634/83‑1 (MQ=255)

                                               cAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCatcatc                                                  <  1:2328736/101‑1 (MQ=255)

                                               cAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCatcatc                                                  >  2:2329017/1‑101 (MQ=255)

                                                aaTGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCAtct                                                 >  1:2329018/1‑101 (MQ=255)

                                                  tGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCAtctc                                                >  2:2329023/1‑100 (MQ=255)

                                                     cgcgCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTtcagcca                                                                                                 <  1:2328522/48‑1 (MQ=255)

                                                     cgcgCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCgg                                            >  2:2329026/1‑101 (MQ=255)

                                                     cgcgCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCgg                                            >  1:2329025/1‑101 (MQ=255)

                                                     cgcgCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCgg                                            <  2:2328654/101‑1 (MQ=255)

                                                      gcgcTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGa                                           >  1:2329028/1‑101 (MQ=255)

                                                       cgcTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGAt                                          >  2:2329031/1‑101 (MQ=255)

                                                        gcTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTa                                                                                    >  1:2329034/1‑58 (MQ=255)

                                                        gcTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGAtt                                         <  2:2328615/101‑1 (MQ=255)

                                                          ttCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTc                                       >  2:2329036/1‑101 (MQ=255)

                                                          ttCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTc                                       <  1:2328852/101‑1 (MQ=255)

                                                          ttCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTc                                       <  2:2328649/101‑1 (MQ=255)

                                                            cGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGt                                     >  2:2329039/1‑101 (MQ=255)

                                                             ggATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCAcgc                                                                          <  2:2328599/63‑1 (MQ=255)

                                                             ggATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtt                                    >  1:2329040/1‑101 (MQ=255)

                                                              gATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGttt                                   >  2:2329043/1‑101 (MQ=255)

                                                               aTAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtttt                                  >  1:2329044/1‑101 (MQ=255)

                                                                tAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTc                                       <  1:2328704/95‑1 (MQ=255)

                                                                tAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTg                                 <  2:2328772/101‑1 (MQ=255)

                                                                tAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTg                                 <  2:2328788/101‑1 (MQ=255)

                                                                tAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTg                                 <  2:2328955/101‑1 (MQ=255)

                                                                 aaaaCGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGc                                >  2:2329047/1‑101 (MQ=255)

                                                                   aaCGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCaa                              <  1:2328590/101‑1 (MQ=255)

                                                                    aCGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAAc                             <  2:2328730/101‑1 (MQ=255)

                                                                    aCGCGGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAAc                             >  2:2329050/1‑101 (MQ=255)

                                                                       cgGCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGc                          <  1:2328575/101‑1 (MQ=255)

                                                                         gCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCgg                        >  1:2329051/1‑101 (MQ=255)

                                                                         gCACCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCgg                        <  2:2328555/101‑1 (MQ=255)

                                                                           aCCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGtt                      <  1:2328720/101‑1 (MQ=255)

                                                                           aCCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGtt                      <  2:2328795/101‑1 (MQ=255)

                                                                           aCCCAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGtt                      <  1:2328768/101‑1 (MQ=255)

                                                                            cccAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGga                                                                       <  1:2328302/51‑1 (MQ=255)

                                                                             ccAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCgg                                            <  1:2328668/77‑1 (MQ=255)

                                                                             ccAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTcc                    >  2:2329054/1‑101 (MQ=255)

                                                                             ccAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTcc                    >  1:2329055/1‑101 (MQ=255)

                                                                              cAGCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCt                   <  1:2328691/101‑1 (MQ=255)

                                                                                gCGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGa                 <  1:2328645/101‑1 (MQ=255)

                                                                                 cGTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGAt                <  1:2328839/101‑1 (MQ=255)

                                                                                  gTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATc               <  1:2328421/101‑1 (MQ=255)

                                                                                  gTCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATc               <  1:2328843/101‑1 (MQ=255)

                                                                                   tCGCGATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCg              <  1:2328809/101‑1 (MQ=255)

                                                                                     gcgATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGt             <  1:2328777/100‑1 (MQ=255)

                                                                                      cgATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCa           <  2:2328863/101‑1 (MQ=255)

                                                                                      cgATCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCa           <  2:2328709/101‑1 (MQ=255)

                                                                                        aTCCCTTCAGCCAGACGGCCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACt                            >  2:2329062/1‑82 (MQ=255)

                                                                                        aTCCCTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCaaa         >  2:2329060/1‑101 (MQ=255)

                                                                                         tCCCGTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATc               <  1:2328641/94‑1 (MQ=255)

                                                                                            cTTCAGCCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCAAACTCa     >  2:2329064/1‑101 (MQ=255)

                                                                                               cagccagACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGc  >  1:2329066/1‑101 (MQ=255)

                                                                                               cagccagACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGc  <  1:2328714/101‑1 (MQ=255)

ACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGC  >  NC_000913/1785337‑1785532

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