BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,980,110	T→G	5.4%	S120A (TCC→GCC) ‡	uspC →	universal stress protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,980,110	0	T	G	5.4%	358.0 / 7.3	122	S120A (TCC→GCC) ‡	uspC	universal stress protein
Reads supporting (aligned to +/- strand): ref base T (51/64); new base G (2/5); total (53/69)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 123 total aligned reads displayed.

TAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTT  >  NC_000913/1979962‑1980259

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

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

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

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

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

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

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

           aGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGa                                                                                                                                           >  1:796857/1‑150 (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)

                       aTCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTg                                                                                                                               <  1:1124252/150‑1 (MQ=255)

                         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)

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

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

                                          tATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACTAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGATCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtgt                                                                                                            >  1:158925/1‑150 (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)

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

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

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

                                                       gAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCa                                                                                               <  2:124488/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                 aTGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACg                                                                     <  1:1337686/150‑1 (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)

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

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

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

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

                                                                                                                 gCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTATGGAAAGGTAGCAACTTTATCgcgc                                     >  1:907123/1‑150 (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)

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

                                                                                                                        tcacaatcacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGGCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGc                              <  1:1190144/150‑1 (MQ=255)

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

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

                                                                                                                             atcacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGGCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTcc                         <  1:1470240/150‑1 (MQ=255)

                                                                                                                               cacaGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACACAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgc                       >  1:1688797/1‑150 (MQ=255)

                                                                                                                               caataactCTTTTCGCCGAGCAGGATGCTC‑GCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCcgc                       >  2:1104586/8‑150 (MQ=38)

                                                                                                                                 caGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCgcgc                     >  2:889686/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:511868/150‑1 (MQ=255)

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

                                                                                                                                          ttCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGaa            <  1:1184861/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:1302308/150‑4 (MQ=255)

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

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

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

TAAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCG‑CGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTT  >  NC_000913/1979962‑1980259

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