BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	644,561	A→T	5.3%	L148Q (CTG→CAG)	rna ←	ribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	644,561	0	A	T	5.3%	295.3 / 4.8	101	L148Q (CTG→CAG)	rna	ribonuclease I
Reads supporting (aligned to +/- strand): ref base A (43/52); new base T (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 = 3.40e-01

Only 100 of 111 total aligned reads displayed.

GCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAGCGCATC  >  NC_000913/644413‑644710

gcggcTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCcag                                                                                                                                                      <  1:437133/150‑1 (MQ=255)

gcggcTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCcag                                                                                                                                                      >  1:817218/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)

               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)

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

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

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

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

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

                              aaGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGAc                                                                                                                        <  1:782423/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                                                                                                                              >  2:87564/1‑137 (MQ=255)

                                     tttCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTgg                                                                                                                              <  1:87564/137‑1 (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:689635/150‑1 (MQ=255)

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

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

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

                                               cGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCgg                                                                                                       <  1:810110/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)

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

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

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

                                                                                     ccGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAacac                                                                 >  2:1844252/1‑150 (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)

                                                                                               cGTCCGGATCAACACCAAAGCAGGCACCGTGGTTGGCATATTCGTAGCGGTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTgg                                                       <  2:638774/150‑1 (MQ=255)

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

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

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

                                                                                                    ggAGCAAACCCAAAGCAGGCACCGTGATTGGCATATTCGTAGCGTTCCAGGCAGGCACGTCCGCCAGCTCCTGGCATGACTTCACTAAGTTTAGCGGCCGCTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGgcgcgc                                                  <  1:49049/150‑1 (MQ=255)

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

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

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

                                                                                                         aaaaCCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTATCCGGCGATGAACACATTCGGCTGGCGCGCGCtgc                                             >  1:999884/1‑148 (MQ=255)

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

                                                                                                          aaaCCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCt                                            >  1:1473218/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)

                                                                                                                           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)

                                                                                                                                tGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAg                      <  1:1319336/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:1707793/150‑1 (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)

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

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

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

                                                                                                                                                    aggcaggAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAGCGCATc  >  2:1519632/1‑150 (MQ=255)

GCGGCTCACTGTTTTACCGTAATTATCCGCAAGGAATTTTCCGGCTTCGCTCTCTTTAATTTCTTGATTCAGGCGTACCATCGTACCGAAGTATGCGTCCGGATCAAAACCAAAGCAGGCACCGTGTTTGGCATATTCGTAGCGTTCCAGGCAGGAACGTCCGCCAGCTCCTGGCATGACTTCACTTAGTTTAGCGGCCGTTTCCAGTGATAATCCGGTTTCCGGCGATGAACACATTCGGCTGGCGCGCGCTTCTGGTAGATTCGGGATTGGGCGAGTAGCGCAACCGAAGCGCATC  >  NC_000913/644413‑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: `‑`