BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,229,177	A→C	12.5%	intergenic (+112/+261)	yohP → / ← dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,229,177	0	A	C	12.5%	606.3 / 12.9	295	intergenic (+112/+261)	yohP/dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand): ref base A (64/164); new base C (27/11); total (101/194)
Fisher's exact test for biased strand distribution p-value = 6.18e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326

   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAggg                                                                                                                                                     >  2:9768145/1‑150 (MQ=255)

      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GCgggggg                                                                                                                                                  >  2:6856558/1‑149 (MQ=255)

      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggc                                                                                                                                                  >  2:3613672/1‑150 (MQ=255)

       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcg                                                                                                                                                 <  2:2321307/150‑1 (MQ=255)

        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcgg                                                                                                                                                >  2:2475492/1‑150 (MQ=255)

        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcgg                                                                                                                                                >  2:3479821/1‑150 (MQ=255)

        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcgg                                                                                                                                                <  1:2537865/150‑1 (MQ=255)

        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAAGGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAACCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGCAGCGC‑G‑CAggggggcg                                                                                                                                                >  2:1694881/1‑147 (MQ=255)

         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcggc                                                                                                                                               <  1:536918/150‑1 (MQ=255)

           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGAGC‑G‑CAGGGGGCGGCCa                                                                                                                                             <  1:6807285/150‑1 (MQ=255)

            cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCaa                                                                                                                                            <  2:10519040/150‑1 (MQ=255)

            cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCaa                                                                                                                                            <  1:189761/150‑1 (MQ=255)

              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATc                                                                                                                                          <  2:8566334/150‑1 (MQ=255)

               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GCGGGGGCGGccccacc                                                                                                                                         >  2:4971205/1‑145 (MQ=255)

                cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcgc                                                                                                                                        <  2:1828897/150‑1 (MQ=255)

                 gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcgcc                                                                                                                                       <  1:7690220/150‑1 (MQ=255)

                 gTAATTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGAGC‑G‑CAGGGGGCGGCCAATcgcc                                                                                                                                       <  2:6811032/150‑1 (MQ=255)

                  tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcgccg                                                                                                                                      <  1:10573468/150‑1 (MQ=255)

                   aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcgccgc                                                                                                                                     <  2:4512716/150‑1 (MQ=255)

                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGGATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCGATTGgccgcc‑ccct                                                                                                                               <  1:7222049/150‑2 (MQ=39)

                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCGGGGGGGGCCAAACcccccc‑cccc                                                                                                                               >  2:2161307/1‑150 (MQ=255)

                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGAGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑Gcccc                                                                                                                              <  2:6649383/150‑1 (MQ=255)

                         gatgatATATTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑Gcccc                                                                                                                              <  1:5819501/150‑1 (MQ=255)

                          atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑Gccccc                                                                                                                             <  1:9029201/150‑1 (MQ=255)

                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGGGGGCAAAAGCcccc‑ccccccc                                                                                                                            >  1:10662899/1‑149 (MQ=255)

                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAACcccccc‑ccccccag                                                                                                                           >  1:10069770/1‑148 (MQ=255)

                            gatATTTTAAAAATAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGACTCTCCTCCTGTTGCGCTCCTGTCATCTGCATTGCTTCATATAAACACTGGCGCAAGGAGCGC‑G‑CACGGGGCGCCCAAAACCCGcc‑ccccccag                                                                                                                           >  2:6972315/1‑148 (MQ=255)

                              tatTGTAAAATTAATTAATGTCCTCAGGTCCGCACATAACGAGAATATGTCAGTCTCACATCCCGTTGCGCTCCTGTCAGGTGCATTGCTTCATATAATCACCGCCGCAAGGAGCGC‑G‑CAGGGGCCGGCCAATCGCCGCC‑GCCCCCTGCa                                                                                                                         <  2:299375/150‑1 (MQ=37)

                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCGTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCAccc                                                                                                                      <  2:2839091/150‑1 (MQ=255)

                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CTGGGGGCGGCCAAAAGCCCCC‑CCCCCCCCCAccc                                                                                                                      >  1:8141772/1‑150 (MQ=38)

                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAACACCCCC‑CCCCCCAGCAccc                                                                                                                      >  1:2749005/1‑150 (MQ=38)

                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCAccccc                                                                                                                    <  2:8476987/150‑1 (MQ=255)

                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCAccccc                                                                                                                    <  2:1804315/150‑1 (MQ=255)

                                      aattaattACGGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCggg                                                                                                                 <  2:7607697/150‑1 (MQ=255)

                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CGGGGGGGGGCCAAAC‑CCCCC‑CCCCCCCCCAcccccccg                                                                                                                 >  1:9125281/1‑147 (MQ=255)

                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CGCGGGGCGGGCCCTCGCCCCC‑CCCCCCCGCccccccccg                                                                                                                 >  2:9917656/1‑147 (MQ=38)

                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCCAAAGCCGCC‑CCCCCCACCAcccccccg                                                                                                                 >  2:8856406/1‑147 (MQ=38)

                                          aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGctct                                                                                                             <  2:5725440/150‑1 (MQ=255)

                                           attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCCGGGGCGGCCCAACCCCCCC‑GCCCCCCGCACCCCCCGGCGCTg                                                                                                            >  2:6792807/1‑150 (MQ=255)

                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGc                                                                                                          <  2:1082978/150‑1 (MQ=255)

                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CGGGGGGGGGCCCCACGCCCCC‑CCCCCCACCACCCCCCGGCTCAGGCgca                                                                                                       >  1:10903781/1‑148 (MQ=255)

                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACa                                                                                                     <  2:7249515/150‑1 (MQ=255)

                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAACGCCCCC‑GCCCCCTGCACCCCCCGGCTCAGGCGAACa                                                                                                     >  2:3702132/1‑150 (MQ=255)

                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCCATCCCCGCC‑CCCCCCAGCACCCCCCGGCTCAGGCGAACaaaa                                                                                                  >  2:5494726/1‑150 (MQ=25)

                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACaaaa                                                                                                  <  2:7172592/150‑1 (MQ=255)

                                                         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                              <  2:10875029/150‑1 (MQ=255)

                                                             gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCAGCCCCCCCGGGCGCTGGCGAACAAAAAcgccgcg                                                                                          >  2:6562675/1‑149 (MQ=255)

                                                                 aTAACTATAATATTTCATTCTCTAATCCTTTTGCGCTCCTTTATTTTGAATTGATTAATATAATCAGTGGCGAATGGTGCGC‑G‑ATGGGGGCGGCCATTCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCTAACAAAATCGCCGCTgcgc                                                                                      <  2:1918275/150‑1 (MQ=12)

                                                                    aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CGGGGGGCGGCCCAAACCCCCC‑CCCCCCACCACCCCCCCGCGCAGGCGCACAAAAAAGCCCCGGCGCCGt                                                                                   >  2:1591411/1‑150 (MQ=255)

                                                                      gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCCAACGCCCCC‑CCCCCCAGCCCCCCCCGGCACCGGCCCACAAAAACCCCCCGGCGCGGTGc                                                                                 >  2:7135809/1‑150 (MQ=21)

                                                                         aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAACGCCCCC‑CCCCCCGGCACCCCCCGGCACAGGCGAACAAAATCGCCGCGGCGCGGGGccca                                                                              >  2:1635336/1‑149 (MQ=21)

                                                                          atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCCCTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCTCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:561425/150‑1 (MQ=255)

                                                                          atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:642116/150‑1 (MQ=255)

                                                                           tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCCGGGGGGGCCAAAAGCCCCC‑CCCCCCCCCCCCCCCCGGCGCAGGGGGACAAAAAAGCCCCAGCGCCGGGCCcgcg                                                                            >  2:4367319/1‑150 (MQ=255)

                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGAAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:5271917/150‑1 (MQ=255)

                                                                             tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                          <  2:8319083/150‑1 (MQ=255)

                                                                               tCAGTCTCTCATCCTTTTGGGCTCCTTTCTTGTGCTTTGCTACATATAATCCCGGGCGCATGGCGCGC‑G‑CTGGGGGCGGCCATTCGCCGCC‑GCCCCGTGCACCCCCGGGCTCTGGCTAACAAAATCGCCGCTGCGCGGTGCCCTCGGCtt                                                                        <  1:837701/150‑1 (MQ=17)

                                                                                cAGTCTCTCATCGTGTTGCGCTCCTGTCTTGTGCTTTGCTTCATATAATCAGGGGCGAAAGGAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                       <  1:9917656/150‑1 (MQ=17)

                                                                                   tctctcATCCTGTTGCGCTCCTGTCAGGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  2:1934874/150‑1 (MQ=255)

                                                                                    ctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAACGCCGCC‑GCCCCCAGCCCCCCCCGGCTCAGGCGAACAAAATCCCCGCTGCGCGGTGCCCACCGCGTATccc                                                                   >  2:6397576/1‑150 (MQ=18)

                                                                                     tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCCAAACCCGCC‑CCCCCCAGCACCCCCCGGCTCAGGCGAAAAAAAACGCCCCAGCGCGGTGCCCCACGCGTATccca                                                                  >  2:2618373/1‑149 (MQ=18)

                                                                                       tcATCCTGTTGCGCTCCTGTCTTGTGAATTGATTCATATAATCACTGGCGCACGGTGCGC‑G‑CTGGGGGGGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:9297135/150‑1 (MQ=18)

                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCCAAAGCCCCC‑CCCCCCACACCCCCCCGGCACTGGCGAAAAAAAAAGCCCCCGCGCGGGGCCCTCGGCATATACCATa                                                                >  2:6003817/1‑150 (MQ=16)

                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:9225492/150‑1 (MQ=255)

                                                                                         aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  1:7809283/150‑1 (MQ=255)

                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGGGGCCCATCGCCCCC‑CCCCCCCGCCCCCCCCGGCGCGGGGGGACAAAATCCCCCCCGCGCGGGGCCCCCGGCGTATACCATACGGCGAcccg                                                      >  1:492530/1‑148 (MQ=12)

                                                                                                  cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCCGGGGCGGCCCAAAGCCCCC‑CCCCCCAGCACCCCCCGGCGCAGGCGCACAAAAACGCCCCGGCGCGGTGCCCTCCGCGTCTCCCTGAAGGCTACCcgg                                                     >  2:10071802/1‑150 (MQ=12)

                                                                                                    cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  1:427120/150‑1 (MQ=255)

                                                                                                    cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CATGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  1:3861370/150‑1 (MQ=255)

                                                                                                    cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  2:5897175/150‑1 (MQ=255)

                                                                                                    cTCCTGTCATGTGCATTGCTTCATATAACCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  1:1168809/150‑1 (MQ=255)

                                                                                                     tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GACGGGGCCGGCAATCGCCCCC‑CCCCCCGGCCCCCCCCGGCTCGGGCGCACAAAAAAGCCCGGGCGCGGGGCCCGCGGCGTATCCCCTACCGCGACCCGGGCg                                                  >  2:8906378/1‑150 (MQ=12)

                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:3467551/150‑1 (MQ=255)

                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:2087289/150‑1 (MQ=255)

                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  1:1313537/150‑1 (MQ=255)

                                                                                                               tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                        <  2:5949237/150‑1 (MQ=255)

                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:9351181/150‑1 (MQ=255)

                                                                                                                 aaTTGCTTAATATAATCACTGGCGAAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCTCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  1:1436800/149‑1 (MQ=255)

                                                                                                                  aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAAc                                     <  2:8103444/150‑1 (MQ=255)

                                                                                                                  aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAAc                                     <  2:10662899/150‑1 (MQ=255)

                                                                                                                  aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAAc                                     <  2:9028900/150‑1 (MQ=255)

                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  1:728930/150‑1 (MQ=40)

                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCCCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  2:9557941/150‑1 (MQ=255)

                                                                                                                      cTTCATATAATCACGGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  1:3999427/150‑1 (MQ=255)

                                                                                                                       gtCACATCATCACTGGCGCCCGGACCGC‑G‑CACGGCGCGGCCCCCCGCCGCC‑GCCCCCAGCACCCCCGGGCTCGGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCATACGGCGACCGGGTCGGGCGCGCGGAAACAAccc                                <  2:7291084/149‑1 (MQ=11)

                                                                                                                        tCATATAATCACTGGCGCACGGAGCGC‑G‑CCGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:4636132/150‑1 (MQ=255)

                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:1967155/150‑1 (MQ=39)

                                                                                                                         cATTTAATCACTGGGGCAAGGTTGGG‑G‑GAGGGGGGGGCCAATCGCCGCCGGCCCCCGGAACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:10268070/150‑1 (MQ=11)

                                                                                                                            ataATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaa                           <  1:7559046/150‑1 (MQ=35)

                                                                                                                              aaTCACTGGCGCACGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:8817566/150‑1 (MQ=21)

                                                                                                                               aTCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCAATAGCCCCC‑GCCCCCGGCCCCCCCGGGCGCTGGCGAACAAAATCGCCGCGGCGCGGGGCCCGCGGCTGATCCCGTACCGCTACCCGGTCGGGCGCGCGGGAACAACCCAGTAAAAc                        >  2:9105715/1‑150 (MQ=11)

                                                                                                                               aTCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  2:7806729/150‑1 (MQ=33)

                                                                                                                               aTCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:1822479/150‑1 (MQ=33)

                                                                                                                                 cACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:3107801/150‑1 (MQ=32)

                                                                                                                                   cTGGCGCATGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:3724869/150‑1 (MQ=17)

                                                                                                                                          aaGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  1:3520378/150‑1 (MQ=32)

                                                                                                                                              aGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac         <  2:1389625/150‑1 (MQ=32)

                                                                                                                                               gcgc‑g‑cAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca        <  1:6249011/150‑1 (MQ=32)

                                                                                                                                                      cAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  2:5120468/150‑1 (MQ=32)

                                                                                                                                                      cAGGGGGAGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCTCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  2:383061/150‑1 (MQ=17)

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCC‑GCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326

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