BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,702,505	C→G	5.8%	G91G (GGC→GGG)	add →	adenosine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,702,505	0	C	G	5.8%	473.1 / 14.0	138	G91G (GGC→GGG)	add	adenosine deaminase
Reads supporting (aligned to +/- strand): ref base C (60/70); new base G (6/2); total (66/72)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

Only 100 of 142 total aligned reads displayed.

ACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGA  >  NC_000913/1702356‑1702653

acacTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGc                                                                                                                                                      <  1:2063257/150‑1 (MQ=255)

acacTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGc                                                                                                                                                      <  1:1148557/150‑1 (MQ=255)

   cTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTg                                                                                                                                                   <  1:1527557/150‑1 (MQ=255)

     gATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCa                                                                                                                                                 >  2:358329/1‑150 (MQ=255)

      aTTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCAc                                                                                                                                                <  2:1275763/150‑1 (MQ=255)

        tCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGGCCGTGACGGGCTGCatct                                                                                                                                              <  2:132377/150‑4 (MQ=255)

        tCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTa                                                                                                                                              <  1:2433062/150‑1 (MQ=255)

         ccccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTAt                                                                                                                                             >  1:2004991/1‑150 (MQ=255)

          cccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATg                                                                                                                                            >  1:515131/1‑150 (MQ=255)

          cccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATg                                                                                                                                            >  1:1597083/1‑150 (MQ=255)

          cccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATg                                                                                                                                            >  2:1181068/1‑150 (MQ=255)

           ccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGt                                                                                                                                           >  1:1520588/1‑150 (MQ=255)

           ccACGTTCAGGTCATTGCCAACGAACCCGATCTGCTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGt                                                                                                                                           >  2:2083779/1‑150 (MQ=255)

            cACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTc                                                                                                                                          >  1:581362/1‑150 (MQ=255)

              cGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGa                                                                                                                                        <  1:54017/150‑1 (MQ=255)

                ttCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGc                                                                                                                                      <  1:35515/150‑1 (MQ=255)

                  cAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTg                                                                                                                                    <  2:628540/150‑1 (MQ=255)

                   aGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGc                                                                                                                                   <  1:620975/150‑1 (MQ=255)

                    ggTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAACCTTGCCTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGGCACGGCCTGCACTATGTCGAGCTGCg                                                                                                                                  <  2:1276882/150‑1 (MQ=255)

                     gTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGt                                                                                                                                 <  2:15462/150‑1 (MQ=255)

                      tCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGtt                                                                                                                                >  1:1234917/1‑150 (MQ=255)

                       cATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGttt                                                                                                                               <  1:672749/150‑1 (MQ=255)

                       cATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGttt                                                                                                                               >  2:1567678/1‑150 (MQ=255)

                        aTTTCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGtttt                                                                                                                              <  2:2320896/150‑1 (MQ=255)

                        aTTGCCAACGAACCCGATCTGGTGAGCTTTCTGCCCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGtttt                                                                                                                              <  2:1788729/150‑1 (MQ=255)

                          tGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTTGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTc                                                                                                                            <  2:1870090/150‑1 (MQ=255)

                                  aaCCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCt                                                                                                                    >  2:1440142/1‑150 (MQ=255)

                                  aaCCCGATCTGGTGAGCTTTATGAACACACTTTACTGGGGCGTTAAAGTTCTCGCCTCTATTGATGCCTGTCGCCGCGTTGCATTTGAAAACATTGAATATGCAGACCGGCAAGTCCTGCACTATGGCTAGCTTCTTTTTTCACCATGCt                                                                                                                    <  2:2424566/150‑1 (MQ=255)

                                       gATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatg                                                                                                               <  1:817141/150‑1 (MQ=255)

                                        aTCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatgg                                                                                                              <  1:1435482/150‑1 (MQ=255)

                                        aTCTGGTGAGCTTTCTGACCAAACCTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatgg                                                                                                              <  1:325347/150‑1 (MQ=255)

                                            ggTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaat                                                                                                          >  2:2123267/1‑150 (MQ=255)

                                              tGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatgg                                                                                                        <  1:2027090/150‑1 (MQ=255)

                                              tGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatgg                                                                                                        <  1:2083779/150‑1 (MQ=255)

                                               gAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatggc                                                                                                       <  2:1355561/150‑1 (MQ=255)

                                               gAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatggc                                                                                                       >  1:1311432/1‑150 (MQ=255)

                                                     tCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCGGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCa                                                                                                 <  2:1218870/150‑1 (MQ=255)

                                                        gACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCATGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATCGCAATGGCACATCAGCt                                                                                              >  2:2193888/1‑150 (MQ=255)

                                                         aCCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTg                                                                                             <  1:1181068/150‑1 (MQ=255)

                                                          ccAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGc                                                                                            <  1:946904/150‑1 (MQ=255)

                                                           cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGcc                                                                                           <  1:1309656/150‑1 (MQ=255)

                                                           cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCACTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGcc                                                                                           <  1:2340563/150‑1 (MQ=255)

                                                            aaaCTTGACTGGGGCGGTAAAGTTCTCGCCTCTCTGGATGCCTGTCGCCGCGTTGCATTTGAAAACATTGAAGATGCAGCCCGTCACGCCCTGCACTATGTCGAGCTGCGGTTTTCACGAGGCTCCATGGCAATGGCCCATCAGCTGCCt                                                                                          >  2:2168844/1‑150 (MQ=255)

                                                                    cTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCggg                                                                                  <  1:444839/150‑1 (MQ=255)

                                                                    cTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCggg                                                                                  <  2:863311/150‑1 (MQ=255)

                                                                    cTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCggg                                                                                  <  1:1648755/150‑1 (MQ=255)

                                                                    cTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCggg                                                                                  >  1:1847631/1‑150 (MQ=255)

                                                                         gCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGtgttg                                                                             >  1:1694386/1‑150 (MQ=255)

                                                                           gTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTc                                                                           <  1:1688815/150‑1 (MQ=255)

                                                                             tAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGa                                                                         >  1:1848157/1‑150 (MQ=255)

                                                                             tAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGa                                                                         <  2:474272/150‑1 (MQ=255)

                                                                             tAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGa                                                                         >  1:1244196/1‑150 (MQ=255)

                                                                                  ttCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTTCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCgg                                                                    <  2:837624/150‑1 (MQ=255)

                                                                                   tctcGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGt                                                                   <  1:2425249/150‑1 (MQ=255)

                                                                                      cGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGGATTATGTATAGCTGCGTTTTTCACGAGGCTACATGTCAATGTCACTTTTGCTGCATTTATCGGGTTTTGTCGTAGCGTTGAt                                                                <  2:841018/150‑1 (MQ=255)

                                                                                       gCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATc                                                               <  2:1220155/150‑1 (MQ=255)

                                                                                           ctctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATg                                                           <  2:1877805/150‑1 (MQ=255)

                                                                                             ctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGc                                                         >  1:1259013/1‑150 (MQ=255)

                                                                                               tGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGt                                                       <  2:1171560/150‑1 (MQ=255)

                                                                                                 aTGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTAc                                                     >  2:1866537/1‑150 (MQ=255)

                                                                                                   gCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGt                                                   <  2:309118/150‑1 (MQ=255)

                                                                                                       gTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAg                                               >  1:1312926/1‑150 (MQ=255)

                                                                                                       gTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAg                                               >  1:142599/1‑150 (MQ=255)

                                                                                                           ccgcGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTg                                           >  2:709845/1‑150 (MQ=255)

                                                                                                             gcgTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGcc                                         <  2:879792/150‑1 (MQ=255)

                                                                                                              cgTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCg                                        <  2:6085/150‑1 (MQ=255)

                                                                                                              cgTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCg                                        <  2:150360/150‑1 (MQ=255)

                                                                                                               gTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGc                                       >  2:2426253/1‑150 (MQ=255)

                                                                                                               gTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGCGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGc                                       >  2:2418084/1‑150 (MQ=255)

                                                                                                                 ggCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTACCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCAc                                     >  1:834977/1‑150 (MQ=255)

                                                                                                                    aTTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCCAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGCTTGCCGCACCtt                                  >  2:867702/1‑150 (MQ=255)

                                                                                                                     tttGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCttt                                 <  1:2484727/150‑1 (MQ=255)

                                                                                                                       tGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTgg                               <  2:69756/150‑1 (MQ=255)

                                                                                                                       tGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTgg                               <  2:253270/150‑1 (MQ=255)

                                                                                                                        gAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGTAGGTTGCCGGAGCTTTGgt                              <  2:515131/150‑1 (MQ=255)

                                                                                                                         aaaaCATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGgtg                             >  1:1601998/1‑150 (MQ=255)

                                                                                                                         aaaaCATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGgtg                             >  1:1646641/1‑150 (MQ=255)

                                                                                                                          aaaCATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGgtgt                            <  2:1478411/150‑1 (MQ=255)

                                                                                                                            aCATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGc                          >  2:1689177/1‑150 (MQ=255)

                                                                                                                               gacatagtGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAgg                       >  1:29613/8‑150 (MQ=255)

                                                                                                                                  aaGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGa                    >  2:596459/1‑150 (MQ=255)

                                                                                                                                   tagtGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGaa                   >  1:276029/4‑150 (MQ=255)

                                                                                                                                    gATGCAGCCCGTCACGGCCTGCTCTATGTCGAGCTACGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAg                  <  1:2460409/150‑1 (MQ=255)

                                                                                                                                    gATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAg                  >  2:620652/1‑150 (MQ=255)

                                                                                                                                    gATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAg                  >  1:2450037/1‑150 (MQ=255)

                                                                                                                                    gATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAg                  <  1:620652/150‑1 (MQ=255)

                                                                                                                                    agtGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAg                  >  2:1254981/3‑150 (MQ=255)

                                                                                                                                     gtGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGc                 >  2:788339/2‑150 (MQ=255)

                                                                                                                                     gtGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGc                 >  2:881382/2‑150 (MQ=255)

                                                                                                                                     aTGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGc                 <  1:2013223/150‑1 (MQ=255)

                                                                                                                                     aTGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGc                 <  2:581362/150‑1 (MQ=255)

                                                                                                                                     aTGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCCTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGc                 >  1:1479534/1‑150 (MQ=255)

                                                                                                                                      tGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCt                >  1:1753119/1‑150 (MQ=255)

                                                                                                                                      tGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCt                >  2:1744887/1‑150 (MQ=255)

                                                                                                                                      tGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCt                >  2:1520285/1‑150 (MQ=255)

                                                                                                                                           cccGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCg           <  1:561655/150‑1 (MQ=255)

                                                                                                                                            ccGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCgg          >  2:1187686/1‑150 (MQ=255)

                                                                                                                                              gTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCa        >  2:1385972/1‑150 (MQ=255)

                                                                                                                                                 aCGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTa     <  1:2311135/150‑1 (MQ=255)

                                                                                                                                                  cGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATAAt    >  2:105577/1‑150 (MQ=255)

ACACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGA  >  NC_000913/1702356‑1702653

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