BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,563,585	(C)_5→6	coding (1320/1404 nt)	b2451 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,563,580	1	.	C	97.1%	643.5 / 10.0	204	coding (1325/1404 nt)	b2451	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (4/2); new base C (106/92); total (110/94)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 237 total aligned reads displayed.

GCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑CC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAA  >  NC_000913/2563488‑2563667

   aaTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑cc‑c                                                                                       <  1:1361700/93‑1 (MQ=255)

   aaTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑cc‑c                                                                                       >  1:5528322/1‑93 (MQ=255)

     tACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑cc‑cc‑c                                                                                    >  1:5124740/1‑93 (MQ=255)

      aCCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑cc‑cc‑cc                                                                                   <  1:2753573/93‑2 (MQ=255)

      aCCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑cc‑cc‑cc                                                                                   <  1:4726002/93‑2 (MQ=255)

       ccGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑cc‑cc‑ccg                                                                                  <  2:3020984/93‑3 (MQ=255)

       ccGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑cc‑cc‑ccg                                                                                  >  1:293655/1‑91 (MQ=255)

              cAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGAc                                                                            <  2:2471821/93‑1 (MQ=255)

              cAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGAc                                                                            >  1:3368856/1‑93 (MQ=255)

              cAGGAAGGAAATGCGAGCGATTTCACCGTCACCGGCACAACCGATCCGCCACAAAGAGGCGTACCAATGTCGATACAG‑T‑CCACC‑CGCGCggcc                                                                            >  2:4055067/1‑91 (MQ=255)

               aggaaggaAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACa                                                                           >  2:5063377/1‑93 (MQ=255)

               aggaaggaAATGCGAGAGATTTCACCGTCACCCGCACAACCGATCCGCCTAAAAGAGGCGTACCAATGGCGATATAT‑TGCC‑CC‑CCCGCGGACa                                                                           >  1:1055134/1‑93 (MQ=255)

                ggaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACaa                                                                          >  2:5515377/1‑93 (MQ=255)

                ggaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACaa                                                                          <  2:1141960/93‑1 (MQ=255)

                 gaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAAt                                                                         <  2:3800207/93‑1 (MQ=255)

                 gaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAAt                                                                         <  1:4020609/93‑1 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGGGCGGACATTCa                                                                       >  1:4076319/1‑93 (MQ=255)

                     gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACt                                                                     <  1:2979975/93‑1 (MQ=255)

                     gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACt                                                                     <  1:3543714/93‑1 (MQ=255)

                      aaaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACtt                                                                    <  2:4571007/93‑1 (MQ=255)

                       aaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTc                                                                   >  2:4198006/1‑93 (MQ=255)

                        aTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACCTCg                                                                  >  1:3462794/1‑93 (MQ=255)

                         tGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGGCCATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGt                                                                 >  2:4841514/1‑93 (MQ=255)

                          gCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTc                                                                <  1:3969082/93‑1 (MQ=255)

                          gCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTc                                                                >  2:3233602/1‑93 (MQ=255)

                            gAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCaa                                                              <  2:2393226/93‑1 (MQ=255)

                            gAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCaa                                                              <  1:5027342/93‑1 (MQ=255)

                            gAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCaa                                                              <  2:593678/93‑1 (MQ=255)

                             aGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAAt                                                             <  1:2457332/93‑1 (MQ=255)

                              gTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATg                                                            >  2:4610534/1‑93 (MQ=255)

                                gATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGAc                                                          <  1:2265899/93‑1 (MQ=255)

                                gATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGAc                                                          >  1:4308777/1‑93 (MQ=255)

                                gATTTCACCGTCACCGGCACAACCGAGCCGCCAAAAAGAGGCGTACCAATGTGGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGAc                                                          <  1:2176850/93‑1 (MQ=255)

                                 aTTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACt                                                         <  2:6124409/93‑1 (MQ=255)

                                     caccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCa                                                     >  1:87583/1‑93 (MQ=255)

                                      accgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCaa                                                    >  2:2340490/1‑93 (MQ=255)

                                      accgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCaa                                                    <  2:5853133/93‑1 (MQ=255)

                                       ccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAAc                                                   <  1:3082791/93‑1 (MQ=255)

                                       ccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAAc                                                   <  1:5841277/93‑1 (MQ=255)

                                        cgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACg                                                  <  2:3643334/93‑1 (MQ=255)

                                        cgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACg                                                  >  2:6497570/1‑93 (MQ=255)

                                         gtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACgg                                                 >  1:830015/1‑93 (MQ=255)

                                           caccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGa                                               >  2:1498544/1‑93 (MQ=255)

                                            gcctgCACAGCCGATCCGCCAAAAAGGGGCGTTCCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAGTCTCTTCGTCAATGACTGCCAACGGGAg                                              <  2:2116361/89‑1 (MQ=255)

                                            cccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAg                                              <  1:6504870/92‑1 (MQ=255)

                                             ccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgt                                             <  1:4615750/93‑1 (MQ=255)

                                               ggCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttg                                           <  2:5992826/93‑1 (MQ=255)

                                                gCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgt                                          <  1:2035277/93‑1 (MQ=255)

                                                gCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgt                                          <  1:4336220/93‑1 (MQ=255)

                                                gCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgt                                          >  2:1821408/1‑93 (MQ=255)

                                                gCACAACCGAACCGCCAAAAAGAGGCGTACTAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgt                                          <  1:5174600/93‑1 (MQ=255)

                                                 cacaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgtt                                         <  1:1828110/93‑1 (MQ=255)

                                                   caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                       >  2:2151764/1‑93 (MQ=255)

                                                   caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑CC‑CCACGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                       >  1:2111614/1‑93 (MQ=255)

                                                    aaCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTa                                      >  1:1199663/1‑93 (MQ=255)

                                                     aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTT‑CC‑CC‑CGCGCGGACAATCACTTCGTCAATGAATGCCAACGGGAGTTGTTGTAg                                     >  1:6683624/1‑93 (MQ=255)

                                                     aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                     >  2:5736629/1‑93 (MQ=255)

                                                     aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                     >  1:5299382/1‑93 (MQ=255)

                                                     aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                     >  1:5431188/1‑93 (MQ=255)

                                                      ccGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGc                                    >  2:192939/1‑93 (MQ=255)

                                                       cGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCt                                   >  2:2940990/1‑93 (MQ=255)

                                                        gATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                  >  2:5130256/1‑93 (MQ=255)

                                                        gATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                  >  1:737223/1‑93 (MQ=255)

                                                         aTCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgt                                 <  1:1226658/93‑1 (MQ=255)

                                                         aTCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgt                                 >  1:3054017/1‑93 (MQ=255)

                                                         aTCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgt                                 >  2:3964994/1‑93 (MQ=255)

                                                         aTCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgt                                 <  2:3368926/93‑1 (MQ=255)

                                                          tCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgtg                                <  2:3735823/93‑1 (MQ=255)

                                                          tCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgtg                                <  1:5744563/93‑1 (MQ=255)

                                                           ccGCCAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                               >  2:155437/1‑93 (MQ=255)

                                                             gcaaaaaaGAGGCGCACCAATGTCGATATAGCC‑CC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgc                             <  1:2550012/90‑1 (MQ=255)

                                                               cAAAAAGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgcgc                           >  1:3774227/1‑93 (MQ=255)

                                                                aaaaaGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGca                          >  2:5355384/1‑93 (MQ=255)

                                                                 aaaaGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaa                         >  2:28717/1‑93 (MQ=255)

                                                                 aaaaGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaa                         >  1:3539190/1‑93 (MQ=255)

                                                                   aaGAGGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaaca                       >  2:843173/1‑93 (MQ=255)

                                                                    agagGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaacaa                      <  1:1961171/93‑1 (MQ=255)

                                                                    agagGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaacaa                      <  2:1248697/93‑1 (MQ=255)

                                                                     gagGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCaacaac                     >  1:169996/1‑93 (MQ=255)

                                                                     gagGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCaacaac                     >  2:1761600/1‑93 (MQ=255)

                                                                     gagGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCaacaac                     >  1:4278014/1‑93 (MQ=255)

                                                                     gagGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGGTGTTGTAGCTGTGGGCGCaacaac                     >  1:2637436/1‑93 (MQ=255)

                                                                      agGCGTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAacaaca                    <  2:5954737/93‑1 (MQ=255)

                                                                          gTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGcc                >  1:6400540/1‑93 (MQ=255)

                                                                          gTACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGcc                <  2:3053258/93‑1 (MQ=255)

                                                                           tACCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGccc               >  1:5755070/1‑93 (MQ=255)

                                                                            aCCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa              >  1:6241615/1‑93 (MQ=255)

                                                                            aCCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa              <  2:2549029/93‑1 (MQ=255)

                                                                            aCCAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa              >  2:5625435/1‑93 (MQ=255)

                                                                             ccAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCag             >  2:6730692/1‑93 (MQ=255)

                                                                              cAATGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCaga            >  2:2134877/1‑93 (MQ=255)

                                                                               aaTGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCagag           >  2:4277234/1‑93 (MQ=255)

                                                                                aTGTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGc          <  2:5903482/93‑1 (MQ=255)

                                                                                  gTCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGAtt        >  2:1459482/1‑93 (MQ=255)

                                                                                   tCGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCttt       >  1:3864521/1‑93 (MQ=255)

                                                                                    cGATATAGCC‑CC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTa      <  1:6208127/93‑1 (MQ=255)

                                                                                    cGATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTa      >  1:1426931/1‑93 (MQ=255)

                                                                                     gATATAG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAc     <  2:2648802/93‑1 (MQ=255)

                                                                                       tataG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCa   <  1:1900580/93‑1 (MQ=255)

                                                                                        ataG‑TCCC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCaa  >  2:2729851/1‑93 (MQ=255)

GCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAG‑T‑CC‑CC‑CGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAA  >  NC_000913/2563488‑2563667

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