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.3%	456.1 / ‑5.1	143	T217A (ACG→GCG)	recB	RecB
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (70/71); total (71/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 152 total aligned reads displayed.

TCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCC  >  NC_000913/2953287‑2953466

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

 cAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtc                                                                                        <  1:4490769/93‑1 (MQ=255)

  aGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtca                                                                                       <  1:4428302/93‑1 (MQ=255)

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

   gcgcATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtcat                                                                                      <  2:3200570/93‑1 (MQ=255)

    cgcATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtcatc                                                                                     <  2:3746048/93‑1 (MQ=255)

      cATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGGAGCCAGCGCTTcatcat                                                                                   >  2:1709446/1‑93 (MQ=255)

      cATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTcatcat                                                                                   >  1:4645663/1‑93 (MQ=255)

      cATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTcatcat                                                                                   >  2:12209/1‑93 (MQ=255)

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

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

        tCCAGTTCCCCCGCTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCg                                                                                 <  2:862541/93‑1 (MQ=255)

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

              tCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGggcggc                                                                           <  2:3463366/93‑1 (MQ=255)

                ccccACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGcggcgg                                                                         <  1:1087820/92‑1 (MQ=255)

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

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

                    aCTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCt                                                                     >  2:1855464/1‑93 (MQ=255)

                    aCTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCt                                                                     >  2:2613192/1‑93 (MQ=255)

                    aCTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGACGGCGGCGCt                                                                     >  1:1501757/1‑93 (MQ=255)

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

                     cTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCtt                                                                    <  1:349185/93‑1 (MQ=255)

                     cTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCtt                                                                    <  2:3148128/93‑1 (MQ=255)

                      tGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCttt                                                                   >  2:846632/1‑93 (MQ=255)

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

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

                         gTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAGGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGAt                                                                >  1:334054/1‑93 (MQ=255)

                          tCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATa                                                               >  2:2354040/1‑93 (MQ=255)

                             cgcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAAcc                                                            <  1:1132704/93‑1 (MQ=255)

                              gcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCg                                                           >  1:1541997/1‑93 (MQ=255)

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

                               ccACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCgg                                                          >  2:2799428/1‑93 (MQ=255)

                               ccACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCgg                                                          >  1:2219209/1‑93 (MQ=255)

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

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

                                  ctgctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGgcgc                                                       <  2:4553789/93‑1 (MQ=255)

                                   tgctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCt                                                      >  1:1982128/1‑93 (MQ=255)

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

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

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

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

                                        ttttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGcc                                                 >  1:2029378/1‑93 (MQ=255)

                                         tttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCt                                                <  2:1037574/93‑1 (MQ=255)

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

                                            tccGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGc                                             <  2:1380786/92‑1 (MQ=255)

                                            aCCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGc                                             <  2:1339405/93‑1 (MQ=255)

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

                                                tATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGAt                                         >  1:3769003/1‑93 (MQ=255)

                                                tATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGAt                                         >  2:3890623/1‑93 (MQ=255)

                                                tATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGAt                                         <  2:1046979/93‑1 (MQ=255)

                                                tATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGGGCTTTGATAACCGGCGCTTCGCCTTGCAGAt                                         >  1:434369/1‑93 (MQ=255)

                                                 aTCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATa                                        <  1:225561/93‑1 (MQ=255)

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

                                                  tCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATaa                                       >  2:2928463/1‑93 (MQ=255)

                                                   cAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAAc                                      >  1:2718638/1‑93 (MQ=255)

                                                   cAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAAc                                      >  1:1879492/1‑93 (MQ=255)

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

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

                                                      tACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGAt                                   >  1:80794/1‑93 (MQ=255)

                                                       aCGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGAtt                                  >  1:3892336/1‑93 (MQ=255)

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

                                                       aCGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGAtt                                  >  1:3594889/1‑93 (MQ=255)

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

                                                         gcgcCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTaa                                <  2:2260914/93‑1 (MQ=255)

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

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

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

                                                           gcCACAATTTGCGCGCGACGGGAAGCCAGCGTTTCAGCATGGGGCGGCGGTTCTTTGATAACCGGCGCTGCGCCTCGCAGAGAACGATTaaca                              >  2:3913977/1‑91 (MQ=255)

                                                           gcCACAAAGTGCGCGTAACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAata                              <  2:1659644/93‑1 (MQ=255)

                                                            ccACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAatat                             >  2:403904/1‑93 (MQ=255)

                                                            ccACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAatat                             <  1:2968664/93‑1 (MQ=255)

                                                               caATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATcgc                          <  2:3740606/93‑1 (MQ=255)

                                                               caATTTGCGCGTCACGGGAAGCCAGCGGTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATcgc                          <  1:1868976/93‑1 (MQ=255)

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

                                                                  tttGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATCACGATTAATATCGCgca                       >  1:120067/1‑93 (MQ=255)

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

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

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

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

                                                                     gcgcGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCCGATCACGATTAATATCGCgcagca                    >  2:2702896/1‑93 (MQ=255)

                                                                     gcgcGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgcagca                    >  1:508305/1‑93 (MQ=255)

                                                                     gcgcGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgcagca                    >  1:2862997/1‑93 (MQ=255)

                                                                     gcgcGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgcagca                    >  2:3688397/1‑93 (MQ=255)

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

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

                                                                       gcgTGAAGGGAGCCCCGCGCTTCATCATCGGAAGACGGTCCTAAGATGACCGACGCTTCGCCTTGCAGAAAACGATTAATATCGCAAAGCAAc                  <  2:4209520/93‑1 (MQ=255)

                                                                        cgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACg                 >  2:2347934/1‑93 (MQ=255)

                                                                          tGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGcc               <  1:3894710/93‑1 (MQ=255)

                                                                          tGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGcc               <  1:4683225/93‑1 (MQ=255)

                                                                           gACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCt              <  2:4749724/93‑1 (MQ=255)

                                                                              gggAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCg           >  2:2146383/1‑93 (MQ=255)

                                                                              gggAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCg           >  1:592713/1‑93 (MQ=255)

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

                                                                                   gCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCt      >  2:3543136/1‑93 (MQ=255)

                                                                                   gCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCt      <  1:1963701/93‑1 (MQ=255)

                                                                                   gCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCt      <  1:257394/93‑1 (MQ=255)

                                                                                     cAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCttt    <  2:134857/93‑1 (MQ=255)

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

                                                                                       gcgcTTCATCATCGGGCGGGGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTcc  <  2:3114878/93‑1 (MQ=255)

                                                                                       gcgcTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTcc  >  2:1883956/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` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`