BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,491,705	A→C	6.9%	intergenic (‑85/+20)	yhfG ← / ← ppiA	DUF2559 domain‑containing protein YhfG/peptidyl‑prolyl cis‑trans isomerase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,491,705	0	A	C	6.9%	323.7 / 16.1	129	intergenic (‑85/+20)	yhfG/ppiA	DUF2559 domain‑containing protein YhfG/peptidyl‑prolyl cis‑trans isomerase A
Reads supporting (aligned to +/- strand): ref base A (55/65); new base C (4/5); total (59/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

Only 100 of 132 total aligned reads displayed.

GAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAAT  >  NC_000913/3491557‑3491854

gAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑At                                                                                                                                                        >  2:710918/1‑150 (MQ=255)

 aaaaTTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGTGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑Att                                                                                                                                                       <  2:1387392/150‑1 (MQ=255)

        cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTTATCGTTTCCATTTGCCACAAGTATAAGCTGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAT‑CTTGCCcgcg                                                                                                                                                <  2:742269/150‑1 (MQ=255)

        cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAT‑CTTGCCcgcg                                                                                                                                                <  2:584273/150‑1 (MQ=255)

        cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAT‑CTTGCCcgcg                                                                                                                                                <  2:1226650/150‑1 (MQ=255)

        cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCcgcg                                                                                                                                                >  1:1579562/1‑150 (MQ=255)

         gattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCcgcgc                                                                                                                                               <  2:2103773/150‑1 (MQ=255)

          attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAT‑CTTGCCcgcgct                                                                                                                                              <  2:1027565/150‑2 (MQ=255)

          attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAT‑CTTGCCcgcgct                                                                                                                                              <  2:1024377/150‑2 (MQ=255)

          attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCcgcgcg                                                                                                                                              <  1:1052558/150‑1 (MQ=255)

            tacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCgag                                                                                                                                            <  2:2339395/150‑1 (MQ=255)

             acgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CCAG‑ATTTCCCGCGCgaga                                                                                                                                           >  2:1886557/1‑150 (MQ=255)

              cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑At                                                                                                                                                        <  1:583990/136‑1 (MQ=255)

              cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑At                                                                                                                                                        >  2:583990/1‑136 (MQ=255)

              cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGaa                                                                                                                                          >  1:1065246/1‑150 (MQ=255)

                cTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAAt                                                                                                                                        >  1:934593/1‑150 (MQ=255)

                 tGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATc                                                                                                                                       <  1:969965/150‑1 (MQ=255)

                    aGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑Att                                                                                                                                   >  2:1532642/1‑150 (MQ=255)

                    aGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑Att                                                                                                                                   >  2:786925/1‑150 (MQ=255)

                      cTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTAc                                                                                                                                 <  1:1026923/150‑1 (MQ=255)

                       tCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑GCAG‑ATTGGCCGCGCGGGGAATC‑ATTaag                                                                                                                                >  1:658236/1‑148 (MQ=255)

                       tCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACg                                                                                                                                >  2:2332707/1‑150 (MQ=255)

                        cccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑GGGG‑GCAG‑AGTTCCCGCCCGCGAAAAC‑AATAAgg                                                                                                                               >  1:580127/1‑150 (MQ=255)

                                cGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACtt                                                                                                                       >  1:1656585/1‑150 (MQ=255)

                                cGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACtt                                                                                                                       <  2:655321/150‑1 (MQ=255)

                                   gACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑GAAG‑ATTTGCCCCGCGCGAAATA‑AATAAGGCCgggcgtttt                                                                                                                    >  1:2455825/1‑143 (MQ=255)

                                   gACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAg                                                                                                                    <  2:1928275/150‑1 (MQ=255)

                                    aCTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGc                                                                                                                   <  2:1276404/150‑1 (MQ=255)

                                     cTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑GGGGCCAAA‑ATTGCCC‑CGCGCGAAATA‑AATTAGGCAGGACgttttag                                                                                                                  >  2:592402/1‑146 (MQ=255)

                                      ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCgg                                                                                                                 <  1:1882647/150‑1 (MQ=255)

                                       tttGCTTGTCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATA‑ATTACGGCAGGACTTTAGCGGa                                                                                                                >  2:46184/1‑150 (MQ=255)

                                         tGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CCGG‑CAAG‑AATGCCCGCCCCAGAAAAC‑AATAAGGCAGGAAATTAGGGGaaa                                                                                                              >  2:2179050/1‑150 (MQ=255)

                                          gCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CCGG‑GGAG‑CTTGGCCGCGCGCGGAATC‑ATTTCGGCCGGCCGTTTGCGGGAAg                                                                                                             >  1:1149538/1‑150 (MQ=255)

                                            ttATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGAGGGG‑CATA‑ATTGCCC‑CCCGCCAAATCAAATACGG‑ACCACTTTTGCGGAAAGG‑a                                                                                                          >  1:2207815/1‑150 (MQ=255)

                                              aTCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATa                                                                                                        <  2:330450/150‑1 (MQ=255)

                                               tCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑AATGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTATCGGAAACG‑aaaa                                                                                                       >  2:1234740/1‑150 (MQ=255)

                                                   tGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTAc                                                                                                   <  1:786925/150‑1 (MQ=255)

                                                      gTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCgg                                                                                                >  2:1835847/1‑150 (MQ=255)

                                                      gTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCgg                                                                                                <  1:1488678/150‑1 (MQ=255)

                                                         tcttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGttt                                                                                             <  2:837044/150‑1 (MQ=255)

                                                         tcttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCCCGAGAAATA‑ATTAAGGCAGGACTTTATGAGAAAGG‑TTAACAAAAGGttt                                                                                             >  2:1285177/1‑150 (MQ=255)

                                                              aCTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACg                                                                                        >  2:1143539/1‑150 (MQ=255)

                                                              aCTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCCGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACg                                                                                        <  1:358680/150‑1 (MQ=255)

                                                                ttGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGc                                                                                      <  1:2321493/150‑1 (MQ=255)

                                                                 tGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGca                                                                                     <  2:271568/150‑1 (MQ=255)

                                                                    tGCCTCCCTGAACGTGTCAATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAA‑ATAGACCCCGCGAGAAATC‑ATTACGGCAGGACTTTAGCAGAAACG‑ATAACAACCGGTTTAGAAGGcacaa                                                                                  >  2:1774181/1‑149 (MQ=255)

                                                                           cTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTa                                                                           <  1:1515401/150‑1 (MQ=255)

                                                                            tGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑GAAG‑ATTGCCCGCGCGAGAAATC‑AATACGGCAGGGCGTTTGGGGGAAGG‑GTTACGACCGGGTTTGACGGGACATTCTGGGAc                                                                          >  1:1552990/1‑150 (MQ=255)

                                                                            tGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGACCGCGCGAGAAAAC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACAACACGTTTTGACGGCACATTCTGGTAc                                                                          >  2:1662320/1‑150 (MQ=255)

                                                                              aaCGTGTCCATTTTCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACgg                                                                        <  1:1215640/150‑1 (MQ=255)

                                                                                cGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTAAGGCAGGACGTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTAAGGAc                                                                      >  1:379052/1‑150 (MQ=255)

                                                                                  tgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCa                                                                    >  2:1479195/1‑150 (MQ=255)

                                                                                  tgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCa                                                                    >  1:1403066/1‑150 (MQ=255)

                                                                                    tCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTAAGGCAGGACATTAGCGGAAAGG‑ATAACTACCGGGTTTGACGGGCCATTCTGGTACGGACCAAc                                                                  >  1:435575/1‑150 (MQ=255)

                                                                                     ccATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTTCCCGCGCGAGAAATA‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTTACGGCACATTATGGTACGGACCAAcc                                                                 >  2:1245069/1‑149 (MQ=255)

                                                                                     ccATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACg                                                                 >  2:1696326/1‑150 (MQ=255)

                                                                                       aTTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAGAATTGACC‑CCCCAGAAAAA‑AATAAGGCAGGACATTAGGGGAAAGGAAAAAAAAACGGGTTTTA‑GGCACAATATGGTTAGGACCAACGTc                                                               >  2:249678/1‑150 (MQ=255)

                                                                                       aTTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTc                                                               >  1:2048223/1‑150 (MQ=255)

                                                                                        tttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCa                                                              <  2:558307/150‑1 (MQ=255)

                                                                                         ttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCAt                                                             <  2:1318587/150‑1 (MQ=255)

                                                                                          tGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                                            <  2:1579562/150‑1 (MQ=255)

                                                                                          tGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                                            <  1:1369476/150‑1 (MQ=255)

                                                                                            ccACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAg                                                          >  2:1430084/1‑150 (MQ=255)

                                                                                             cacaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGt                                                         >  1:1524999/1‑150 (MQ=255)

                                                                                             cacaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTAAGGCAGGACATTTGCGGAAAGG‑ATAACTACAGGGTTTGACGGCACATTCTGGTACGGGCCAAAGGGATGAGt                                                         >  2:1257902/1‑150 (MQ=255)

                                                                                                aaGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCgg                                                      <  1:1532642/150‑1 (MQ=255)

                                                                                                 aGTATAAGCGTCAAATCCGGGTTACGCCGGGATAGCAGAAGCG‑CTGG‑CAAT‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTATCGGAAAGG‑TTAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGc                                                     <  2:1065246/150‑1 (MQ=255)

                                                                                                 aGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATA‑ATTAAGGCAGGACTTTAGCGGGAAGG‑ATAAAAAACGGTTTTGGCGGCACATTGTGGTACGGACCAACGGCATGAGTGGGc                                                     >  2:219936/1‑150 (MQ=255)

                                                                                                      aaGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTg                                                <  1:1231885/150‑1 (MQ=255)

                                                                                                              aaTCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                        <  2:934593/150‑1 (MQ=255)

                                                                                                                  cGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                    <  1:304506/150‑1 (MQ=255)

                                                                                                                    ggTTACGCCGGGAGAGCAGAAGCG‑CGGG‑GAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                                  >  2:636102/1‑150 (MQ=255)

                                                                                                                    ggTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGATGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                                  <  1:1141981/150‑1 (MQ=255)

                                                                                                                    ggTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                                  <  1:796495/150‑1 (MQ=255)

                                                                                                                     gTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGc                                 <  1:528026/150‑1 (MQ=255)

                                                                                                                        aCGCCGGGAGAGCAGAAGCG‑CGGG‑CCAT‑ATTGGCCGCGCGAGAAATC‑ATTTAGGCAGGACGTTATCGGGAAGG‑ATAACTACCGGTTTTTACGGCACAATGGGGGAAGGACCCACGTCATGGGGGCGCCCCTGGGGAATCTTATAGGCAAc                              >  2:2159770/1‑150 (MQ=255)

                                                                                                                         cGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACg                             <  2:1417789/150‑1 (MQ=255)

                                                                                                                           ccGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTc                           >  2:1371264/1‑150 (MQ=255)

                                                                                                                            cGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTcc                          <  2:138481/150‑1 (MQ=255)

                                                                                                                                 gagCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGcc                     <  1:1088863/150‑1 (MQ=255)

                                                                                                                                  tgCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCt                    <  2:2100426/149‑1 (MQ=255)

                                                                                                                                     aGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc                 <  1:700411/150‑1 (MQ=255)

                                                                                                                                     aGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc                 <  2:443823/150‑1 (MQ=255)

                                                                                                                                      gAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca                <  1:1886557/150‑1 (MQ=255)

                                                                                                                                      gAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca                >  1:192864/1‑150 (MQ=255)

                                                                                                                                      ctcgcg‑cGGG‑CAAT‑CTTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca                >  2:1106850/4‑150 (MQ=255)

                                                                                                                                        cgcg‑cGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacc              >  1:2385683/2‑150 (MQ=255)

                                                                                                                                        aGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacc              <  2:2048223/150‑1 (MQ=255)

                                                                                                                                        aGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacc              >  2:893237/1‑150 (MQ=255)

                                                                                                                                         gcg‑cGGG‑CAAT‑CTTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacca             >  2:2143954/1‑150 (MQ=255)

                                                                                                                                          cg‑cgGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcaccac            <  1:2323377/150‑1 (MQ=255)

                                                                                                                                          cg‑cgGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcaccac            >  2:290158/1‑150 (MQ=255)

                                                                                                                                              ggg‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACttt         <  1:1348793/150‑1 (MQ=255)

                                                                                                                                               gg‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa        >  2:1519071/1‑150 (MQ=255)

                                                                                                                                                g‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc       <  1:636102/150‑1 (MQ=255)

                                                                                                                                                  cAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACTGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCTGCAACGTCCATGCCTTTCACCACTTTAcc      <  2:599138/150‑1 (MQ=255)

                                                                                                                                                   aaG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCa     <  1:1715386/150‑1 (MQ=255)

                                                                                                                                                    aG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaa    <  1:1243150/150‑1 (MQ=255)

                                                                                                                                                       aTTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAAt  <  1:1143539/150‑1 (MQ=255)

                                                                                                                                                       aTTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAAt  <  1:2188349/150‑1 (MQ=255)

GAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCG‑CGGG‑CAAG‑ATTGCCCGCGCGAGAAATC‑ATTACGGCAGGACTTTAGCGGAAAGG‑ATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAAT  >  NC_000913/3491557‑3491854

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