BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,099,217	A→C	5.3%	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,099,217	0	A	C	5.3%	300.8 / 6.3	97	intergenic (+238/+54)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base A (70/21); new base C (3/2); total (74/23)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

Only 100 of 105 total aligned reads displayed.

TCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTA  >  NC_002947/5099068‑5099366

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

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

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGaa                                                                                                                                                       >  1:1098971/1‑150 (MQ=255)

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGaa                                                                                                                                                       >  2:2935683/1‑150 (MQ=255)

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATTCGTGCGGGCAGGGGACAGGG‑TAGCCTGaa                                                                                                                                                       >  1:4241837/1‑150 (MQ=255)

   gATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCGGAACTg                                                                                                                                                    >  2:2690628/1‑150 (MQ=255)

   gATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCGCCTGCGGGAGGTCCGGGGCCCGCA‑TAGCCAGAAcgg                                                                                                                                                    >  1:3499278/1‑148 (MQ=255)

    aTGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGa                                                                                                                                                   <  1:4322287/150‑1 (MQ=255)

    aTGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGa                                                                                                                                                   <  1:2277226/150‑1 (MQ=255)

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

      gCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGGCCGCGAAAGCGGGCGGTCAGGGGGCAGCG‑TAGCCGGAAcataac                                                                                                                                                 >  1:1998430/1‑145 (MQ=255)

        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCACGGGACAGCT‑TAGCCGGAACTGACGCa                                                                                                                                               >  2:2521189/1‑150 (MQ=255)

        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCCAAATCGTGAGGTAAGGGGACAGCT‑TAGCCAGAACAGACGCa                                                                                                                                               >  1:2889969/1‑150 (MQ=255)

        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TTGCCTGAACTGACGCa                                                                                                                                               >  2:4371589/1‑150 (MQ=255)

        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATGCGTGAGGTATGGGGTCAGGT‑TTACCTGTACAGAAGCa                                                                                                                                               >  1:3939740/1‑150 (MQ=255)

        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGGAACAATGCAAAAGCGGGCTTGCCCGCGCATGCGTGAGGTCGGGGGACAGCG‑TAGCCGGGACTGACGcg                                                                                                                                               >  2:2635662/1‑149 (MQ=255)

         ccGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCAt                                                                                                                                              <  2:3033259/150‑1 (MQ=255)

         ccGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCAt                                                                                                                                              <  2:3033331/150‑1 (MQ=255)

          cGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGCATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCagt                                                                                                                                             >  2:3285381/1‑148 (MQ=255)

          cGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCAtt                                                                                                                                             >  2:3285375/1‑150 (MQ=255)

             aaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCA‑TAGCCTGAACTGACGCATTagc                                                                                                                                          >  1:1815962/1‑150 (MQ=255)

              aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgcg                                                                                                                                         >  2:3850138/1‑150 (MQ=255)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                        <  1:3111256/150‑1 (MQ=255)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                        <  2:4226408/150‑1 (MQ=255)

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

                 ggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTTCCCGCGCATGCGGGAGGTCCGGGGACAGCT‑TAGCCTGAACTGACGCATTAGCGGGc                                                                                                                                      >  1:4213106/1‑150 (MQ=255)

                   ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGGCAGCT‑TAGCCGGGACTGTCGCATTCGCGGGcca                                                                                                                                    >  1:4005280/1‑148 (MQ=255)

                    cGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCGAg                                                                                                                                   >  2:3475970/1‑150 (MQ=255)

                     ggAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGc                                                                                                                                  <  1:4343698/150‑1 (MQ=255)

                       aaTTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGccc                                                                                                                                >  2:1698864/1‑150 (MQ=255)

                           gTTTCCTGGCGCCCCGAACCGGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCtt                                                                                                                            <  2:3152516/150‑2 (MQ=18)

                            tttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTGAGGGGACAGCT‑TAGCCTGAACGGACGCATTCGCGGGCAAGCCCGCGcc                                                                                                                           >  1:302277/1‑150 (MQ=255)

                            tttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCTTGCCCGcgaa                                                                                                                           <  1:3927725/150‑4 (MQ=255)

                             ttCCTGGCGCCCCGAACCTGGGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAAGGCAAACGCGGGCTTGCCCGCTAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCtttt                                                                                                                          <  2:512336/150‑4 (MQ=17)

                              tCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGAGTGAGGTGAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTg                                                                                                                         >  2:4312093/1‑150 (MQ=255)

                               tcaccggcCCCGCACAGGTTCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGc                                                                                                                        >  2:282165/7‑150 (MQ=39)

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

                                   gcgcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAggg                                                                                                                    >  2:3297975/1‑150 (MQ=255)

                                   gcgcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAggg                                                                                                                    >  2:1460867/1‑150 (MQ=255)

                                      ccccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGt                                                                                                                 >  2:3057568/1‑150 (MQ=255)

                                       cccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGCCAGCA‑TAGCCTGAACTGTCGCATTCGCGGGCCAGCCCGCTCCTGCCGGGGGtg                                                                                                                >  2:3585178/1‑149 (MQ=255)

                                       cccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGGCGCATTAGCGGGCAAGCCCGCTCCTGCAGGGGGTa                                                                                                                >  2:2141093/1‑150 (MQ=255)

                                       cccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTa                                                                                                                >  2:4463592/1‑150 (MQ=255)

                                       cccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAAAGCGGGGGGTCAGGGGACAGCG‑TAGCCGGAAATGACGCAATAGCGGGCACCCCCGCGCCGGCCgggggga                                                                                                                >  1:3624451/1‑148 (MQ=255)

                                        ccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgcg                                                                                                                                         <  1:3613533/124‑1 (MQ=255)

                                        ccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgcg                                                                                                                                         >  2:3613533/1‑124 (MQ=255)

                                        ccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGGGAGGTCAGGGGCCCGCG‑TAGCCGGAACTGACGCAGTCGCGGGCCAGCCCGCGCCGGCCGGGGGGAt                                                                                                               >  2:3089726/1‑150 (MQ=255)

                                            cccTGGGGGGGGCCGGTGACAGAACTGTACCCGATGATGGTCGGGTAACACGGCAACCGGGGGCTGGCCCGCGATTGGGTGTGGTCCGGGGACCGCT‑TCGCCTGAACTGCCGCCTGCGGGGGCAAGCCCGCTCCTGCAGGGGGTATTGcc                                                                                                           <  2:1089208/149‑1 (MQ=255)

                                                        ccGGTGACAGAACCGTACCCGTTGAGTGTCGGGTAACACTGCAAACGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTa                                                                                               <  2:4433370/150‑1 (MQ=39)

                                                         cGGTGACAGAACAGTAACCGATGATTGTCGGGTAACCATGCAAACGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TCGCCTGAACTGACGCATGCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTAc                                                                                              <  2:2574373/150‑1 (MQ=255)

                                                              aCAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGCATGCGTGAGGTCCGGGGACAGCG‑TAGCCTGTACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGa                                                                                         >  2:3525830/1‑150 (MQ=255)

                                                                aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGGTTGGCCGGGTATGGGTTTGGTGATGGGTGATCT‑TTTGCTGTAATGTCGCATTCGCGGGGAAGCCCGCTGCTGGAGGGGGTATTGGGGGGTAGGGTTTTCGCCGAAg                                                                                       >  2:1587581/1‑150 (MQ=255)

                                                                aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCAGGGGACAGCT‑TATCCTGAACTGACGCATTAGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAg                                                                                       >  2:1269354/1‑150 (MQ=255)

                                                                         aCCGATGATTGTCGGGTAACAATGCAAAAGCGGGCATGGCCGCGTAATCGTGTGGTGAGGGGACATCT‑TTGCCGGTGATGTCGCATTAGACGGCAAGAACGCACCTGCAGGGGGTATTGCCGGGAAGGGGTTACTACCAAGACACGTTCa                                                                              >  2:1590741/1‑150 (MQ=255)

                                                                          ccGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCCCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCGGAACTGACGCATTCGCGGGCCAGCCCGCTCCAGCAGGGGGGATTTCCGGGGAGGGTTTACGCCGAAGCCAGGTTCaa                                                                             >  1:326159/1‑149 (MQ=255)

                                                                                    tCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATTCGTGTGGTCAGGGGGCACCT‑TAGCCTGGACTGACGCAATCGCGGGCAAGCCCCCTCCTGCAGGGGGTATTTCCGGGAACGGTTTAAGCCCAAGCCACGTTCAACGCCGTGTgg                                                                   >  1:1416448/1‑150 (MQ=255)

                                                                                      gggTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGGGGTGAGGGGACCGCT‑TAGCCTGAACGGACGCAATAGCGGGCAAGCCCGCTCCGGCAGGGGGTATTGGCGGGAAGGGTTTACGCCGGAGCCAGGTTCATAGCCGTTTGGGt                                                                 >  1:4031300/1‑150 (MQ=255)

                                                                                        gTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGGCCGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCCCGCCTGCCGGGGGGATTGCCGGGCAGGGGTTACGCCGAAGCCAGGGTCATCGCCGTGTGGGtgt                                                               >  1:3333626/1‑148 (MQ=255)

                                                                                             aaTGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCGGAACTGACGCATTAGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGCAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATc                                                          >  2:280507/1‑150 (MQ=255)

                                                                                                 cAAAAGCGGGCTTGCCCGCCAAAGCGTGAGGTCAGGGGCCACCG‑TTGCCGGGACTGACGCATTAGCGGGCAAGCCCCCTCCGGCACGGGGGATTGCCGGGACGGGTTTACGCCCAAGCCCGGGTCATCGCCATGTGGGTATAGGTAAGGt                                                      >  2:2417960/1‑150 (MQ=255)

                                                                                                   aacgggggCTTGCCCGCTACTGCGTGTGTTCTGGGGACAGCT‑TCGCCTGAACTGACGCATTCGGGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGc                                                    <  2:681919/145‑1 (MQ=255)

                                                                                                       ggggCTGGCCCGCTATTGGGTGTGGCCGGGGTACTGTT‑TCGCCTCAACGCCCGCTTGCGGGGCCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTa                                                <  2:1998430/149‑1 (MQ=255)

                                                                                                       cGGGCTTGCCCGCGAATGCGTGAGGTCCGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTa                                                >  1:1173471/1‑150 (MQ=255)

                                                                                                        gggCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTAc                                               >  1:4505534/1‑150 (MQ=255)

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

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

                                                                                                               cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                                        >  1:3858199/1‑150 (MQ=255)

                                                                                                               cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTAATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                                        >  1:775556/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                                                                 cgcgAATGCGTGAGGTCCGGGGACAGCT‑TCGCCTGAACTGACGCATGCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg                                      <  2:4507294/150‑1 (MQ=255)

                                                                                                                  gcgAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg                                     >  2:3095308/1‑150 (MQ=255)

                                                                                                                  gcgAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg                                     >  2:3140252/1‑150 (MQ=21)

                                                                                                                      aTGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAg                                 >  2:3711183/1‑150 (MQ=255)

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

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

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

                                                                                                                           tGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAg                            >  1:901065/1‑150 (MQ=255)

                                                                                                                           tGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAg                            >  1:901041/1‑150 (MQ=255)

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

                                                                                                                               gTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCg                        >  2:45309/1‑150 (MQ=255)

                                                                                                                                tctggttaCAGCT‑TAGCCTGAACTGACGATTTCGGGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                       <  2:4241837/143‑1 (MQ=255)

                                                                                                                                tCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                       >  1:3943235/1‑150 (MQ=255)

                                                                                                                                  aGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAACCAAGGTCgaga                     >  1:2322367/1‑150 (MQ=255)

                                                                                                                                   ggggACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGat                    >  1:429435/1‑150 (MQ=255)

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

                                                                                                                                      gACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTAGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATc                 >  1:3801169/1‑150 (MQ=255)

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

                                                                                                                                                gCCTGAACTGCCGCATTCGGGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGc        <  2:2889969/150‑1 (MQ=255)

                                                                                                                                                  cTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgt      >  1:3511947/1‑150 (MQ=255)

                                                                                                                                                  cTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgt      >  1:1798756/1‑150 (MQ=255)

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

                                                                                                                                                    gAACTGCCGCTTGCGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtcg    <  2:2884450/150‑1 (MQ=255)

                                                                                                                                                      aCTGAGGCATTCGCGGGCATGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTa  <  2:1815962/150‑1 (MQ=255)

                                                                                                                                                      aCTGACGCCTTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTa  <  1:1424855/150‑1 (MQ=255)

                                                                                                                                                      aCTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTa  <  2:2418374/150‑1 (MQ=255)

TCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTA  >  NC_002947/5099068‑5099366

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