BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,980,112	C→A	5.0%	S120S (TCC→TCA) ‡	uspC →	universal stress protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,980,112	0	C	A	5.0%	358.0 / 7.0	120	S120S (TCC→TCA) ‡	uspC	universal stress protein
Reads supporting (aligned to +/- strand): ref base C (50/64); new base A (2/4); total (52/68)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 126 total aligned reads displayed.

AAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTT  >  NC_000913/1979963‑1980260

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

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

  gTTAATTCAGGATGCGGGTTATCACGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTTCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTg                                                                                                                                                    <  2:1413692/150‑1 (MQ=255)

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

     aaTTCAGGATTCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGTTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTc                                                                                                                                                 <  2:1538546/150‑1 (MQ=255)

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

      aTTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGATCATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCg                                                                                                                                                >  2:30971/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                                                                                                                                             >  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)

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

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

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

                          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                                                                                                                      <  1:1494279/150‑1 (MQ=255)

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

                                  aGACGTTTATTGCCTACGTTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                                    <  2:19392/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‑CGCGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCt                                                                                                          <  1:946747/150‑1 (MQ=255)

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

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

                                             gCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTg                                                                                                         >  1:770723/1‑150 (MQ=255)

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

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

                                               cTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGtt                                                                                                       <  1:390181/150‑1 (MQ=255)

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

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

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

                                                        gTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAATTCCACt                                                                                              >  1:284906/1‑150 (MQ=255)

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

                                                            aGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACg                                                                                          >  2:575341/1‑150 (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)

                                                                                         gCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTg                                                             >  2:743942/1‑150 (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)

                                                                                                    aTTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCa                                                  >  1:892370/1‑150 (MQ=255)

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

                                                                                                          tGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTa                                            <  2:191050/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)

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

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

                                                                                                                        cacaatcacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCt                              <  1:708844/150‑1 (MQ=255)

                                                                                                                              cacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACACAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgc                        >  1:1688797/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)

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

                                                                                                                                        tttCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGa              <  1:511868/150‑1 (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‑cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTa          >  2:176085/1‑150 (MQ=255)

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

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

                                                                                                                                              cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGAAGCGCGGaaagcc         <  1:1488242/150‑5 (MQ=255)

                                                                                                                                                   aTCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCTCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCt    <  2:675533/150‑1 (MQ=255)

                                                                                                                                                    tCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCtt   <  1:1393207/150‑1 (MQ=255)

                                                                                                                                                     ccTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCttt  <  2:530846/150‑1 (MQ=255)

                                                                                                                                                     ccTGCTCGGCGAAAAGAGTTAGTGCCTCAAGTGAGGGCGATGTGCTGTGAGTTCCACTTAAGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCttt  <  2:833868/150‑1 (MQ=255)

                                                                                                                                                     ccTGCTCGGCGAAAAGAGTTAATGCCTCAAGTGAGGTCGATGAGCTGTTAGTTCCACTTACGTGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTTCTGGCTTTCCGCTCTTCGCTGAACTATCTGCttt  >  2:750967/1‑150 (MQ=255)

AAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTT  >  NC_000913/1979963‑1980260

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