BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	644,563	G→C	5.5%	C147W (TGC→TGG)	rna ←	ribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	644,563	0	G	C	5.5%	323.3 / 7.6	101	C147W (TGC→TGG)	rna	ribonuclease I
Reads supporting (aligned to +/- strand): ref base G (43/52); new base C (2/3); total (46/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

Only 100 of 109 total aligned reads displayed.

CGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAGCGCATC  >  NC_000913/644414‑644710

cggcTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCcagg                                                                                                                                                     >  2:1398967/1‑150 (MQ=255)

cggcTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCcagg                                                                                                                                                     <  1:1787088/150‑1 (MQ=255)

 ggcTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCcaggc                                                                                                                                                    <  1:179882/150‑1 (MQ=255)

    ccACGGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCcaggcagg                                                                                                                                                 <  2:1725251/149‑1 (MQ=255)

           tttACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCTGCCTGGAACgcta                                                                                                                                          <  1:1103882/150‑4 (MQ=255)

            ttACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAAACAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCg                                                                                                                                         <  2:573527/150‑1 (MQ=255)

            atcccttAATTCTCCGCAAGGAATTTTCCGGCTTCGCTCGCTTTAATTTCTTGAATCAGGCGCACCACCCTACCGTAGTATGAGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTCGCGTTCCAGGCAGGAACGTCCg                                                                                                                                         <  2:271333/144‑1 (MQ=255)

             tACCGTAATTATCCGCCAGGAAATTACCGGCTTCGCTCTCTATAATTTCTTGAGTCAGGCGTACCATCGTCCCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCTGCCTGGAACgctacg                                                                                                                                        <  2:270730/150‑6 (MQ=255)

              aCCGTAATTATCCGCAAGGAATTTTCCGGCTTGGCTCTCTTTAATTTATTTATTCAGGCTTACCATCGTACCGAAGTATCCGTCCGGATCAAAACAAAAGCAGGCACCTTGTTTGGCAGATTCGTAGCGTTCAAGACAGCAACTTCCGcc                                                                                                                                       >  1:243517/1‑150 (MQ=255)

              aCCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCTGCCTGGAACgctacga                                                                                                                                       <  1:867147/150‑7 (MQ=255)

              aCCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGcc                                                                                                                                       <  1:1267619/150‑1 (MQ=255)

                 gTAATTATCCTCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGc                                                                                                                                    <  1:709109/150‑1 (MQ=255)

                 gTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATATCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGc                                                                                                                                    <  1:971620/150‑1 (MQ=255)

                        tCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGc                                                                                                                             >  1:1284814/1‑150 (MQ=255)

                          cGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCAt                                                                                                                           >  1:1708572/1‑150 (MQ=255)

                          cGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGAATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCAt                                                                                                                           >  2:601516/1‑150 (MQ=255)

                            cAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGa                                                                                                                         >  1:558095/1‑150 (MQ=255)

                             aaGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGAc                                                                                                                        <  2:1702233/150‑1 (MQ=255)

                             aaGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGAc                                                                                                                        <  1:782423/150‑1 (MQ=255)

                             aaGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGAc                                                                                                                        <  1:183717/150‑1 (MQ=255)

                             aaGGAATTTTCCGGCGTCGCTGTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGCAGTATGCGTCCGGATCAAAACCACAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATgtc                                                                                                                        >  1:893584/1‑148 (MQ=255)

                               ggAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACtt                                                                                                                      >  2:812824/1‑150 (MQ=255)

                                  aTTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCAc                                                                                                                   <  2:800644/150‑1 (MQ=255)

                                  aTTTCCCGGCTTCGCTCTCATTAATATCTTGATTCAGGAGTACCATCGTACCGAAGTCTGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATCTTCGTAGCGTTCCAGGCCGGAACGTCCGCCAGCTCCTGGCATGACTTCAc                                                                                                                   <  2:180781/150‑1 (MQ=255)

                                    tttCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTgg                                                                                                                              <  1:87564/137‑1 (MQ=255)

                                    tttCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTgg                                                                                                                              >  2:87564/1‑137 (MQ=255)

                                    tttCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACtt                                                                                                                 >  1:648932/1‑150 (MQ=255)

                                    tttCCGGCTTCGCTCTCTTTAATATCTTGATTCAGGCGTACCATCGTACCGTATTATGCGTCCGGATCAAAACCAAAGAAGGCACCATGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGACCGCCAGCTCCTGGCATGACTTCACtt                                                                                                                 >  1:511281/1‑150 (MQ=255)

                                       ccGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCCTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGt                                                                                                              <  2:1512187/150‑1 (MQ=255)

                                        cGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGtt                                                                                                             <  2:1831945/150‑1 (MQ=255)

                                        cGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGAGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGtt                                                                                                             >  1:832416/1‑150 (MQ=255)

                                            ttCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGc                                                                                                         <  1:1326759/150‑1 (MQ=255)

                                            ttCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGc                                                                                                         <  1:689635/150‑1 (MQ=255)

                                            tgcGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCCGGCACCGTGGTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGc                                                                                                         <  1:1349951/148‑1 (MQ=255)

                                             tCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCg                                                                                                        <  2:75345/150‑1 (MQ=255)

                                               gCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGc                                                                                                      >  1:1321673/1‑150 (MQ=255)

                                                 tctctTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCAGGCATGACTTCACTTAGTTTAGCGGCCg                                                                                                    >  1:1351612/1‑150 (MQ=255)

                                                       tAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCa                                                                                              >  1:1313691/1‑150 (MQ=255)

                                                        aaTTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATTTTCGTAGCGTTCCGGGCAGGAACGTCCGCGAGCTGCTGGCATGAGTTCACTTAGGTTAGCGGCCGTTTCCAg                                                                                             >  2:599013/1‑150 (MQ=255)

                                                        aaTTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAg                                                                                             <  1:1252603/150‑1 (MQ=255)

                                                          tttCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTg                                                                                           <  1:1340459/150‑1 (MQ=255)

                                                           ttCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGa                                                                                          >  2:781563/1‑150 (MQ=255)

                                                            tCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGAt                                                                                         >  2:728231/1‑150 (MQ=255)

                                                             attGATTCAGGCTTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATa                                                                                        <  2:1008613/149‑1 (MQ=255)

                                                                  ttCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCg                                                                                   <  2:623283/150‑1 (MQ=255)

                                                                   tCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCgg                                                                                  <  2:1210621/150‑1 (MQ=255)

                                                                    cAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGt                                                                                 >  2:1250322/1‑150 (MQ=255)

                                                                      ggCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGttt                                                                               >  2:1670610/1‑150 (MQ=255)

                                                                      ggCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGttt                                                                               >  1:1359638/1‑150 (MQ=255)

                                                                               tCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGACGATg                                                                      <  1:603734/150‑1 (MQ=255)

                                                                               tCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGACGATg                                                                      <  1:597471/150‑1 (MQ=255)

                                                                                  tACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAac                                                                   >  1:1425228/1‑150 (MQ=255)

                                                                                  tACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAac                                                                   <  2:1389049/150‑1 (MQ=255)

                                                                                   aCCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGTCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAaca                                                                  <  2:799973/150‑1 (MQ=255)

                                                                                   aCCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAAGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAaca                                                                  <  1:829857/150‑1 (MQ=255)

                                                                                      gAAGTATGCGTCCGGATCAACACCAAAGCAGGCACCGTGTTTGGCATATTCGTATCGTTCCAGGCAGGAACGTCCGCCAGCTACTGGCATGACTTCACTTAGTTTAGCGGCCGTTGCCAGTGATAATCCGGTTTCCGGCGATGAACACAt                                                               >  1:1225395/1‑150 (MQ=255)

                                                                                           aTGCGTCCGGATCATAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGCTACTCCGGTTTCCGGCGATGAACACATTCGGc                                                          >  1:221727/1‑150 (MQ=255)

                                                                                            tGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCTTAGCGTTCCAGGCAGGAACGGCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCt                                                         <  1:818074/150‑1 (MQ=255)

                                                                                            tGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCt                                                         >  2:14005/1‑150 (MQ=255)

                                                                                              cGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTgg                                                       >  1:1044060/1‑150 (MQ=255)

                                                                                              cGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGGCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTgg                                                       <  1:215946/150‑1 (MQ=255)

                                                                                                 ccGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCGAGGCAGGAACGTCCGCCCGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGgcgc                                                    >  2:644525/1‑150 (MQ=255)

                                                                                                 ccGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGgcgc                                                    >  2:259586/1‑150 (MQ=255)

                                                                                                    gATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCAAGGCAGTAACTTACGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGGTTCCAGTCATAATACGGTTTCCGGCGATGAACAGATTCGGCTGgcgcgcg                                                 >  2:155311/1‑150 (MQ=255)

                                                                                                      tCAAAACCACAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCt                                               <  1:730/150‑1 (MQ=255)

                                                                                                       cAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCtt                                              >  1:978774/1‑150 (MQ=255)

                                                                                                         aaaCCAACGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCAGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCt                                            <  1:1170363/150‑1 (MQ=255)

                                                                                                         aaaCCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCt                                            >  1:1473218/1‑150 (MQ=255)

                                                                                                          aaCCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGCTAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTg                                           >  1:128733/1‑150 (MQ=255)

                                                                                                           aCCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTgg                                          >  2:123436/1‑150 (MQ=255)

                                                                                                             cAAAGCAGGCACCGTGTTGGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTa                                        <  2:897321/150‑1 (MQ=255)

                                                                                                                 gCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCCTGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATACTCCGGTTTCCGGCGATGAACACCTTCGGCTGGCGCGCGCTTCTGGTAGAtt                                    >  1:1144543/1‑150 (MQ=255)

                                                                                                                  cAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTc                                   >  1:422656/1‑150 (MQ=255)

                                                                                                                     gCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGCCTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCggg                                <  2:890550/150‑1 (MQ=255)

                                                                                                                     gCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCggg                                <  2:1810239/150‑1 (MQ=255)

                                                                                                                     gCACCGTGGTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGGTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCggg                                <  1:1032021/150‑1 (MQ=255)

                                                                                                                      cACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGa                               >  1:514675/1‑150 (MQ=255)

                                                                                                                         cGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTg                            <  1:812824/150‑1 (MQ=255)

                                                                                                                          gtgtTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTgg                           <  1:17894/150‑1 (MQ=255)

                                                                                                                           tgtTTGGCATCTTCGTAGCGTGCCAGGCCGGAACGTCCGCCAGCTCCCGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTggg                          <  2:1473784/150‑1 (MQ=255)

                                                                                                                           tgtTTCGCATATTCGTCGCTTTCCAGGCAGGAACGGCCGCCAGCTCCTTGCATGACTTCTCTTTGTTTAGCTGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAAATCATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTgtg                          <  2:511281/150‑3 (MQ=255)

                                                                                                                            ggttGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGGTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGc                         <  1:1398967/148‑1 (MQ=255)

                                                                                                                              ttGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGa                       <  1:601516/150‑1 (MQ=255)

                                                                                                                               tGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATACGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAg                      >  2:1644071/1‑150 (MQ=255)

                                                                                                                                 gCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTAAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTa                    >  1:905686/1‑150 (MQ=255)

                                                                                                                                  cATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCGGTTTCCAGGGATATTCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTGCTGGTAGTTTCGGTATTGGGCGAGTAg                   >  1:1208437/1‑150 (MQ=255)

                                                                                                                                  cATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCGTCTGGTAGATTCGGGATTGGGCGAGTAg                   >  2:230382/1‑150 (MQ=255)

                                                                                                                                        gcggacGTTCCTGCCTGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAAc             >  2:1799341/6‑150 (MQ=255)

                                                                                                                                         gTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAAcc            >  1:45066/1‑150 (MQ=255)

                                                                                                                                         gTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGAAGCGCAAcc            <  2:1624232/150‑1 (MQ=255)

                                                                                                                                         gTAGCGGACCAGGCAGGAACGGCCGCCAGCTCCTGGCATGACTTCACTTTGTATCGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACAGATTCGGCTGGCGCGCTCTTCTGGTAGAATCGGGATTGGGCGAGTAGCGCAAcc            <  2:1345051/150‑1 (MQ=255)

                                                                                                                                          tAGCGTTCCAGGCAGGAACGGCCGCCAGCTCCTGGCATGACGTCACTTCGTTTAGCGGCCGGTTCCAGTGAGAATCCGGTTGCCGGCGATGAACACATTCGGCTGGTGCGCGCTTCTGGTAGATTCGCGATTGGGCGAGTAGCGCAACCg           <  1:305050/150‑1 (MQ=255)

                                                                                                                                           gacGTTCCTGCCTGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGa          >  1:944243/3‑150 (MQ=255)

                                                                                                                                           aGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGa          >  2:1869686/1‑150 (MQ=255)

                                                                                                                                            gCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGaa         >  1:1428377/1‑150 (MQ=255)

                                                                                                                                             cGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCAATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAg        <  2:48776/150‑1 (MQ=255)

                                                                                                                                               ttCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGTACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAgcg      <  2:1708572/150‑1 (MQ=255)

                                                                                                                                                  caggcaggAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAGCGCAt   >  2:791017/1‑150 (MQ=255)

                                                                                                                                                  caggcaggAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAGCGCAt   >  2:1045517/1‑150 (MQ=255)

                                                                                                                                                  caggcaggAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCAGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAGCGCAt   >  1:1675544/1‑150 (MQ=255)

CGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAGCGCATC  >  NC_000913/644414‑644710

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