BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,297,315	T→G	11.6%	L536R (CTG→CGG)	PP_2893 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,297,315	0	T	G	11.6%	668.8 / 25.3	267	L536R (CTG→CGG)	PP_2893	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (70/149); new base G (32/4); total (114/153)
Fisher's exact test for biased strand distribution p-value = 8.17e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 278 total aligned reads displayed.

GCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGACA  >  NC_002947/3297166‑3297464

  tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCtgt                                                                                                                                                      <  1:5563660/150‑1 (MQ=255)

  tCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCggt                                                                                                                                                      >  1:5730164/1‑150 (MQ=255)

   cAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCtgtg                                                                                                                                                     <  2:3513993/150‑1 (MQ=255)

    caGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTgg                                                                                                                                                    <  2:4077759/149‑1 (MQ=255)

     aGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGc                                                                                                                                                   <  1:3557648/150‑1 (MQ=255)

       gCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACCGCCACCTGTGGccc                                                                                                                                                 >  1:4898110/1‑149 (MQ=255)

         cGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCAGCTGTGGCCGAg                                                                                                                                               >  2:5390128/1‑150 (MQ=255)

         cGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACGGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACAGCTGGCCGAg                                                                                                                                               <  2:6344754/150‑1 (MQ=255)

          gagaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAgg                                                                                                                                              <  1:2325157/150‑1 (MQ=255)

            gaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCGGTGGCCGAGGCt                                                                                                                                            >  1:939014/1‑150 (MQ=255)

            gaCCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCt                                                                                                                                            <  1:1721263/150‑1 (MQ=255)

               ccAGGGGGCAGGGCACAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCGGGGGCCGCGGCGAcc                                                                                                                                         >  2:1132300/1‑149 (MQ=255)

               ccAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACa                                                                                                                                         >  2:4738581/1‑150 (MQ=255)

               ccAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACa                                                                                                                                         >  2:4541688/1‑150 (MQ=255)

                   ggggCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAATGc                                                                                                                                     >  2:4756808/1‑150 (MQ=255)

                     ggcagggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCg                                                                                                                                   >  2:5697348/1‑150 (MQ=255)

                       cagggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCggg                                                                                                                                 <  2:282607/150‑1 (MQ=255)

                          ggcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCa                                                                                                                              >  2:2213447/1‑150 (MQ=255)

                           gcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCaa                                                                                                                             <  1:4342882/150‑1 (MQ=255)

                           gcaAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGAATTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCaa                                                                                                                             <  1:5873716/150‑1 (MQ=255)

                            caAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAg                                                                                                                            <  2:3941316/150‑1 (MQ=255)

                            caAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAg                                                                                                                            <  2:376130/150‑1 (MQ=255)

                             aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAgg                                                                                                                           <  2:428128/150‑1 (MQ=255)

                             aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAgg                                                                                                                           <  1:2784723/150‑1 (MQ=255)

                             aaaaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAgg                                                                                                                           <  2:1239255/150‑1 (MQ=255)

                               aaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGc                                                                                                                         <  1:5257624/150‑1 (MQ=255)

                               aaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCAGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGc                                                                                                                         >  2:6307600/1‑150 (MQ=255)

                               aaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCGGTGGGCGCGGCGGCCACGCCGGGTCAAGGGc                                                                                                                         >  1:1785268/1‑150 (MQ=255)

                               aaCTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCGGTGGCCGAGGCGACAACGCCGGGTCAAGGGc                                                                                                                         >  1:3613694/1‑150 (MQ=255)

                                    gACCCTGGGCACACTGGCCATCCAGCTGGGT‑ACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTg                                                                                                                   <  1:1553191/150‑1 (MQ=255)

                                      cccTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTgg                                                                                                                  <  2:3651056/150‑1 (MQ=255)

                                       ccTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGa                                                                                                                 <  1:2166831/150‑1 (MQ=255)

                                        cTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCCGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGAt                                                                                                                >  2:3717527/1‑150 (MQ=255)

                                            gCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCg                                                                                                            <  1:2466039/150‑1 (MQ=255)

                                               cacTGGCCATCCAGCTGGGTGACACCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCg                                                                                                         <  1:696724/150‑1 (MQ=255)

                                                acTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCgg                                                                                                        <  1:1176293/150‑1 (MQ=255)

                                                    gCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCGGTGGCCGAGGCGACAACGCCGGGTCAACGGCCGGCGGGATGTCGCCCGCAcc                                                                                                    >  1:2025754/1‑150 (MQ=255)

                                                      cATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCtt                                                                                                  <  2:2634889/150‑1 (MQ=255)

                                                      cATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCtt                                                                                                  <  2:5200287/150‑1 (MQ=255)

                                                      cATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACCTGTGGCCGAGGCGACAATGCCGGGTCAAGGGCCCGCTGGATTTCGCCCGccccat                                                                                                  >  1:1023661/1‑148 (MQ=255)

                                                       aTCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTg                                                                                                 <  2:3005124/150‑1 (MQ=255)

                                                       aTCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTg                                                                                                 <  2:1312590/150‑1 (MQ=255)

                                                       aTCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTg                                                                                                 >  2:487177/1‑150 (MQ=255)

                                                        tCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGc                                                                                                <  2:4619828/150‑1 (MQ=255)

                                                        tCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGc                                                                                                <  2:5128462/150‑1 (MQ=255)

                                                              tGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAACGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGccc                                                                                          >  2:2664249/1‑150 (MQ=255)

                                                                 gTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAATGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCATGCAGGCCCCGc                                                                                       >  2:854224/1‑150 (MQ=255)

                                                                      aaCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGGCTACGGCAGCGCTGTGCCAGGCTACCAGGCCGGGTGGATGGCCGGCTGGGTGTCTCCGGCCCCGTGCGTGCCCCGccacgt                                                                                  >  1:6120599/1‑145 (MQ=255)

                                                                      aaCCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCACCGGTGGCCGAGGCTACAAGGCCCGGTCAAGGGCCGGCTGGATGTCGCCGGCACCATGCAGGCCCCCCAAGGt                                                                                  >  2:2483435/1‑150 (MQ=255)

                                                                        ccGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCTGTGGCCGCGGCGACAAGGCCGGGTCAAGGGCCGGCTGGCTGTCGCCGGCACCGTGCAGGCCCCCCCAGGTAc                                                                                >  2:3203742/1‑150 (MQ=255)

                                                                         cGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCGGTGGCCGAGGCGACAAGGCCGGGTCAAGGGCCGGCGGGATGTCGCCGGCACCGTGCAGGCCCCGCAAGGTACg                                                                               >  2:5323982/1‑150 (MQ=255)

                                                                         cGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCAGCTGTGGCCGAGGCTACAACGCCCGGTCAAGGGCCGGCTGGATTTCGCCGGCACCGTGCAGGCCCCCCAAGGTACg                                                                               >  1:6703626/1‑150 (MQ=255)

                                                                           cATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCt                                                                             >  2:3372658/1‑150 (MQ=255)

                                                                              cAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCAc                                                                          <  1:4628378/150‑1 (MQ=255)

                                                                               aaCGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCAcc                                                                         <  1:3064412/150‑1 (MQ=255)

                                                                               aaCGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCACGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCCGCTGGATGTCGCCGGCACCATGCAGGCCCCCCCAGGTACGCTCAcc                                                                         >  2:2041560/1‑150 (MQ=255)

                                                                                       tGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCATGCAGGCCCCGCAAGGTACGCTCACCCAGCAggg                                                                 >  1:3992255/1‑150 (MQ=255)

                                                                                         gCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGcc                                                               >  2:4448503/1‑150 (MQ=255)

                                                                                         gCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCAGTTGCCGAGGCGACAAGGCCGGGTCAAGGGCCGGCGGGATGTCGCCGGCACCATGCAGGCCCCCCAAGGTACGCTCACCCAGCAGGGcc                                                               >  1:2107733/1‑150 (MQ=255)

                                                                                          cAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCTGCAGGGTGAGCGTACCTTGCGGGGcct                                                              <  1:6230123/150‑2 (MQ=255)

                                                                                            gCTTGCAACAGCGCCTGGCTGGGCGACTCGCCCTTGCCGGGCCGCGCCTCGGCCAGCGGGGGCCGAGGCTACCCGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑Agc                                                           <  1:6584310/150‑1 (MQ=255)

                                                                                              ttGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCTGGCCCTGCAGGGTGAGCGTACCTTGCGGGGCC‑Tgcaa                                                         <  1:4977491/150‑3 (MQ=37)

                                                                                               tGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCTGGCCCTGCAGGGTGAGCGTACCTTGCGGGGCC‑Tgcaag                                                        <  2:6805535/150‑4 (MQ=37)

                                                                                                 cAACAGCGCCTGGCTGGGCGACTCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTGGCCCTGCAGGGTG‑Agcgtacc                                                      <  1:1972147/150‑5 (MQ=255)

                                                                                                 cAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTg                                                      <  1:2628176/150‑1 (MQ=255)

                                                                                                 cAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTg                                                      <  2:6196860/150‑1 (MQ=255)

                                                                                                  aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTgg                                                     <  1:5074348/150‑1 (MQ=255)

                                                                                                  aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCGGGCTGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑CGCGCTTgg                                                     >  2:4992014/1‑150 (MQ=255)

                                                                                                  aaCAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCGCGGCCCCCAGTGGCCGAGGCTACAAGGCCGGGGCAAGGGCCGGCTGGATTTAGCCCGCACCATGCAGGCCCCCCACGGTACGCACACCCGGCAGGGCC‑CCCGCtttg                                                     >  1:4410006/1‑148 (MQ=255)

                                                                                                      gcgcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCCAAGCGC‑TGGccct                                                 <  1:3467132/150‑2 (MQ=255)

                                                                                                      gcgcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCa                                                 <  1:469122/150‑1 (MQ=255)

                                                                                                       cgcCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAg                                                <  1:6500390/150‑1 (MQ=255)

                                                                                                         cccggctggGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑CGCGCTTGGCCCAGGc                                              <  2:1816764/147‑1 (MQ=255)

                                                                                                         ccTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCTGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGc                                              <  2:6622995/150‑1 (MQ=255)

                                                                                                             gctggGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAATGCCGGGTCAACGGCCGGCTGGGTTTCGCCGGCCCCTTGCAGGCCCCGCAAGGTACGCTCACCCGGCAGGGCC‑AGCGCGTGGCCCAGGCCGaa                                          >  2:152083/1‑150 (MQ=255)

                                                                                                                gggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAAc                                       <  1:3759235/150‑1 (MQ=255)

                                                                                                                 ggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCCAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAAcc                                      <  2:2412233/150‑1 (MQ=255)

                                                                                                                 ggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAAcc                                      <  2:3433241/150‑1 (MQ=255)

                                                                                                                 ggCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCGTGCAGGCCCCGCAAGGTACGCTCACCCAGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAAcc                                      >  2:6646797/1‑150 (MQ=255)

                                                                                                                   cGAATCGACCTTGACTTGCCGCGCCGCGGCCCGCTGTGGCCGAGGCTACAAGGCCCGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGTAGGCCCCGCAAGGTACGCTCACCCAGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGc                                    >  2:1196267/1‑150 (MQ=255)

                                                                                                                      agcGACCGTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTg                                 <  1:326223/148‑1 (MQ=255)

                                                                                                                           ccTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTgcagca                            <  2:3594373/150‑1 (MQ=255)

                                                                                                                           ccTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTgcagca                            <  2:2797118/150‑1 (MQ=255)

                                                                                                                             ttGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGt                          <  1:981607/150‑1 (MQ=255)

                                                                                                                              tGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGtt                         <  2:1635318/150‑1 (MQ=255)

                                                                                                                                    gCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATTTCGCCTGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGATTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCt                   <  1:1392463/150‑1 (MQ=255)

                                                                                                                                    gCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCt                   <  1:5616573/150‑1 (MQ=255)

                                                                                                                                      cgcgCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTgg                 <  1:503197/150‑1 (MQ=255)

                                                                                                                                       gcgcCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGa                <  1:5390128/150‑1 (MQ=255)

                                                                                                                                       gcgcCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGAAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGa                >  2:4373714/1‑150 (MQ=255)

                                                                                                                                          ccTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGc             <  1:415045/150‑1 (MQ=255)

                                                                                                                                          ccTCGGCCACAGCTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGc             >  2:2348253/1‑150 (MQ=255)

                                                                                                                                            tCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGccc           >  1:2407963/1‑150 (MQ=255)

                                                                                                                                            tCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGGTGGACCTGGATGccc           >  2:5678978/1‑150 (MQ=255)

                                                                                                                                               gCCAGGTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGcc        <  1:1486599/150‑1 (MQ=255)

                                                                                                                                               gCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGcc        >  2:5580031/1‑150 (MQ=255)

                                                                                                                                                ccAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCt       <  1:3093285/150‑1 (MQ=255)

                                                                                                                                                ccAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCGTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCt       <  1:1131763/150‑1 (MQ=255)

                                                                                                                                                     tgtgGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGACa  >  1:2741411/1‑150 (MQ=255)

                                                                                                                                                     tgtgGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGACa  >  2:884785/1‑150 (MQ=255)

GCTCAAGGCCGAGACCCAGGGGGCAGGGCAAAACTGGACCCTGGGCACACTGGCCATCCAGCTGGGTGACAACCGCATCAACGGCAGTGGCAGCTTGCAACAGCGCCTGGCTGGGCGAATCGACCTTGACTTGCCGCGCCTCGGCCAGCTGTGGCCGAGGCTACAAGGCCGGGTCAAGGGCCGGCTGGATGTCGCCGGCACCTTGCAGGCCCCGCAAGGTACGCTCACCCTGCAGGGCC‑AGCGCTTGGCCCAGGCCGAAAACCGCCTGCAGCAGTTGGACCTGGATGCCCGCCTGGACA  >  NC_002947/3297166‑3297464

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