BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,686,486	G→T	5.3%	G284V (GGC→GTC)	PP_3251 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,686,486	0	G	T	5.3%	322.5 / 2.7	102	G284V (GGC→GTC)	PP_3251	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (38/58); new base T (2/4); total (40/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

CGGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGTC  >  NC_002947/3686337‑3686627

cGGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCCTGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTCCCGCCTGAATGACCGGGGTATCGCCCCGACAGCCTTCAC‑‑CCGCCAGTTCGACCTGCAAgg                                                                                                                                               <  2:1553132/150‑1 (MQ=255)

cGGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGTGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAgg                                                                                                                                               <  1:228759/150‑1 (MQ=255)

cGGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAgg                                                                                                                                               <  1:1690388/150‑1 (MQ=255)

cGGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAgg                                                                                                                                               >  2:442573/1‑150 (MQ=255)

 ggCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGc                                                                                                                                              >  1:380998/1‑150 (MQ=255)

    cACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGaa                                                                                                                                           <  1:1116635/150‑1 (MQ=255)

     aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGCCTTGCAGGTCGAAc                                                                                                                                          <  2:853859/150‑1 (MQ=255)

     aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAAc                                                                                                                                          >  1:210239/1‑150 (MQ=255)

      cgcccgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGCCTTGCAGGTCGAACt                                                                                                                                         <  2:855232/150‑1 (MQ=255)

         ccgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGCCTTGCAGGTCGAACTggc                                                                                                                                      <  2:1773035/150‑3 (MQ=255)

         ccgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGtt                                                                                                                                      >  2:331867/1‑150 (MQ=255)

          cgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGttt                                                                                                                                     <  2:1383155/150‑1 (MQ=255)

          cgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGttt                                                                                                                                     <  2:1618720/150‑1 (MQ=255)

           gccTGGCCATTCAGTCGATGAACCAGCGAAGCTCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGTCCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTg                                                                                                                                    <  2:109794/150‑1 (MQ=255)

           gccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTg                                                                                                                                    >  2:431849/1‑150 (MQ=255)

             cTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCa                                                                                                                                  <  1:1341278/150‑1 (MQ=255)

                  aaTTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCt                                                                                                                             >  2:1426584/1‑150 (MQ=255)

                    ttCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTg                                                                                                                           <  1:1175628/150‑1 (MQ=255)

                       aGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGGTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAAAGAAAGCAGTAACGTTA‑TCCACA‑ATGATGGGCGTCAGTTCGAACTGCCAGGCGAACTGTTTGCAGCCCTTGAAc                                                                                                                        >  1:223858/1‑150 (MQ=255)

                        ggCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACg                                                                                                                       >  2:936962/1‑150 (MQ=255)

                         gcgcTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTCGTACGCCAGGGCGCGCTTAGGCTACCACCTGAATGACCGGGGCATCGCCACGCAAGACCTCAT‑‑CCGCCAGTTCGACCTGCTAGGCGAACTGTTTGCAGCCCTTGAACGc                                                                                                                      >  2:737268/1‑150 (MQ=255)

                           gcTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCaa                                                                                                                    <  1:1143657/150‑1 (MQ=255)

                               aaCCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCgg                                                                                                                <  2:544901/150‑1 (MQ=255)

                                 ccAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCaa                                                                                                                    <  1:1111028/144‑1 (MQ=255)

                                 ccAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCaa                                                                                                                    >  2:1111028/1‑144 (MQ=255)

                                 ccAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGt                                                                                                              <  1:852403/150‑1 (MQ=255)

                                 ccAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGt                                                                                                              >  1:1212543/1‑150 (MQ=255)

                                 ccAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGt                                                                                                              <  1:852505/150‑1 (MQ=255)

                                 ccAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGt                                                                                                              >  2:338499/1‑150 (MQ=255)

                                   aGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGTCACGCCACACCTCAT‑‑CCGCCCGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCt                                                                                                            <  2:1055552/150‑1 (MQ=255)

                                     cTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGTCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTAt                                                                                                          >  1:17625/1‑150 (MQ=255)

                                       aaGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGAGCAAACCCCTTGAACGCAAGCGGGTCTATgg                                                                                                        <  1:528272/150‑1 (MQ=255)

                                        aGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGc                                                                                                       <  2:1924647/150‑1 (MQ=255)

                                            ccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGCCTTGCAGGTCGAACTggcg                                                                                                                                     <  1:564181/116‑4 (MQ=255)

                                            ccGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGCCTTGCAGGTCGAACTggcg                                                                                                                                     >  2:564181/1‑113 (MQ=255)

                                             cgcgCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAAtt                                                                                                  <  1:1365289/150‑1 (MQ=255)

                                               cgcATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGa                                                                                                >  1:1639296/1‑150 (MQ=255)

                                                gcATGATCCGCAGGCGTTCTACTCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                               <  2:954885/150‑1 (MQ=255)

                                                gcATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                               <  1:296154/150‑1 (MQ=255)

                                                gcATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                               <  2:1772916/150‑1 (MQ=255)

                                                gcATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                                               <  2:281298/150‑1 (MQ=255)

                                                 cATGATCCGCAGGCGTTCTCCGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGatg                                                                                              <  1:659226/150‑1 (MQ=255)

                                                   tGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGatgca                                                                                            >  1:177093/1‑150 (MQ=255)

                                                    gATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGatgcat                                                                                           <  2:179783/150‑1 (MQ=255)

                                                    gATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGatgcat                                                                                           <  1:869421/150‑1 (MQ=255)

                                                       ccGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGcc                                                                                        <  2:1212543/150‑1 (MQ=255)

                                                                ttCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTgc                                                                               >  2:1484861/1‑150 (MQ=255)

                                                                   tACGCCAGGGCGCGCTTAGGCTACCGCCTGTATGACCGGGGTATCGGCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGACCGCAAGCGGGTCTATTGCCAATTGATGCATGCCGACCATTgccgc                                                                            <  2:413779/150‑1 (MQ=255)

                                                                    aCGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTgccgcc                                                                           <  1:165460/150‑1 (MQ=255)

                                                                       ccAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTgc                                                                        >  1:269724/1‑150 (MQ=255)

                                                                        cAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTgca                                                                       <  2:975130/150‑1 (MQ=255)

                                                                            gcgcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCt                                                                                                                             <  1:154308/92‑1 (MQ=255)

                                                                            gcgcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCt                                                                                                                             >  2:154308/1‑92 (MQ=255)

                                                                            gcgcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCTCCTGcagcag                                                                   <  2:1261758/150‑1 (MQ=255)

                                                                             cgcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                                                  <  1:1428875/150‑1 (MQ=255)

                                                                                 ctttgagtACCGCCGGAATGATCGGGGTAGCGCCACGCCAGACCTCAT‑‑CCGCCCGTTCGACCTGCGAGGCGAACTGTTTGCAGCCCTTGACCGCAAGCGGGTCTATGGCCAATTGGTGCATGCCGACCATTGCCGCCTGCAGCAGCGCAc                                                              <  2:106715/143‑1 (MQ=255)

                                                                                 cTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTAGCGGCTGTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCAc                                                              >  1:118490/1‑150 (MQ=255)

                                                                                   tAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGa                                                            <  2:1525404/150‑1 (MQ=255)

                                                                                      gCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGatca                                                         >  2:725937/1‑150 (MQ=255)

                                                                                        tACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCACTTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGatcaat                                                       >  1:1373641/1‑150 (MQ=255)

                                                                                        tACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGatcaat                                                       <  1:936962/150‑1 (MQ=255)

                                                                                        tACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGatcaat                                                       >  1:1123960/1‑150 (MQ=255)

                                                                                        tACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGatcaat                                                       <  1:1840092/150‑1 (MQ=255)

                                                                                          ccGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTAt                                                                                                          >  1:925795/1‑97 (MQ=255)

                                                                                          ccGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTAt                                                                                                          <  2:925795/97‑1 (MQ=255)

                                                                                            gCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                   <  1:1161874/150‑1 (MQ=255)

                                                                                            gCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                   <  1:1161915/150‑1 (MQ=255)

                                                                                            gCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                   <  1:1836987/150‑1 (MQ=255)

                                                                                            gCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAAGTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATcaatcaat                                                   >  1:1057382/1‑150 (MQ=255)

                                                                                             ccTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtg                                                  >  2:147133/1‑150 (MQ=255)

                                                                                                 aaTGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGatgca                                                                                            >  2:1760146/1‑104 (MQ=255)

                                                                                                 aaTGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGatgca                                                                                            <  1:1760146/104‑1 (MQ=255)

                                                                                                  aTGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtggctgg                                             <  1:442573/150‑1 (MQ=255)

                                                                                                  aTGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtggctgg                                             <  1:319602/150‑1 (MQ=255)

                                                                                                           gTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGcc                                           >  2:1902769/1‑143 (MQ=255)

                                                                                                           gTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGcc                                           <  1:1902769/143‑1 (MQ=255)

                                                                                                                gccacgccaGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGcc                               <  2:1639296/150‑1 (MQ=255)

                                                                                                                gccacgccaGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGcc                               <  1:1233043/150‑1 (MQ=255)

                                                                                                                 ccacgccaGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCa                              <  1:1426584/150‑1 (MQ=255)

                                                                                                                    cgccaGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGAACAACGGCCCCTGCAGCAGCACGCGGTCAATCAATGGCTGGCCGAACCCGGCGCCATGc                           <  1:331867/150‑1 (MQ=255)

                                                                                                                      ccaGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGccc                         >  1:1620122/1‑150 (MQ=255)

                                                                                                                        aGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGgc                                  <  2:1698902/139‑1 (MQ=255)

                                                                                                                        aGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGgc                                  >  1:1698902/1‑139 (MQ=255)

                                                                                                                        aGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGccc                         <  2:1897646/148‑1 (MQ=255)

                                                                                                                        aGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGccc                         >  1:1897646/1‑148 (MQ=255)

                                                                                                                           ccTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGGCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCGCGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCCtt                    >  1:17965/1‑150 (MQ=255)

                                                                                                                             tCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCGATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCg                  >  1:1535643/1‑150 (MQ=255)

                                                                                                                             tCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCg                  >  1:1822836/1‑150 (MQ=255)

                                                                                                                                t‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTg               <  1:737268/150‑1 (MQ=255)

                                                                                                                                  cccGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTg               <  2:1602047/149‑1 (MQ=255)

                                                                                                                                     gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATTCCGACCATTgc                                                                               >  1:1076284/1‑83 (MQ=255)

                                                                                                                                     gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTgc                                                                               <  2:1076284/83‑1 (MQ=255)

                                                                                                                                     gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCg            >  2:538008/1‑150 (MQ=255)

                                                                                                                                     gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCg            <  1:1579513/150‑1 (MQ=255)

                                                                                                                                     gCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCg            >  1:956349/1‑150 (MQ=255)

                                                                                                                                         ggtCGACCTGCAAGGCGAACTGGATGCAGGCCTTGAACTCAAGCCGGCCTATGGCCAATTGATGCATGCCGACCATTCCCGCCTGCAGCAGCTCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCAtt        <  1:510841/148‑1 (MQ=255)

                                                                                                                                           tCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTAc      <  2:181426/150‑1 (MQ=255)

                                                                                                                                               ccTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGtc  <  1:538008/150‑1 (MQ=255)

                                                                                                                                                 tGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGACGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGcc                                           <  2:1563972/107‑1 (MQ=255)

                                                                                                                                                 tGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGACGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGcc                                           >  1:1563972/1‑107 (MQ=255)

CGGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCAT‑‑CCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCCGCATTCGTTGCCGCATTACAGTC  >  NC_002947/3686337‑3686627

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