BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,443,500	T→A	6.4%	F6I (TTT→ATT) ‡	PP_2140 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,443,500	0	T	A	6.4%	301.9 / 5.3	99	F6I (TTT→ATT) ‡	PP_2140	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (45/47); new base A (5/2); total (50/49)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GGTGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAA  >  NC_002947/2443352‑2443643

ggTGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCt                                                                                                                                                  >  2:6986/1‑149 (MQ=255)

  tGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCttt                                                                                                                                                >  4:113293/1‑149 (MQ=255)

    aTCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTTCAGTCCTATGTCGATGGCAAGCCGTCCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGc                                                                                                                                              <  2:123451/149‑1 (MQ=255)

     tCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGc                                                                                                                                              <  4:127942/148‑1 (MQ=255)

       cGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGACGACGATGGCAAGCCGACCGGCTGGAAGGCATTCGACGACGGCAAAGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGc                                                                                                                                           >  8:282383/1‑149 (MQ=255)

                 gCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCAtt                                                                                                                                 <  3:63977/149‑1 (MQ=255)

                   aaGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTa                                                                                                                               >  4:183521/1‑149 (MQ=255)

                      acgaaGGAGCAGTATGTGCAGTCCGTGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTa                                                                                                                               <  6:207003/143‑1 (MQ=255)

                      aCCAAGGAGCAGTAAGTGCAGTCCGAGGTCGATGGCAAGCCGACCCGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTCTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTa                                                                                                                               >  5:207003/1‑146 (MQ=255)

                        cAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTatgatg                                                                                                                          >  8:51556/1‑149 (MQ=255)

                        cAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTatgatg                                                                                                                          >  3:97893/1‑149 (MQ=255)

                          aGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCa                                                                                                                        <  3:183521/149‑1 (MQ=255)

                                  tATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGc                                                                                                                         >  1:229740/1‑140 (MQ=255)

                                  tATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGc                                                                                                                         <  2:229740/140‑1 (MQ=255)

                                  tATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCg                                                                                                                >  2:179660/1‑149 (MQ=255)

                                   aTGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCAAAGCGGccc                                                                                                                                         >  4:334718/1‑123 (MQ=255)

                                     gtgCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCt                                                                                                             >  4:258553/1‑149 (MQ=255)

                                      tgCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAATATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCtt                                                                                                            <  2:239067/149‑1 (MQ=255)

                                      tgCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCtt                                                                                                            <  5:217481/149‑1 (MQ=255)

                                       gCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTg                                                                                                           >  8:268055/1‑149 (MQ=255)

                                        cAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTg                                                                                                           >  1:320334/1‑148 (MQ=255)

                                         aGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCTGCCCCAATCATTTATGATGCAGATATCCGCCTTGTa                                                                                                         <  2:89198/149‑1 (MQ=255)

                                         aGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTa                                                                                                         <  7:86463/149‑1 (MQ=255)

                                          gTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAg                                                                                                        >  8:168470/1‑149 (MQ=255)

                                           tCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGc                                                                                                       >  7:261859/1‑149 (MQ=255)

                                                  tCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAATGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATgg                                                                                                <  5:185729/149‑1 (MQ=255)

                                                    gATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATgg                                                                                                >  2:347552/1‑147 (MQ=255)

                                                    gATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATgg                                                                                                >  4:177636/1‑147 (MQ=255)

                                                     aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATTTCCGCCTTGTAGCCTTATGGGAg                                                                                             >  1:44904/1‑149 (MQ=255)

                                                     aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATgg                                                                                                >  7:26015/1‑146 (MQ=255)

                                                     aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGa                                                                                              <  8:26015/148‑1 (MQ=255)

                                                     aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAg                                                                                             <  1:6986/149‑1 (MQ=255)

                                                      tGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAgg                                                                                            <  2:320334/149‑1 (MQ=255)

                                                       ggCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCAAATCATTTATGATGCac                                                                                                                       >  2:119507/1‑120 (MQ=255)

                                                        gCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAggg                                                                                           >  5:318481/1‑148 (MQ=255)

                                                           aGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACt                                                                                       <  8:261859/149‑1 (MQ=255)

                                                                aCCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTCGTAGCCTTATGGGAGGGCACTGAAAt                                                                                  <  1:179660/149‑1 (MQ=255)

                                                                 ccGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATg                                                                                 <  7:168470/149‑1 (MQ=255)

                                                                 ccGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATg                                                                                 <  3:177636/149‑1 (MQ=255)

                                                                      tGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCa                                                                            <  3:113293/149‑1 (MQ=255)

                                                                        gAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑g                                                                          >  5:79426/1‑148 (MQ=255)

                                                                          aGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGaa                                                                       <  4:97893/149‑1 (MQ=255)

                                                                           ggCATTCTACGACGGCAATGCCTGGAAGGTTAAAGATAAGCGCTGATCCAATGGTGCTGTGATTTGGGGCCGCATTGCGGACCCAATAATTAAAAATGAAGAAAAACGCCTAGTAGCCATATGGAAGGGCAAGGAAATGGCCCA‑GGaa                                                                       >  7:221189/1‑148 (MQ=255)

                                                                             cATTCTACGACGGCAATGCCTGGAATTTTGAAGATAAGCGGTGATCCATTTGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGATGGCACTGAAATGTCCCA‑GGAActc                                                                    <  2:44904/149‑1 (MQ=255)

                                                                              aTTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAActc                                                                    >  2:91916/1‑148 (MQ=255)

                                                                                tCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTa                                                                  >  2:4075/1‑148 (MQ=255)

                                                                                  tACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTacac                                                               >  1:133058/1‑149 (MQ=255)

                                                                                   acgacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAcc                                                              >  7:200296/1‑149 (MQ=255)

                                                                                    cgacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAccc                                                             <  8:200296/149‑1 (MQ=255)

                                                                                    cgacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAccc                                                             <  6:318481/149‑1 (MQ=255)

                                                                                    cgacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAccc                                                             >  2:330959/1‑149 (MQ=255)

                                                                                     gacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGcc                                                                                                      <  3:237446/108‑1 (MQ=255)

                                                                                     gacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGcc                                                                                                      >  4:237446/1‑108 (MQ=255)

                                                                                     gacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCg                                                            >  2:167428/1‑149 (MQ=255)

                                                                                       cgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTa                                                          >  1:102386/1‑149 (MQ=255)

                                                                                       cgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTa                                                          <  1:167428/149‑1 (MQ=255)

                                                                                               ccTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTaccaacca                                                   <  1:347552/148‑1 (MQ=255)

                                                                                                cTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTac                                                                 >  8:71727/1‑133 (MQ=255)

                                                                                                cTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTac                                                                 <  7:71727/133‑1 (MQ=255)

                                                                                                 tGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGa                                                 >  3:250477/1‑148 (MQ=255)

                                                                                                 tGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGa                                                 >  7:171161/1‑148 (MQ=255)

                                                                                                  ggAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGaa                                                >  8:275305/1‑148 (MQ=255)

                                                                                                   gAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGATCTCTACACCCGTACCAACCAGaaa                                               <  7:51556/148‑1 (MQ=255)

                                                                                                   gAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGaaa                                               <  4:250477/148‑1 (MQ=255)

                                                                                                   gAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGaaa                                               <  8:171161/148‑1 (MQ=255)

                                                                                                   gAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAAc                                              >  1:62610/1‑149 (MQ=255)

                                                                                                       gTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGt                                           >  7:286214/1‑148 (MQ=255)

                                                                                                         tGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTaccaacca                                                   <  6:18584/138‑1 (MQ=255)

                                                                                                         tGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACAACCGTACCAACCAGagatcggg                                          >  5:18584/1‑140 (MQ=255)

                                                                                                             gATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGc                                    >  4:14450/1‑149 (MQ=255)

                                                                                                               tAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCg                                  <  2:62610/149‑1 (MQ=255)

                                                                                                               tAAGCGCTGATCCAATGGTGCTGTGAATTGGGGCCGCTTTGCGGCCCCAATCATTTAAGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGAACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTATACGCCg                                  >  7:335207/1‑149 (MQ=255)

                                                                                                                  gcgcTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTaccaacca                                                   <  2:135549/129‑1 (MQ=255)

                                                                                                                  gcgcTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTaccaacca                                                   >  1:135549/1‑129 (MQ=255)

                                                                                                                      tGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCg                                  >  3:180145/1‑142 (MQ=255)

                                                                                                                      tGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCg                                  <  4:180145/142‑1 (MQ=255)

                                                                                                                           cATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                       <  3:14450/148‑1 (MQ=255)

                                                                                                                           cATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                       <  2:102386/148‑1 (MQ=255)

                                                                                                                            aTTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                       >  6:63150/1‑147 (MQ=255)

                                                                                                                            aTTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  6:79426/149‑1 (MQ=255)

                                                                                                                            aTTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  7:275305/149‑1 (MQ=255)

                                                                                                                            aTTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  2:133058/149‑1 (MQ=255)

                                                                                                                             ttGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  1:91916/149‑1 (MQ=255)

                                                                                                                                gTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTc                  <  1:330959/148‑1 (MQ=255)

                                                                                                                                 tGCTGTGATTTGGTGTCGCTTTGCGGCCCAAATTATTTATGTTGGAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGTAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAAGTGTTCTTCGCGGGCCTCGCCCTCGAGTcc                 <  8:335207/148‑1 (MQ=255)

                                                                                                                                    tgtgATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATgg             >  6:162402/1‑149 (MQ=255)

                                                                                                                                    aaatgattGGGGCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     >  1:225098/7‑141 (MQ=255)

                                                                                                                                    aaatgattGGGGCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  2:225098/135‑1 (MQ=255)

                                                                                                                                      tgATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGc            <  3:258553/148‑1 (MQ=255)

                                                                                                                                        aTTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATgg             >  1:79042/1‑145 (MQ=255)

                                                                                                                                        aTTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATgg             <  2:79042/145‑1 (MQ=255)

                                                                                                                                           tGGGGCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGc      >  2:104446/1‑149 (MQ=255)

                                                                                                                                             gggCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCt     <  8:286214/148‑1 (MQ=255)

                                                                                                                                             gggCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     >  2:93034/1‑132 (MQ=255)

                                                                                                                                             gggCCGCAAAGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGa                     <  1:93034/132‑1 (MQ=255)

                                                                                                                                              ggCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTg    >  2:4124/1‑148 (MQ=255)

                                                                                                                                               gCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGaa  >  5:297590/1‑149 (MQ=255)

                                                                                                                                                  gCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                       <  6:153422/125‑1 (MQ=255)

                                                                                                                                                   cTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑gg                                                                         >  2:330136/1‑74 (MQ=255)

                                                                                                                                                   cTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑gg                                                                         <  1:330136/74‑1 (MQ=255)

GGTGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCCAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCA‑GGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAA  >  NC_002947/2443352‑2443643

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