BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,850,430	C→T	5.9%	T717T (ACC→ACT)	PP_3397 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,850,430	0	C	T	5.9%	302.0 / 6.6	103	T717T (ACC→ACT)	PP_3397	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (50/46); new base T (3/3); total (54/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 104 total aligned reads displayed.

GACAGCCACGGTCACGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATGAAACGGCTGTACTC  >  NC_002947/3850283‑3850579

gACAGCCACGGTCACGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCgg                                                                                                                                                     >  2:2175375/1‑150 (MQ=255)

 aCAGCCACGGTCACGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGc                                                                                                                                                    >  1:1595432/1‑150 (MQ=255)

     ccACGGTCACGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCt                                                                                                                                                >  2:125698/1‑150 (MQ=255)

     ccACGGTCACGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCt                                                                                                                                                >  1:1570564/1‑150 (MQ=255)

      cACGGTCACGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCCATAGCGTAACCGTGCGGGTAAACACTACCACGTGCGTGCAGACCTTGGGGACGGCCACGGTAGCGCTGACGGGCCGCTg                                                                                                                                               >  1:1142945/1‑150 (MQ=255)

       actgggaCGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGCCGGCAACGGTAGCGCTGACCGGCCGCTgg                                                                                                                                              >  1:695852/7‑150 (MQ=255)

          gTCACGGCGCGCTCCAACCACCCTTACAACTCGCACATCAGTCGCACCTCCTCCGTGCCCACCGGCAATACCGTAACCGTTCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGgcgc                                                                                                                                           <  1:2238433/150‑1 (MQ=255)

            cACGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCt                                                                                                                                         <  2:1073451/150‑1 (MQ=255)

             aCGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCtg                                                                                                                                        <  2:1154846/150‑1 (MQ=255)

             aCGGAGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCtg                                                                                                                                        >  1:1418360/1‑150 (MQ=255)

              cGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCtgg                                                                                                                                       <  1:1387145/150‑1 (MQ=255)

               ggCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCtggt                                                                                                                                      >  1:1083694/1‑150 (MQ=255)

                        caacaaCCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCaa                                                                                                                             >  1:966120/1‑150 (MQ=255)

                        caacaaCCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCaa                                                                                                                             <  1:2175375/150‑1 (MQ=255)

                           caacCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCAACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAg                                                                                                                          >  1:160982/1‑150 (MQ=255)

                              ccGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcac                                                                                                                       >  2:1646325/1‑150 (MQ=255)

                                gTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcacca                                                                                                                     <  2:1418360/150‑1 (MQ=255)

                                gTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcacca                                                                                                                     <  2:1984996/150‑1 (MQ=255)

                                gTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcacca                                                                                                                     >  2:655814/1‑150 (MQ=255)

                                 ttACAACTGGTACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcacca                                                                                                                     <  1:486522/149‑1 (MQ=255)

                                 ttACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCg                                                                                                                                                      >  2:794927/1‑116 (MQ=255)

                                 ttACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCg                                                                                                                                                      <  1:794927/116‑1 (MQ=255)

                                 ttACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcacca                                                                                                                     >  2:486522/1‑149 (MQ=255)

                                   aCAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCCTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGa                                                                                                                  <  1:2294613/150‑1 (MQ=255)

                                        tGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTCACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTGGCGCGGACCGGCCGCTGGCGCCTGGTGGTAAGCAATCGCACCACGctcatc                                                                                                             >  1:148936/1‑150 (MQ=255)

                                        tGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGatcatc                                                                                                             >  2:234718/1‑150 (MQ=255)

                                        tGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGatc                                                                                                                >  1:1824762/1‑147 (MQ=255)

                                        tGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGatc                                                                                                                <  2:1824762/147‑1 (MQ=255)

                                              aTCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACt                                                                                                       >  2:625702/1‑150 (MQ=255)

                                               tCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTg                                                                                                      >  2:225126/1‑150 (MQ=255)

                                                 aGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCg                                                                                                    <  2:1142945/150‑1 (MQ=255)

                                                  gTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACCACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCgg                                                                                                   >  2:912668/1‑150 (MQ=255)

                                                   tGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcacc                                                                                                                      >  1:620548/1‑130 (MQ=255)

                                                   tGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcacc                                                                                                                      <  2:620548/130‑1 (MQ=255)

                                                     gCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCGATAGCACCACGATCATCCCGACTGCGGCAc                                                                                                >  2:2350256/1‑150 (MQ=255)

                                                        ccTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGACagggc                                                                                             >  1:1242403/1‑146 (MQ=255)

                                                         cTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCGATAGCACCACGATCATCCCGACTGCGGCACCGAc                                                                                            >  1:1733876/1‑150 (MQ=255)

                                                           cgtcgGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACgg                                                                                          <  1:109904/150‑1 (MQ=255)

                                                             tcgGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGcg                                                                                        >  1:8094/1‑150 (MQ=255)

                                                              cgGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGcga                                                                                       <  2:1083694/150‑1 (MQ=255)

                                                              cgGTGACCACCGGAAATACCGTAAACGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGcga                                                                                       <  2:1229745/150‑1 (MQ=255)

                                                               ggTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACTGGCCGCTGGCGCCTGGTGGTaa                                                                                                                                 >  2:627884/1‑107 (MQ=255)

                                                               ggTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACTGGCCGCTGGCGCCTGGTGGTaa                                                                                                                                 >  2:628413/1‑107 (MQ=255)

                                                               ggTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACTGGCCGCTGGCGCCTGGTGGTaa                                                                                                                                 <  1:627884/107‑1 (MQ=255)

                                                               ggTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACTGGCCGCTGGCGCCTGGTGGTaa                                                                                                                                 <  1:628413/107‑1 (MQ=255)

                                                               ggTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGcgac                                                                                      >  1:1301983/1‑150 (MQ=255)

                                                               ggTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGcgac                                                                                      <  1:1646325/150‑1 (MQ=255)

                                                                gTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACTGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACg                                                                                                                   >  1:464450/1‑120 (MQ=255)

                                                                gTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACTGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACg                                                                                                                   <  2:464450/120‑1 (MQ=255)

                                                                   accaccGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATc                                                                                  <  1:125698/150‑1 (MQ=255)

                                                                     caccGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATcac                                                                                >  1:75801/1‑150 (MQ=255)

                                                                        cGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCtggt                                                                                                                                      <  1:1491345/93‑1 (MQ=255)

                                                                        cGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCtggt                                                                                                                                      >  2:1491345/1‑93 (MQ=255)

                                                                          gCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCAccacca                                                                           <  1:704176/150‑1 (MQ=255)

                                                                            aaTACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCg                                                                                                    <  2:846462/123‑1 (MQ=255)

                                                                            aaTACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCg                                                                                                    >  1:846462/1‑123 (MQ=255)

                                                                               aCCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCaa                                                                      >  2:1365075/1‑150 (MQ=255)

                                                                               aCCGTAACAGTGCGGGTACTCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCaa                                                                      <  2:1157259/150‑1 (MQ=255)

                                                                                  gTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCAACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGc                                                                   >  2:885951/1‑150 (MQ=255)

                                                                                      ccGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTaa                                                                                                                                 <  1:2391799/84‑1 (MQ=255)

                                                                                      ccGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTaa                                                                                                                                 >  2:2391799/1‑84 (MQ=255)

                                                                                            gggTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCa                                                         >  1:9570/1‑150 (MQ=255)

                                                                                             ggTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGCGCTTCAGCACCAGCCAAGGCACCACGCGCAc                                                        <  1:1872340/150‑1 (MQ=255)

                                                                                              gtgaCCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCAcc                                                       >  1:2404642/4‑150 (MQ=255)

                                                                                              gTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCAcgc                                                            >  2:1502448/1‑145 (MQ=255)

                                                                                              gTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCAcgc                                                            <  1:1502448/145‑1 (MQ=255)

                                                                                                  ccACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcacca                                                                                                                     <  2:2042205/84‑1 (MQ=255)

                                                                                                  ccACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGcacca                                                                                                                     >  1:2042205/1‑84 (MQ=255)

                                                                                                  ccACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCa                                                   >  1:1079924/1‑150 (MQ=255)

                                                                                                   cACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCaa                                                  >  2:199718/1‑150 (MQ=255)

                                                                                                    aCTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCGAGGCACCACGCGCACCGTCAAc                                                 <  1:625702/150‑1 (MQ=255)

                                                                                                    aCTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAAc                                                 <  1:1651470/150‑1 (MQ=255)

                                                                                                    aCTACCACGGGCGTCAAGACCTTGGGGACGGCTACGGTAGCTCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAAc                                                 <  2:193741/150‑1 (MQ=255)

                                                                                                      tACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGt                                               >  1:345521/1‑150 (MQ=255)

                                                                                                           cGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCg                                          >  2:309106/1‑150 (MQ=255)

                                                                                                            gggCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGt                                         >  1:1534611/1‑150 (MQ=255)

                                                                                                               cGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCCTCTCTACTTCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCa                                      <  2:580672/150‑1 (MQ=255)

                                                                                                               cGTCCAGACCTTGGGCACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGGGGCCAAAACGGGCATCAGCAGGATCGAGCCTGCGGCACCGACGGCGACGATCTCCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCa                                      <  1:2350256/150‑1 (MQ=255)

                                                                                                                 tCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGCCTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACTACGCGCACCGTCAACGTAGTCGTGCAGt                                    <  2:1151608/150‑1 (MQ=255)

                                                                                                                    aGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGGGCCTGGGGGGAAGCAATAGCACCACGATCATCGCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTaaa                                 >  2:757189/1‑150 (MQ=255)

                                                                                                                    aGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTaaa                                 >  1:2263632/1‑150 (MQ=255)

                                                                                                                      aCCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAg                               >  2:1506281/1‑150 (MQ=255)

                                                                                                                      aCCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAg                               >  1:2175323/1‑150 (MQ=255)

                                                                                                                       ccTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGc                              <  2:2263632/150‑1 (MQ=255)

                                                                                                                           ggggACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGa                          <  2:1029908/150‑1 (MQ=255)

                                                                                                                           ggggACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGa                          <  1:234718/150‑1 (MQ=255)

                                                                                                                                 ggCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTc                    <  2:160982/150‑1 (MQ=255)

                                                                                                                                 ggCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCAACCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAATCGTGCGGTAAAAGGCGGGGGCTTc                    >  2:418704/1‑150 (MQ=255)

                                                                                                                                 ggCCACGGTAGCGCGGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCTCTCCTAGCCAACCCAACACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTc                    <  2:695852/150‑1 (MQ=255)

                                                                                                                                    cACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCaa                                                                      >  1:774451/1‑97 (MQ=255)

                                                                                                                                    cACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCaa                                                                      <  2:774451/97‑1 (MQ=255)

                                                                                                                                     aCGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATg                <  2:2175323/150‑1 (MQ=255)

                                                                                                                                            cgcTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATGAAACGGc         <  1:1057500/150‑1 (MQ=255)

                                                                                                                                             gcTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATGAAACGGCt        <  2:1595432/150‑1 (MQ=255)

                                                                                                                                             gcTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATGAAACGGCt        <  1:225126/150‑1 (MQ=255)

                                                                                                                                             gcTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCCGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATGAAACGGCt        <  1:2148191/150‑1 (MQ=255)

                                                                                                                                               tGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATGaaa             >  2:1850324/1‑143 (MQ=255)

                                                                                                                                               tGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATGaaa             <  1:1850324/143‑1 (MQ=255)

                                                                                                                                               tGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATGAAACGGCTGt      >  2:608400/1‑150 (MQ=255)

                                                                                                                                                   cGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCGTCGATGAAACGGCTGTACTc  >  1:125508/1‑150 (MQ=255)

GACAGCCACGGTCACGGCGCGCTCCAACAACCGTTACAACTGGGACATCAGTGGCACCTCGTCGGTGACCACCGGCAATACCGTAACCGTGCGGGTAACCACTACCACGGGCGTCCAGACCTTGGGGACGGCCACGGTAGCGCTGACCGGCCGCTGGCGCCTGGTGGTAAGCAATAGCACCACGATCATCCCGACTGCGGCACCGACGGCGACGATCACCACCAGCCAAGGCACCACGCGCACCGTCAACGTAGTCGTGCAGTAAAAGCCGGAGGCTTCGATGAAACGGCTGTACTC  >  NC_002947/3850283‑3850579

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