BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,057,685	G→C	5.7%	L140V (CTG→GTG)	PP_3572 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,057,685	0	G	C	5.7%	831.1 / 29.7	244	L140V (CTG→GTG)	PP_3572	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (123/107); new base C (9/5); total (132/112)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 250 total aligned reads displayed.

GATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGT  >  NC_002947/4057536‑4057834

gATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAg                                                                                                                                                       <  1:859728/150‑1 (MQ=255)

 aTCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAgg                                                                                                                                                      <  1:2949036/150‑1 (MQ=255)

     ggTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGgc                                                                                                                                                  >  2:5507540/1‑150 (MQ=255)

      gTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGgcc                                                                                                                                                 >  1:274370/1‑150 (MQ=255)

        cGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGgccgc                                                                                                                                               >  2:3298877/1‑150 (MQ=255)

         gCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCCGGCGGCCACCTGgccgac                                                                                                                                              >  2:3450461/1‑148 (MQ=255)

         gCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCCGGCGGCCACCTGgccgac                                                                                                                                              >  2:637704/1‑148 (MQ=255)

           tGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCgg                                                                                                                                            <  1:1415202/150‑1 (MQ=255)

            gTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGc                                                                                                                                           <  2:5471420/150‑1 (MQ=255)

            gTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGc                                                                                                                                           <  2:4089786/150‑1 (MQ=255)

               ggCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCAATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTg                                                                                                                                        >  2:4678278/1‑150 (MQ=255)

                gCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGt                                                                                                                                       >  1:1018909/1‑150 (MQ=255)

                gCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGt                                                                                                                                       >  2:161356/1‑150 (MQ=255)

                gCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGt                                                                                                                                       >  1:1818627/1‑150 (MQ=255)

                gCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGt                                                                                                                                       >  1:1018936/1‑150 (MQ=255)

                 cccAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGtt                                                                                                                                      >  1:2725701/1‑150 (MQ=255)

                  ccAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCCGGCGGCCACCTGGCCGACCGCCTgctg                                                                                                                                     <  1:637704/150‑4 (MQ=255)

                    aaTCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGttggc                                                                                                                                   <  2:5064844/150‑4 (MQ=255)

                    aaTCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaa                                                                                                                                   >  2:254097/1‑150 (MQ=255)

                     aTCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaa                                                                                                                                  <  2:814775/150‑1 (MQ=255)

                      tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                 >  1:4976794/1‑150 (MQ=255)

                      tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                 >  2:1501042/1‑150 (MQ=255)

                      tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                 <  1:3209187/150‑1 (MQ=255)

                      tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                 >  1:281519/1‑150 (MQ=255)

                      tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                 <  2:3982030/150‑1 (MQ=255)

                      tCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                 <  1:1867113/150‑1 (MQ=255)

                         aCGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGcc                                                                                                                              <  2:3521525/150‑1 (MQ=255)

                           gcgcAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAg                                                                                                                            >  1:2614626/1‑150 (MQ=255)

                              cAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCAc                                                                                                                         <  1:5507540/150‑1 (MQ=255)

                                 gTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCaaaa                                                                                                                                 >  2:1292060/1‑139 (MQ=255)

                                 gTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCAt                                                                                                                      >  2:1908977/1‑150 (MQ=255)

                                 gTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCAt                                                                                                                      >  2:361314/1‑150 (MQ=255)

                                  tCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATg                                                                                                                     >  1:742110/1‑150 (MQ=255)

                                     accaGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCt                                                                                                                  >  2:1344798/1‑150 (MQ=255)

                                        aGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCg                                                                                                               <  1:3088539/150‑1 (MQ=255)

                                          cGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGAt                                                                                                             <  2:5247430/150‑1 (MQ=255)

                                            gCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCAt                                                                                                                      >  1:98348/1‑139 (MQ=255)

                                            gCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCa                                                                                                           >  2:3187763/1‑150 (MQ=255)

                                             ccAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCAc                                                                                                          >  2:113181/1‑150 (MQ=255)

                                             ccAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCAc                                                                                                          >  2:112801/1‑150 (MQ=255)

                                                ggTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTg                                                                                                       <  1:4607271/150‑1 (MQ=255)

                                                 gTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTgg                                                                                                      >  1:3949860/1‑150 (MQ=255)

                                                    aGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTc                                                                                                   >  2:4813285/1‑150 (MQ=255)

                                                    aGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTc                                                                                                   >  1:1484541/1‑150 (MQ=255)

                                                      cTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGa                                                                                                 <  1:2293384/150‑1 (MQ=255)

                                                       tCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAg                                                                                                <  1:2397898/150‑1 (MQ=255)

                                                       tCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAg                                                                                                <  2:742110/150‑1 (MQ=255)

                                                       tCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAg                                                                                                >  2:2048461/1‑150 (MQ=255)

                                                         gCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGAt                                                                                                             <  1:570327/135‑1 (MQ=255)

                                                         gCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGa                                                                                              <  2:4162370/150‑1 (MQ=255)

                                                           cAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAca                                                                                            <  2:281519/150‑1 (MQ=255)

                                                            aGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcag                                                                                           >  1:3817956/1‑150 (MQ=255)

                                                            aGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcag                                                                                           <  1:2142762/150‑1 (MQ=255)

                                                            aGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcag                                                                                           >  1:4816522/1‑150 (MQ=255)

                                                              cAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcagca                                                                                         <  2:2222773/150‑1 (MQ=255)

                                                              cAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcagca                                                                                         <  2:2509496/150‑1 (MQ=255)

                                                                 tCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGt                                                                                      >  1:4081281/1‑150 (MQ=255)

                                                                    tCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCAGGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGc                                                                                   >  1:5568474/1‑150 (MQ=255)

                                                                     cAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGa                                                                                              <  1:1906704/138‑1 (MQ=255)

                                                                     cAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGcc                                                                                  >  1:3377930/1‑150 (MQ=255)

                                                                      aGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGtt                                                                                     <  2:238087/146‑1 (MQ=255)

                                                                         gccgccGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAgcc                                                                              <  1:3556633/150‑1 (MQ=255)

                                                                              cgAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCcgccgc                                                                         <  1:1121329/150‑1 (MQ=255)

                                                                              cgAGTGCACGGGCCCGGCTTCACCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCcgccgc                                                                         <  1:5150078/150‑1 (MQ=255)

                                                                               gAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCcgcg                                                                        >  2:5427884/1‑150 (MQ=255)

                                                                                        ggCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCg                                                               >  2:5077869/1‑150 (MQ=255)

                                                                                           ccGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCg                                                            >  1:1066298/1‑150 (MQ=255)

                                                                                            cGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCgg                                                           >  2:4209401/1‑150 (MQ=255)

                                                                                              gCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGcc                                                         >  1:2040621/1‑150 (MQ=255)

                                                                                              gCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGcc                                                         >  1:2040305/1‑150 (MQ=255)

                                                                                                 tCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAACAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCAtt                                                      <  2:3796939/150‑1 (MQ=255)

                                                                                                    ccACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAACCGCCGCGCATTGTTCGCCGGCCATTGCa                                                   >  2:1745811/1‑150 (MQ=255)

                                                                                                      aCCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATg                                                 <  2:4857070/150‑1 (MQ=255)

                                                                                                        ccGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGgt                                               >  1:4837726/1‑150 (MQ=255)

                                                                                                              ccATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGc                                         <  1:1754376/150‑1 (MQ=255)

                                                                                                                aTCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGccc                                       >  2:1544966/1‑150 (MQ=255)

                                                                                                                 cgcccgcGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCa                                      >  1:4541417/4‑150 (MQ=255)

                                                                                                                   gccgcGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGt                                    <  1:4118314/150‑1 (MQ=255)

                                                                                                                    ccgcGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTg                                   <  1:3467905/150‑1 (MQ=255)

                                                                                                                            aTGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTcccgc                           <  1:4402968/150‑1 (MQ=255)

                                                                                                                            aTGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTcccgc                           >  2:2663011/1‑150 (MQ=255)

                                                                                                                                ttttgaacagGCCGGCGGCCACCTGGCCGACCGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCgg                       >  1:3521786/8‑150 (MQ=38)

                                                                                                                                cTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCgg                       <  2:3654377/150‑1 (MQ=255)

                                                                                                                                 tGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGt                      <  2:875704/150‑1 (MQ=255)

                                                                                                                                 tGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGt                      <  2:875837/150‑1 (MQ=255)

                                                                                                                                 tGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGt                      >  1:3341193/1‑150 (MQ=255)

                                                                                                                                    cagcagGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAca                                                                                            <  2:4616342/77‑1 (MQ=255)

                                                                                                                                     aacagGCCGGCGGCCACCTGGCCGACCGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAAc                  >  1:430505/3‑150 (MQ=255)

                                                                                                                                        agGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAgg               <  2:1002870/150‑1 (MQ=255)

                                                                                                                                         ggCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAggg              <  1:4923173/150‑1 (MQ=255)

                                                                                                                                          gCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTgg                                                                                                      <  2:2990427/61‑1 (MQ=255)

                                                                                                                                               cGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCt        >  1:3001186/1‑150 (MQ=255)

                                                                                                                                                 gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCcgccgc                                                                         <  2:3814308/83‑1 (MQ=255)

                                                                                                                                                 gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTgg      >  2:3671954/1‑150 (MQ=255)

                                                                                                                                                 gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTgg      >  2:71332/1‑150 (MQ=255)

                                                                                                                                                 gCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTgg      >  1:4883729/1‑150 (MQ=255)

                                                                                                                                                  ccAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTgg      <  2:3356291/149‑1 (MQ=255)

                                                                                                                                                   cAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGTCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGcc    >  1:3775909/1‑150 (MQ=255)

                                                                                                                                                   cAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGcc    <  1:3654581/150‑1 (MQ=255)

                                                                                                                                                     ggTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGt  >  2:3377719/1‑150 (MQ=255)

GATCGGGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCAGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGGGCCACCTGGCCGT  >  NC_002947/4057536‑4057834

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