BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,102,244	A→C	11.5%	intergenic (+238/+54)	PP_mr66 → / ← acsA‑I	P6/acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,102,244	0	A	C	11.5%	393.7 / 31.6	133	intergenic (+238/+54)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base A (92/25); new base C (8/8); total (100/33)
Fisher's exact test for biased strand distribution p-value = 2.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 137 total aligned reads displayed.

TCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGT  >  NC_002947/5102095‑5102392

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGCATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGaa                                                                                                                                                      >  2:1591392/1‑150 (MQ=255)

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGaa                                                                                                                                                      <  2:5112093/150‑1 (MQ=255)

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGaa                                                                                                                                                      <  2:5112114/150‑1 (MQ=255)

 ccGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAAc                                                                                                                                                     <  1:1516805/150‑1 (MQ=255)

  cGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTAAGGGGACAGCT‑TAGCCTGAACt                                                                                                                                                    >  1:2359182/1‑150 (MQ=255)

   gATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCGGAAcgg                                                                                                                                                   >  2:3533623/1‑148 (MQ=255)

     tGCACCGAAACGGGCCTGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  1:379159/150‑1 (MQ=255)

     tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  1:3426581/150‑1 (MQ=255)

     tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  2:4801807/150‑1 (MQ=255)

           gAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGATGCATTc                                                                                                                                           <  2:3570267/150‑1 (MQ=255)

           gAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTc                                                                                                                                           <  2:6545201/150‑1 (MQ=255)

             aaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgc                                                                                                                                         <  1:6129184/150‑1 (MQ=255)

              aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgcg                                                                                                                                        <  1:811573/150‑1 (MQ=255)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCTATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                       <  1:4217622/150‑1 (MQ=255)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCGGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGACTGCGTGAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCgg                                                                                                                                       <  2:1666942/150‑1 (MQ=21)

                gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCggg                                                                                                                                      <  2:6048415/150‑1 (MQ=255)

                   ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGTCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCaa                                                                                                                                   <  1:1464818/150‑1 (MQ=255)

                   ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGcca                                                                                                                                   >  2:5653923/1‑148 (MQ=255)

                    cGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCAGGGGCCAGCT‑TTGCCGGTCCGGACGCATTCGCGGGCacg                                                                                                                                  >  1:6745363/1‑148 (MQ=255)

                        tttGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACTTTAACCGATGATTGTGGGTTAACAATGCAAACGGGGGCTTGCCCGCGATTGCGTGCGGTCAGGGGACAGCT‑TAGCCTTAACTGACGCATTCGCGGGCAAGCCCg                                                                                                                              <  2:2579606/149‑1 (MQ=255)

                            tttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTCACAATGCAAAAGCGGGCGTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACGGCCGCCGTCGCGGGCCAGcccccgcc                                                                                                                          >  2:3938585/1‑145 (MQ=255)

                            tttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGTGCAAGccaccacc                                                                                                                          >  2:1082530/1‑144 (MQ=18)

                                cTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCGTTCGCGGGc                                                                                                                                     >  1:575147/1‑135 (MQ=255)

                                cTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCGTTCGCGGGc                                                                                                                                     <  2:575147/135‑1 (MQ=255)

                                 tGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAg                                                                                                                     >  2:4894390/1‑150 (MQ=255)

                                  ggCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCTGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTTTTGCatt                                                                                                                    <  2:1193785/150‑3 (MQ=17)

                                     gcCCCGCACAGGTTCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGCATGCGGGGGGTCAGGGGCCAGCG‑TAGCCGGCACGGACCCAATCGCGGGCCAGCCCCCGCCGGCAggggg                                                                                                                 >  1:2109142/1‑150 (MQ=255)

                                     gcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGGTTGCCCGCGGATGCGTGTGGTGAGGGGTCAGCG‑TAGCCGGAACGGTCGCATTCGCGGGCCACCCCGCGCCGGCgggggg                                                                                                                 >  1:843405/1‑150 (MQ=255)

                                     gcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTGAGGGGGCACCG‑TAGCCAGGACTGACGCAGTCGCGGGCAAGCCCGCGCCGGCAggggg                                                                                                                 >  1:2669248/1‑150 (MQ=255)

                                     gcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGCCAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAggggg                                                                                                                 >  2:3935647/1‑150 (MQ=255)

                                       cccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGGCCCCGAATGCGTGAGGTCAGGGGGCAGCT‑TTGCCTGTACAGACGCCATAGCGGGGCAGCCCCCGCCAGCAgggggga                                                                                                               >  1:4033123/1‑148 (MQ=255)

                                           acccggtgGGGGGCGGGTGACACAACATTACCCTTTGATTGTCGGTTAAAAATGAAAACGCGGGCTTGCCCGCGATTGCGTGTGGTCCGGGGAAAACT‑TAGCCTCAACGGCCGCTTTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGc                                                                                                           <  2:5387405/145‑1 (MQ=255)

                                                     gggCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATGCGGGCGGTGGGGGGACCCCG‑TAGCCGGAACGGGCGCAGTCGCGGGCCAGCCCCCGCCGGCCGGGGGGAGTGGCGGGGAgggg                                                                                                 >  1:6563764/1‑149 (MQ=255)

                                                          ggTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTACCCCGAAATCGTGTGGTGAGGGGACAGCA‑TATCCTGTAAGGAAGCAATAGCGGGCAACCCCGCTCCTGCAGGGGGGGTTTGCGGGGAGGGTTTaag                                                                                            >  1:4585621/1‑148 (MQ=255)

                                                           gTGACAGAACAGTAAGCGATGATTGTCGGGTAACAATGCAAATGCGGGCTTGGGCGCGTATGGGTGAGGTGGGGGGTCAGGT‑TTGCCTGTAGTGACGGATTGGGGGGCAAGGCCGCTGCTGGAGGGGGTTTTGGGGGGGAGGGttttagg                                                                                           >  2:1041355/1‑146 (MQ=255)

                                                           gTGACAGAACAGTAAGCGATGATTGTCGGGTAACAATGCAAAAGCGGGGTTTCCCGCGAATGCGTGTGGTGAGGGGACAGGT‑TTGGGTGTAGTGTCGCATTGGCGGGCAAGCCCGCTCGTGGAGGGGGTTTTGGCGGGTAGGGTTTAAGc                                                                                           >  2:1041336/1‑150 (MQ=255)

                                                                        aaCCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGCGGTCCGGGGCCAGCT‑TAGCCGGAACGGACCCATTCGCGGGCAAGCCCGCTCCTGCCGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCCCGTTc                                                                              >  1:4211405/1‑150 (MQ=255)

                                                                        aaCCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTTCCCCCGCAAGCGTGAGGTCAGGGGACAGCA‑TAGCCGGAACAGACCCAATAGCGGGCAAGCCCGCACCGGCAGGGGGTATTTCCGGGAAGGGTTTAACCCGAACCCACGttg                                                                              >  1:2883454/1‑149 (MQ=255)

                                                                        aaCCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTTCCCCCGAATGCGTGAGGTGAGGGGCCAGCT‑TAGCCGGAACTGACGCAATCGCGGGCAAGCCCCCTCCGGCAGGGGGTATTTCCGGGAAGGGTTTACGCCGAAGCCAGGTTc                                                                              >  1:4211498/1‑150 (MQ=255)

                                                                          ccGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGCAAGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCACCAGCAGGGGGGATTGCCGGGAAGGGTTTACGCCGAAGCCAGGttaaa                                                                            >  1:1732435/1‑147 (MQ=255)

                                                                            gatgatTGTCGGGTAACAATGCAAATGCGGGCTTGCCCGCGAATTGGTGTGGTGAGGGGTCAGGT‑TAGGGTGAAATGACGCATTCGCGGGCAAGCCCGCTGCTGGAGGGGGTATTGCCGGGTAGGGTTTACGCCGAAGCCAGGTTAAtgg                                                                          >  2:837664/1‑148 (MQ=255)

                                                                            gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCCAAAGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAAGCATTCGCGGGCAAGCCCGCTCCAGCAGGGGGGATTGCCGGGAACGGTTTACGCCGAAGCCAGGTTCATCg                                                                          >  2:3881478/1‑150 (MQ=255)

                                                                            gatgatTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATGCGTGAGGTCAGGGGACCGCG‑TAGCCTGAACTGACCCATTCGCGGGCAAGCCCGCTCCAGCAGGGGGTATTGCCGGGACGGGTTTACGCCGAAGCCAGGTTCATCg                                                                          >  1:5452051/1‑150 (MQ=255)

                                                                                       ggTAACAATGCAAAAGCGGGCTTGCCCCCGCATGCGTGAGGTCAGGGGACAGCG‑TAGCCGGAACTGACGCATTCGCGGGCAAGCCCGCTCCGGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGGTCATCGCCGTGTgggga                                                               >  2:3534686/1‑148 (MQ=255)

                                                                                            cAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCAGGGGACAGCT‑TAGCCAGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGAt                                                          >  2:3711072/1‑150 (MQ=255)

                                                                                               tGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCAGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTCGTGGGTATCGATCAg                                                       >  2:4269643/1‑150 (MQ=255)

                                                                                                gCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGGATCGATCAgg                                                      >  2:6781936/1‑150 (MQ=255)

                                                                                                  attcgCGGGCAAGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                                                    >  1:6185254/5‑150 (MQ=255)

                                                                                                  aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                                                    >  2:7399779/1‑150 (MQ=255)

                                                                                                  aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                                                    >  2:6943031/1‑150 (MQ=255)

                                                                                                  aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                                                    >  2:3665459/1‑150 (MQ=255)

                                                                                                  aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACAGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGACAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                                                    >  2:2230659/1‑150 (MQ=255)

                                                                                                    aaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCt                                                  >  2:5540351/1‑150 (MQ=255)

                                                                                                          gCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACca                                            >  2:4530282/1‑150 (MQ=255)

                                                                                                                ccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                      >  2:5362765/1‑150 (MQ=255)

                                                                                                                ccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                      >  2:2332395/1‑150 (MQ=255)

                                                                                                                ccGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTTCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATc                                                                           >  2:3422029/1‑113 (MQ=11)

                                                                                                                ccGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                      >  1:5944621/1‑150 (MQ=21)

                                                                                                                ccGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                      >  1:673600/1‑150 (MQ=21)

                                                                                                                ccGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                      >  1:5476414/1‑150 (MQ=21)

                                                                                                                 cgcgAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg                                     >  2:468615/1‑150 (MQ=255)

                                                                                                                 cgcgAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg                                     >  1:4631461/1‑150 (MQ=21)

                                                                                                                  gcgAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg                                    <  1:3158651/150‑1 (MQ=255)

                                                                                                                  gcgAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTTTTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg                                    >  2:3477241/1‑150 (MQ=18)

                                                                                                                    gAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTc                                  >  1:3903317/1‑150 (MQ=255)

                                                                                                                         cGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCa                             >  1:1982791/1‑150 (MQ=255)

                                                                                                                          gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                            >  2:254331/1‑150 (MQ=255)

                                                                                                                          gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                            >  2:6599351/1‑150 (MQ=255)

                                                                                                                          gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                            >  1:3237932/1‑150 (MQ=255)

                                                                                                                          gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                            >  1:5146179/1‑150 (MQ=255)

                                                                                                                          gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                            >  1:5146069/1‑150 (MQ=255)

                                                                                                                            gAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAgg                          >  1:114090/1‑150 (MQ=255)

                                                                                                                             aGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGt                         >  1:2778311/1‑150 (MQ=255)

                                                                                                                              ggTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTc                        <  1:6599351/150‑1 (MQ=255)

                                                                                                                                 cAGGGTACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgag                     <  2:4846170/150‑1 (MQ=255)

                                                                                                                                 cAGGGGACCGCT‑TAGCCTTAACTGACGCATTCGGGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgag                     <  2:4846096/150‑1 (MQ=255)

                                                                                                                                  tggggACAGCT‑TCGCCTTAACTGCCGCATTCGGGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgaga                    <  2:4607456/149‑1 (MQ=255)

                                                                                                                                  aGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgaga                    >  1:6082314/1‑150 (MQ=255)

                                                                                                                                  aGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgaga                    >  1:6494487/1‑150 (MQ=255)

                                                                                                                                  aGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgaga                    >  1:6980199/1‑150 (MQ=255)

                                                                                                                                    gggACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGata                  >  1:2891907/1‑150 (MQ=255)

                                                                                                                                    gggACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGata                  >  2:2633704/1‑150 (MQ=255)

                                                                                                                                      gACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATc                >  2:5696826/1‑150 (MQ=255)

                                                                                                                                      gACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATc                >  1:3048648/1‑150 (MQ=255)

                                                                                                                                        ctgCT‑TAGCCTGACCTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGc              <  2:4184813/148‑1 (MQ=255)

                                                                                                                                         tgCT‑TTGCCTGACCTGACGCTTGCGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGcc             <  2:3372035/149‑1 (MQ=255)

                                                                                                                                         cgCT‑TAGCCTGACCGGACGCTTGCGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGcc             <  2:3447400/149‑1 (MQ=255)

                                                                                                                                         aGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGcc             >  2:6211227/1‑150 (MQ=255)

                                                                                                                                           cT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCa           >  1:2190526/1‑150 (MQ=255)

                                                                                                                                           cT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCa           >  1:2189812/1‑150 (MQ=255)

                                                                                                                                            t‑tAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAg          >  1:4816904/1‑150 (MQ=255)

                                                                                                                                            t‑tAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAg          >  1:6294575/1‑150 (MQ=255)

                                                                                                                                            t‑tAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAg          >  1:5550615/1‑150 (MQ=255)

                                                                                                                                            t‑tAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAg          >  1:3612595/1‑150 (MQ=255)

                                                                                                                                                gCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                                                    >  1:5616151/1‑105 (MQ=255)

                                                                                                                                                gCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGc       >  1:1958872/1‑150 (MQ=255)

                                                                                                                                                 ccTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcg      >  2:339830/1‑150 (MQ=255)

                                                                                                                                                   tGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtc    >  2:6235124/1‑150 (MQ=255)

                                                                                                                                                   tGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtc    <  1:6211227/150‑1 (MQ=255)

                                                                                                                                                     aaCTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtcgt  >  2:6448895/1‑150 (MQ=255)

TCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGT  >  NC_002947/5102095‑5102392

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: `‑`