BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,856,334	C→A	7.4%	E199* (GAG→TAG) ‡	ydjH ←	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,856,334	0	C	A	7.4%	320.9 / 11.0	94	E199* (GAG→TAG) ‡	ydjH	putative kinase
Reads supporting (aligned to +/- strand): ref base C (48/39); new base A (4/3); total (52/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTACC  >  NC_000913/1856193‑1856481

ccgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCaaaa                                                                                                                                              <  2:861360/149‑1 (MQ=255)

ccgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCaaaa                                                                                                                                              >  1:609840/1‑149 (MQ=255)

 cgGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAAt                                                                                                                                             <  2:84313/149‑1 (MQ=255)

   gCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTa                                                                                                                                           <  2:149036/149‑1 (MQ=255)

    cACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAg                                                                                                                                          <  2:817249/149‑1 (MQ=255)

      ccTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTactc                                                                                                                                        <  4:144397/149‑4 (MQ=255)

       cTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAgg                                                                                                                                         >  3:860313/1‑147 (MQ=255)

       cTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAgg                                                                                                                                         <  4:860313/147‑1 (MQ=255)

       cTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGaa                                                                                                                                       >  3:868036/1‑149 (MQ=255)

        ttCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCTCCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGaaa                                                                                                                                      <  4:164589/149‑1 (MQ=255)

        ttCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCaaaa                                                                                                                                              >  3:786177/1‑141 (MQ=255)

        ttCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCaaaa                                                                                                                                              <  4:786177/141‑1 (MQ=255)

        ttCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGaaa                                                                                                                                      >  2:115398/1‑149 (MQ=255)

         tCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAAc                                                                                                                                     >  2:632047/1‑149 (MQ=255)

                   gTCACCACGCTTGATCAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGAACGCAAAAAAGCAGGCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                                                                           >  2:268945/1‑149 (MQ=255)

                   gTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                                                                           >  4:657462/1‑149 (MQ=255)

                    tCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGAc                                                                                                                          >  3:812997/1‑149 (MQ=255)

                     caccacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACa                                                                                                                         <  4:57354/149‑1 (MQ=255)

                      accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                                                                        >  3:290612/1‑149 (MQ=255)

                      accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                                                                        >  2:448715/1‑149 (MQ=255)

                      accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                                                                        >  2:540439/1‑149 (MQ=255)

                      accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                                                                        >  1:68102/1‑149 (MQ=255)

                      accacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                                                                        >  1:147941/1‑149 (MQ=255)

                        cacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATaa                                                                                                                      >  1:406918/1‑149 (MQ=255)

                        cacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATaa                                                                                                                      >  1:396494/1‑149 (MQ=255)

                           gCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTc                                                                                                                   >  1:645684/1‑149 (MQ=255)

                            cTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATATTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCa                                                                                                                  >  1:232017/1‑149 (MQ=255)

                             ttGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCG‑TTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                                                                                <  3:751107/149‑1 (MQ=255)

                              tGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                                                                                <  2:68102/149‑1 (MQ=255)

                               gATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATg                                                                                                               <  4:15324/149‑1 (MQ=255)

                               gATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATg                                                                                                               >  2:297108/1‑149 (MQ=255)

                                   aaGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                                                                           >  1:102231/1‑149 (MQ=255)

                                   aaGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                                                                           <  4:219417/149‑1 (MQ=255)

                                   aaGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                                                                           <  1:297108/149‑1 (MQ=255)

                                         ccGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAat                                                                                                     >  2:195469/1‑149 (MQ=255)

                                            tCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatc                                                                                                  >  2:872187/1‑149 (MQ=255)

                                              tttACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcat                                                                                                <  2:336288/149‑1 (MQ=255)

                                              tttACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcat                                                                                                <  1:186411/149‑1 (MQ=255)

                                                  cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAg                                                                                            >  2:521924/1‑149 (MQ=255)

                                                   ccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTAGCAAATATCATCCAGc                                                                                           >  1:97518/1‑149 (MQ=255)

                                                      ttttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGtt                                                                                        <  1:632047/149‑1 (MQ=255)

                                                           aTCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcatt                                                                                   <  1:521924/149‑1 (MQ=255)

                                                              accaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCaa                                                                                <  1:580550/149‑1 (MQ=255)

                                                                 accGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAATcg                                                                             >  3:179662/1‑149 (MQ=255)

                                                                 accGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACcg                                                                             >  4:508863/1‑149 (MQ=255)

                                                                 accGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCAAGCGTTTCATTCAACcg                                                                             >  4:566517/1‑149 (MQ=255)

                                                                  ccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACcgc                                                                            >  3:208676/1‑149 (MQ=255)

                                                                    gTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCgg                                                                          >  2:273090/1‑149 (MQ=255)

                                                                    gTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCgg                                                                          <  4:782987/149‑1 (MQ=255)

                                                                      ttttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGtt                                                                        <  1:273090/149‑1 (MQ=255)

                                                                         tACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                                                                     <  2:830070/149‑1 (MQ=255)

                                                                           cgccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATc                                                                   <  1:115398/149‑1 (MQ=255)

                                                                           cgccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATc                                                                   <  3:481670/149‑1 (MQ=255)

                                                                                caaggaaGAAAGCGGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATCTTCCCTCGGCAAAATTAGGTAACAGATAATCGACATAACTCAATGCTTCGCAAATCGCATCTAGCGCTTCATTCAACCGCGGTTTGATCATATc                                                              >  4:69841/6‑149 (MQ=255)

                                                                                     aaGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGcac                                                         >  4:418787/1‑149 (MQ=255)

                                                                                      aGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGcaca                                                        >  2:754406/1‑149 (MQ=255)

                                                                                             aGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatc                                                 >  1:783145/1‑149 (MQ=255)

                                                                                             aGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatc                                                 <  1:444429/149‑1 (MQ=255)

                                                                                               tCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatcat                                               >  3:675270/1‑149 (MQ=255)

                                                                                               tCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatcat                                               <  1:448715/149‑1 (MQ=255)

                                                                                                 aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCt                                             <  4:290612/149‑1 (MQ=255)

                                                                                                 aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCt                                             >  3:187564/1‑149 (MQ=255)

                                                                                                 aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCt                                             >  4:24512/1‑149 (MQ=255)

                                                                                                  gCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTg                                            >  1:628649/1‑149 (MQ=255)

                                                                                                    aaTTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGcc                                          >  4:349886/1‑149 (MQ=255)

                                                                                                           tCCAGTGTCTCTTTCCCGGTGAGTAATATTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCttt                                   <  2:232017/149‑1 (MQ=255)

                                                                                                            ccAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTa                                  <  3:349886/149‑1 (MQ=255)

                                                                                                            ccAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTa                                  <  3:865096/149‑1 (MQ=255)

                                                                                                             cAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAg                                 >  1:259020/1‑149 (MQ=255)

                                                                                                               gtgtCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCt                               >  2:580927/1‑149 (MQ=255)

                                                                                                                tgtCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCtt                              >  1:44602/1‑149 (MQ=255)

                                                                                                                tgtCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCtt                              <  4:812997/149‑1 (MQ=255)

                                                                                                                  tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                                                                     >  2:826376/1‑108 (MQ=255)

                                                                                                                  tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                                                                     <  1:826376/108‑1 (MQ=255)

                                                                                                                  tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGc                            >  4:742353/1‑149 (MQ=255)

                                                                                                                  tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGc                            <  2:783145/149‑1 (MQ=255)

                                                                                                                  tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGc                            <  3:657462/149‑1 (MQ=255)

                                                                                                                  tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGc                            >  3:120896/1‑149 (MQ=255)

                                                                                                                  tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGc                            <  1:647167/149‑1 (MQ=255)

                                                                                                                  tcgctTTCCCGGTGAGTAAGTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGc                            <  1:268945/146‑1 (MQ=255)

                                                                                                                    tctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACATATAATCGACATAACTc                                                                                                                   <  2:228458/60‑1 (MQ=255)

                                                                                                                    tctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTc                                                                                                                   >  1:228458/1‑60 (MQ=255)

                                                                                                                    tctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGt                          >  3:272045/1‑149 (MQ=255)

                                                                                                                         cccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTaaaaa                     <  3:289779/149‑1 (MQ=255)

                                                                                                                          ccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAAt                    >  2:508327/1‑149 (MQ=255)

                                                                                                                             gTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTc                 <  1:540439/149‑1 (MQ=255)

                                                                                                                              tGAATAATTTTGCCTCGACAAAATTAGGAAACAGATAATCGACATAACTCCATGCTTCGAAAATATCATCCAGCGTTTCATTGAACCGCGGTTTGATCATATCGGCACCGCTAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCt                <  2:319284/149‑1 (MQ=255)

                                                                                                                                cctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGt              >  2:257961/3‑149 (MQ=255)

                                                                                                                                  tAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                                                                     <  2:805125/92‑1 (MQ=255)

                                                                                                                                  tAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                                                                     >  1:805125/1‑92 (MQ=255)

                                                                                                                                  tAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTa            >  3:126347/1‑149 (MQ=255)

                                                                                                                                    aTTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAgc          >  3:381140/1‑149 (MQ=255)

                                                                                                                                     ttttGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAgcg         >  4:878122/1‑149 (MQ=255)

                                                                                                                                            tCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTAcc  <  3:173895/149‑1 (MQ=255)

CCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCGCTTTACC  >  NC_000913/1856193‑1856481

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`