BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,980,103:1	+C	5.1%	coding (351/429 nt)	uspC →	universal stress protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,980,103	1	.	C	5.1%	489.9 / 18.8	118	coding (351/429 nt)	uspC	universal stress protein
Reads supporting (aligned to +/- strand): ref base . (52/60); new base C (2/4); total (54/64)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

Only 100 of 127 total aligned reads displayed.

TTTCTTGATAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTA  >  NC_000913/1979954‑1980252

tttCTTGATAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTcg                                                                                                                                                        <  2:1339623/150‑1 (MQ=255)

 ttCTTGATAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTcg‑c                                                                                                                                                      <  1:997378/150‑1 (MQ=255)

 ttCTTGATAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTcg‑c                                                                                                                                                      <  1:1522215/150‑1 (MQ=255)

   cTTGATAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGa                                                                                                                                                    >  2:841112/1‑150 (MQ=255)

    ttGATAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTcg‑ccga                                                                                                                                                   <  1:1589132/150‑4 (MQ=255)

    ttGATAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAg                                                                                                                                                   <  1:160675/150‑1 (MQ=255)

        tAAGTTAAATCAGGATGCGGGGTATCACGTTCACAAGACGTATATTGCCCACGGTCAGTTAAGCGAACATATACTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTAG‑CGAGCATc                                                                                                                                               >  1:491296/1‑150 (MQ=255)

         aaGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATcc                                                                                                                                              <  1:470917/150‑1 (MQ=255)

          aGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCt                                                                                                                                             >  2:1184861/1‑150 (MQ=255)

            ttAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGc                                                                                                                                           <  1:841112/150‑1 (MQ=255)

              aaTTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACTGTTTCTTTTCG‑CGAGCATCCTGCTc                                                                                                                                         >  1:1692133/1‑150 (MQ=255)

               aTTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCg                                                                                                                                        <  1:1698435/150‑1 (MQ=255)

                 tCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGc                                                                                                                                      >  1:375035/1‑150 (MQ=255)

                  cAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCg                                                                                                                                     <  2:178335/150‑1 (MQ=255)

                  cAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCg                                                                                                                                     <  2:649579/150‑1 (MQ=255)

                  cAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCg                                                                                                                                     >  1:110917/1‑150 (MQ=255)

                  cAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTAGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGct                                                                                                                                     >  2:72787/1‑149 (MQ=255)

                  aaGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCTGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCggag                                                                                                                                     >  1:1607531/2‑148 (MQ=255)

                   aGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGa                                                                                                                                    >  1:796857/1‑150 (MQ=255)

                    ggATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGaa                                                                                                                                   <  2:884083/150‑1 (MQ=255)

                     gATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGaaa                                                                                                                                  <  2:1159159/150‑1 (MQ=255)

                      aTGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTC‑GCGaaaa                                                                                                                                 <  1:1395495/150‑1 (MQ=255)

                       tGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAag                                                                                                                                <  1:480261/150‑1 (MQ=255)

                          ggggTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGt                                                                                                                             >  2:1582102/1‑150 (MQ=255)

                           gggTATCCCGTTGACAAGCCGTTTATTGCCTACGGTGAGTTACGCGAACAAATTCTGGAAGTATGTCACAAGCATCAATTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTC‑GCGAAAagaaac                                                                                                                            <  1:716255/150‑4 (MQ=255)

                           gggTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTC‑GCGAAAagaaac                                                                                                                            <  2:61526/150‑4 (MQ=255)

                            ggTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATATCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCCGAGCAGGATGCTC‑GCGAAAagaaact                                                                                                                           <  1:258599/150‑5 (MQ=38)

                                 cccGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGcc                                                                                                                      >  2:529746/1‑150 (MQ=255)

                                   cGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTc                                                                                                                    <  2:563380/150‑1 (MQ=255)

                                    gTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCa                                                                                                                   <  2:1046724/150‑1 (MQ=255)

                                     ttGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCaa                                                                                                                  >  1:132587/1‑150 (MQ=255)

                                      tGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAg                                                                                                                 <  2:1170495/150‑1 (MQ=255)

                                       gACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGt                                                                                                                >  1:109384/1‑150 (MQ=255)

                                        aCAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAATAGTTATTGCCTCAAGTg                                                                                                               >  2:863940/1‑150 (MQ=255)

                                         cAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGa                                                                                                              >  2:511868/1‑150 (MQ=255)

                                         cAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGa                                                                                                              <  1:1494279/150‑1 (MQ=255)

                                           aGACGTTTATTGCCTACGTTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                            <  2:19392/150‑1 (MQ=255)

                                           aGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTGTGGAAGTATGTCACAAGCATCATATCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGGGAgg                                                                                                            <  2:1460468/150‑1 (MQ=255)

                                            gACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCT‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGACTCAAGTGAGGt                                                                                                           >  1:184994/1‑150 (MQ=255)

                                            gACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGt                                                                                                           <  2:110917/150‑1 (MQ=255)

                                            gACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGt                                                                                                           <  2:1551316/150‑1 (MQ=255)

                                             aCGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGGCACAAGCATCATTTCTATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTc                                                                                                          <  1:878174/150‑1 (MQ=255)

                                              cGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCg                                                                                                         >  1:124488/1‑150 (MQ=255)

                                                  tATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACTAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGATCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtgt                                                                                                     >  1:158925/1‑150 (MQ=255)

                                                    ttGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtgtgg                                                                                                   <  1:465681/150‑2 (MQ=255)

                                                     tGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCt                                                                                                  <  1:725247/150‑1 (MQ=255)

                                                      gCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTg                                                                                                 <  2:1507433/150‑1 (MQ=255)

                                                       ccTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAGTCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGt                                                                                                >  1:716343/1‑150 (MQ=255)

                                                       ccTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGt                                                                                                >  1:974588/1‑150 (MQ=255)

                                                         tACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTa                                                                                              <  1:933138/150‑1 (MQ=255)

                                                             gTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTc                                                                                          >  2:1393207/1‑150 (MQ=255)

                                                             gTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTc                                                                                          <  1:1291937/150‑1 (MQ=255)

                                                              tGAGTTAAGCGAACATATTGTGGAAGTATGTCACAAGCATAATTTCGATTTGGTGATTTGCGGTAAGCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTcc                                                                                         <  1:881283/150‑1 (MQ=255)

                                                                aGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCAc                                                                                       >  1:219755/1‑150 (MQ=255)

                                                                  ttAAGCGAACATATTCTGGAATTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACtt                                                                                     <  1:139184/150‑1 (MQ=255)

                                                                  ttAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACtt                                                                                     <  2:1304582/150‑1 (MQ=255)

                                                                   tAAGCGAACATATTCTGGAAGTATCTCACAAGCATCATTTCGATTTGGTGATTTGCGTTAATCACAATCACAGTTTCTTTTCG‑CGAGCTTCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTa                                                                                    <  1:297494/150‑1 (MQ=255)

                                                                     aGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACg                                                                                  >  1:285370/1‑150 (MQ=255)

                                                                     aGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATGTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACg                                                                                  <  1:54692/150‑1 (MQ=255)

                                                                      gCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACgg                                                                                 >  1:90843/1‑150 (MQ=255)

                                                                      gCGAACATATTCTGGAAGTATGACACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACgg                                                                                 >  1:1051608/1‑150 (MQ=255)

                                                                         aaCATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGgag                                                                              <  2:162891/150‑1 (MQ=255)

                                                                          aCATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGgaga                                                                             <  1:107663/150‑1 (MQ=255)

                                                                           cATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGAt                                                                            <  2:491699/150‑1 (MQ=255)

                                                                              atTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTaa                                                                         >  2:128816/1‑150 (MQ=255)

                                                                                tCTGGAAGTATGTCATAAGCATCATTTCGATTTGGTTATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CTAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                       <  2:491296/150‑1 (MQ=255)

                                                                                tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTTGTTGCACTTACGGGAGATTAAcc                                                                       >  2:300893/1‑150 (MQ=255)

                                                                                tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                       >  1:12658/1‑150 (MQ=255)

                                                                                tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                       <  1:60555/150‑1 (MQ=255)

                                                                                   ggAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCt                                                                    <  2:1438194/150‑1 (MQ=255)

                                                                                        tATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTAc                                                               >  2:915293/1‑150 (MQ=255)

                                                                                        tATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTAc                                                               >  1:214032/1‑150 (MQ=255)

                                                                                             cacaAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAg                                                          >  1:675533/1‑150 (MQ=255)

                                                                                                  gCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTg                                                     >  2:743942/1‑150 (MQ=255)

                                                                                                    atcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTAGGCGAAAAGAGTTATTGCCTAAAGTGAGGTCGATGTGATGTTAGTTCCACTTACTGGAGATTAACCGATACTACGAAATCTTTgta                                                   >  1:1698251/1‑148 (MQ=255)

                                                                                                      catTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGaaa                                                 >  2:1488242/1‑150 (MQ=255)

                                                                                                         ttCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGt                                              >  2:1316860/1‑150 (MQ=255)

                                                                                                          tCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTa                                             <  2:954175/150‑1 (MQ=255)

                                                                                                            gATTTGGTGATTTGAGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGc                                           <  2:865669/150‑1 (MQ=255)

                                                                                                                tGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACt                                       >  2:451706/1‑150 (MQ=255)

                                                                                                                 ggTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACtt                                      <  2:375035/150‑1 (MQ=255)

                                                                                                                   tGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTa                                    <  2:285370/150‑1 (MQ=255)

                                                                                                                         gCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTATGGAAAGGTAGCAACTTTATCgcgc                              >  1:907123/1‑150 (MQ=255)

                                                                                                                          cGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgct                             <  2:1271960/150‑1 (MQ=255)

                                                                                                                           ggTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAGCCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctg                            >  2:185616/1‑150 (MQ=255)

                                                                                                                               atcacaatcacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTACTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTgg                        <  2:284906/150‑1 (MQ=255)

                                                                                                                                tcacaatcacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGc                       <  1:1337838/150‑1 (MQ=255)

                                                                                                                                   caatcacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGAGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCttt                    >  1:738805/1‑150 (MQ=255)

                                                                                                                                         caGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGGGAGGTCGATGTGCTGTTAGTTCCACTTACGGGTGATTAACCGCTCCTACGCATGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCgcgc              <  2:1051608/150‑1 (MQ=255)

                                                                                                                                          aGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCt             >  1:1527327/1‑150 (MQ=255)

                                                                                                                                              tCTTTTCGCCGAGCAGGATGCTC‑GCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcgc         >  1:1137985/1‑150 (MQ=255)

                                                                                                                                               cTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGACTTTCCGCGCTTcgcg        >  2:934575/1‑150 (MQ=255)

                                                                                                                                               attttCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTAACGCGCTGCTGGCTTTCCGCGCTTcgcg        <  2:184994/149‑1 (MQ=255)

                                                                                                                                                 tttCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGa      >  1:1522448/1‑150 (MQ=255)

                                                                                                                                                  ttCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGaa     <  1:1184861/150‑1 (MQ=255)

                                                                                                                                                   tCG‑CGAGCATCCTTCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAAc    <  2:1692133/150‑1 (MQ=255)

                                                                                                                                                   tCG‑CGAGCAGCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAGGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAAc    <  1:1131696/150‑1 (MQ=255)

                                                                                                                                                   acg‑cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGCGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAAc    <  1:1137477/149‑1 (MQ=255)

                                                                                                                                                    cg‑cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGAAGCGCGGaaag   <  1:30971/150‑3 (MQ=255)

                                                                                                                                                     g‑cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGAAGCGCGGaaagc  <  2:1295352/150‑4 (MQ=255)

TTTCTTGATAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTA  >  NC_000913/1979954‑1980252

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