BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,953,377	T→C	T217A (ACG→GCG)	recB ←	RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,953,377	0	T	C	99.4%	970.7 / ‑3.5	311	T217A (ACG→GCG)	recB	RecB
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (155/153); total (156/155)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 322 total aligned reads displayed.

TCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCC  >  NC_000913/2953287‑2953466

tCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGctt                                                                                         >  2:3513097/1‑93 (MQ=255)

tCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGctt                                                                                         >  1:4159567/1‑93 (MQ=255)

tCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGctt                                                                                         <  1:1048645/93‑1 (MQ=255)

  aGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtca                                                                                       >  2:575860/1‑93 (MQ=255)

  aGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtca                                                                                       >  1:2017964/1‑93 (MQ=255)

    cgcATCCAGTCCCCCCACTGCTTCGCCCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtcatc                                                                                     <  1:2909487/93‑1 (MQ=255)

     gcATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtcatca                                                                                    >  2:7105058/1‑93 (MQ=255)

       aTCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCatcatc                                                                                  <  2:3087085/93‑1 (MQ=255)

       aTCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCatcatc                                                                                  >  1:2107239/1‑93 (MQ=255)

       aTCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTCcatcatt                                                                                  <  2:4798602/93‑2 (MQ=255)

        tCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCg                                                                                 >  1:5747086/1‑93 (MQ=255)

         ccAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCgg                                                                                <  2:2946079/93‑1 (MQ=255)

          cAGTTCACCCACTGCGGCGCGCCACTTCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCggg                                                                               >  2:7617806/1‑93 (MQ=255)

           aGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGggc                                                                              <  2:3228957/93‑1 (MQ=255)

            gTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCGTCATCATCGggcg                                                                             >  1:7023718/1‑93 (MQ=255)

            gTTCACCCACTGCGTCGCACCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGggcg                                                                             >  1:458098/1‑93 (MQ=255)

             ttCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGggcgg                                                                            >  2:1777871/1‑93 (MQ=255)

             ttCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGggcgg                                                                            >  1:6925719/1‑93 (MQ=255)

               cACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGgcggcg                                                                          <  2:6507391/93‑1 (MQ=255)

                 cccACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGt                                                                        <  1:2661383/93‑1 (MQ=255)

                 cccACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGt                                                                        >  2:6026847/1‑93 (MQ=255)

                 cccACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGt                                                                        >  1:4367663/1‑93 (MQ=255)

                   cACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGc                                                                      <  1:6132419/93‑1 (MQ=255)

                   cACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGc                                                                      <  1:664537/93‑1 (MQ=255)

                     cTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCtt                                                                    >  1:2780640/1‑93 (MQ=255)

                       gCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTg                                                                  >  2:7893863/1‑93 (MQ=255)

                        cGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGa                                                                 <  1:3677338/93‑1 (MQ=255)

                         gTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGAt                                                                >  2:2269313/1‑93 (MQ=255)

                           cgcgCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATaa                                                              >  2:2859858/1‑93 (MQ=255)

                            gcgcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAAc                                                             <  1:3389016/93‑1 (MQ=255)

                              gcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCg                                                           <  2:3117588/93‑1 (MQ=255)

                              gcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCg                                                           <  2:1238683/93‑1 (MQ=255)

                                cACTGCTGTTTTACCTTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGgc                                                         <  2:3308674/93‑1 (MQ=255)

                                cACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGgc                                                         >  2:7442840/1‑93 (MQ=255)

                                 aCTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCgtag                                                        >  2:5407043/1‑90 (MQ=255)

                                 aCTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGgcg                                                        >  1:467366/1‑93 (MQ=255)

                                  ctgctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGgcgt                                                       >  2:5903162/1‑92 (MQ=255)

                                  ctgctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGgcgc                                                       >  1:4677998/1‑93 (MQ=255)

                                    gctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCtt                                                     >  1:5974549/1‑93 (MQ=255)

                                    gctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCtt                                                     <  2:726778/93‑1 (MQ=255)

                                     ctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTc                                                    <  1:5836254/93‑1 (MQ=255)

                                      tgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCg                                                   >  2:1093100/1‑93 (MQ=255)

                                       gTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGc                                                  >  2:508137/1‑93 (MQ=255)

                                        ttttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGcc                                                 <  1:5432663/93‑1 (MQ=255)

                                        ttttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGcc                                                 <  1:7540643/93‑1 (MQ=255)

                                          ttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGCTAACCGGCGCTTCGCCtt                                               >  2:4162815/1‑93 (MQ=255)

                                          ttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCtt                                               <  2:1911179/93‑1 (MQ=255)

                                             ccGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCa                                            <  2:3166421/93‑1 (MQ=255)

                                             ccGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCa                                            >  1:4827871/1‑93 (MQ=255)

                                             ccGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCa                                            <  2:118511/93‑1 (MQ=255)

                                               gTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGa                                          >  1:3288543/1‑93 (MQ=255)

                                               gTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGa                                          <  2:4957652/93‑1 (MQ=255)

                                                tATCAATACGCGCCACAATTTGCGCATGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGAt                                         >  2:6240272/1‑93 (MQ=255)

                                                  tCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATaa                                       <  1:1865516/93‑1 (MQ=255)

                                                   cAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAAc                                      <  2:1583172/93‑1 (MQ=255)

                                                    aaTACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACg                                     >  2:1481358/1‑93 (MQ=255)

                                                    aaTACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACg                                     >  1:3584802/1‑93 (MQ=255)

                                                    aaTACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACg                                     >  1:4961182/1‑93 (MQ=255)

                                                    aaTACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTGCATCATCGGGCGGCGGGGCTTTGATAACCGGGGCTTCGGCTTGCAGATAACg                                     >  2:2201361/1‑93 (MQ=255)

                                                     aTACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGa                                    <  1:7176810/93‑1 (MQ=255)

                                                     aTACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGa                                    >  2:6934991/1‑93 (MQ=255)

                                                      tACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGAt                                   <  2:6655576/93‑1 (MQ=255)

                                                      tACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGAt                                   <  1:2314855/93‑1 (MQ=255)

                                                       aCGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGAtt                                  <  2:4523835/93‑1 (MQ=255)

                                                        cgcgCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTa                                 <  2:797969/93‑1 (MQ=255)

                                                         gcgcCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTaa                                >  1:2373944/1‑93 (MQ=255)

                                                           gcCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAata                              <  2:537925/93‑1 (MQ=255)

                                                           gcCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAata                              <  1:549551/93‑1 (MQ=255)

                                                           gcCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAata                              >  1:3525559/1‑93 (MQ=255)

                                                             cacaATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATc                            >  1:1516880/1‑93 (MQ=255)

                                                             cacaATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATc                            <  2:3818255/93‑1 (MQ=255)

                                                               caATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGGTTCCCCTTGCAGATAATGCTTAATATcgc                          >  1:6806883/1‑93 (MQ=255)

                                                               caATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCCGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATagc                          >  1:7789676/1‑93 (MQ=255)

                                                                 aTTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCgcgc                        <  1:2450140/93‑1 (MQ=255)

                                                                  tttGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgca                       <  1:3850644/93‑1 (MQ=255)

                                                                  tttGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgca                       >  1:4284063/1‑93 (MQ=255)

                                                                  tgtGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgca                       <  1:7362290/91‑1 (MQ=255)

                                                                   ttGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgcag                      <  2:5815498/93‑1 (MQ=255)

                                                                    tGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgcagc                     <  2:5597300/93‑1 (MQ=255)

                                                                      cgcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCaa                   <  1:21391/93‑1 (MQ=255)

                                                                      cgcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCaa                   <  1:2539195/93‑1 (MQ=255)

                                                                      cgcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCaa                   <  1:3445480/93‑1 (MQ=255)

                                                                      cgcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCaa                   <  1:5110750/93‑1 (MQ=255)

                                                                       gcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTTATATCGCGCAGCAAc                  >  2:432081/1‑93 (MQ=255)

                                                                       gcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAAc                  <  2:5328104/93‑1 (MQ=255)

                                                                       gcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAAc                  <  1:3096820/93‑1 (MQ=255)

                                                                       gcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAAc                  <  2:2425830/93‑1 (MQ=255)

                                                                          tGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGcc               >  2:6431630/1‑93 (MQ=255)

                                                                          tGACGGGAAGCCAGCGCTTCATCATCGGGCGGAGGTGCTTTGATAACCGGCGCTTCGGCTTGCATATAACGATTAATATCGCGCAGCAACGcc               >  2:1285443/1‑93 (MQ=255)

                                                                             cGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGc            <  1:2214820/93‑1 (MQ=255)

                                                                             cGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGc            <  1:1398118/93‑1 (MQ=255)

                                                                              gggAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCg           <  1:1222744/93‑1 (MQ=255)

                                                                               ggAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCgg          >  1:6542183/1‑93 (MQ=255)

                                                                               ggAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCgg          <  1:1650726/93‑1 (MQ=255)

                                                                                gAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGc         >  2:6759924/1‑93 (MQ=255)

                                                                                  aGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGGCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGccc       >  2:388015/1‑93 (MQ=255)

                                                                                     cAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCttt    >  2:4656807/1‑93 (MQ=255)

                                                                                      aGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTc   <  2:2065194/93‑1 (MQ=255)

                                                                                      aGCGCTTCATCATCGGGCGGCGATGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTc   >  1:1340565/1‑93 (MQ=255)

                                                                                       gcgcTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTcc  >  1:7222596/1‑93 (MQ=255)

TCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCC  >  NC_000913/2953287‑2953466

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